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CHAPTER | 


Mathematical methods and 
differential equations 


|. Introduction 


In the first two chapters of this course we need to lay the groundwork for 
the more advanced material that comes later, so here the emphasis is on 
consolidation of material that you may have seen before. We also need to 
be sure that you are familiar with the notation that we shall use later. 


In this chapter we concentrate on areas that are essential components of an 
applied mathematician’s tool-kit: functions, complex numbers and differen- 
tial equations. In addition to reminding you of some useful techniques, we 
also provide you with a number of exercises that are preparatory to material 
that will be discussed in later chapters. 


1.2 Some initial definitions and 
results 


1.2.1 Functions 


A function f is a mapping from a set A known as the domain, to a set B 
known as the codomain, in which every element of A maps to a unique 
element of B. The notation 


y=f(x) for tEA 


is used to encapsulate such a relationship, where x is an element of A, y is 
an element of B, and the symbol € is an abbreviation for ‘is a member of’. 


Often A and B are sets of real numbers, in which case f is said to be a 
real function. For example, f(x) = x? for x € R, where R is the set of real 
numbers, is a real function. The set R is often referred to as the real line. 
A real function f is even if f(x) = f(—x), and odd if f(z) = —f(—a). For 
example, y = x” is even, while y = 2° is odd. 


A function g is an inverse of the function f if f(g(z)) = x, and it is common 
practice to denote an inverse of a function f by f~'. For example, if 
f(x) = Vz for x > 0 (where \/z denotes the positive square root of x), then 


f-'(z) = «*, because V22 = x. 
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Later we shall meet functions for which the domain is a set of real numbers 
and the codomain is a set of complex numbers. Such functions are known 
as complex functions of a real variable. ~~ 


We shall also meet functions for which the domain is itself a set of functions. 
Such functions are often known as operators or transforms. This may 
at first seem a rather strange concept, but you have in fact met the idea 
before. Whenever you differentiate a function, you are actually applying 
the differential operator to a function to obtain another function, namely 
its derivative. So differentiation is in fact an operator that maps functions 
to their derivatives. We shall not discuss operators in any great depth, but 
the concept will occur again when we introduce the Fourier transform later 
in the course. 


The domains and codomains of specific functions are often defined implicitly. 
For example, for the real function 


1 
1" -# 


it is implicitly assumed that the domain is the set —1 < x < 1 since this is 
the largest set of real numbers on which the definition is sensible. Strictly 
speaking, we should distinguish between a function f and a function value 


f(x), but generally we shall make no such distinction and, for example, refer 
1 


V1 — 22 

Functions themselves are sometimes defined implicitly. For example, the 
equation (of a circle) 2* + y? = 1 determines two functions, y = V1 — x? 
and y = —V1—<2?, and we require more information in order to discover 


which of these functions is intended. In either case the domain would be the 
set -—L<a<l. 


to the function y = 


Smooth functions 


Often we shall need a function to be ‘sufficiently smooth’. The minimum 
requirement for this is that the function should be continuous at the point 
in question. In intuitive terms, a function is continuous if its graph has no 
sudden breaks, but we can be a little more precise. A function f is said 
to be continuous at a point x = a if Tim f(x) = f(a). In other words, as x 


approaches a from above (x > a) or below (x < a), the function value f(x) 
approaches the value f(a). 


In order for the function to be ‘smooth’, we also require that it should be 


differentiable at the point in question. This means that for a function f It can be shown that a 
f(x) — f(a) function that is 
to be smooth at the point x = a, we require that the limit lim —~————— differentiable at a point 


La a ae s 
-.. <n a must be continuous there. 
exists. Again in intuitive terms, a smooth function is one for which the 


graph has a tangent at the point in question. 


d 
The derivative of a function y = f(x) is often denoted by ce although we 


dx 
d 
also denote the derived function by f’, so f’(x) = =: 
2 
The second derivative is usually written as f(x) or aa? while the higher 
x 


a” 
derivatives are written as f(x) or a bor Bs 
= 


When we require a function to be ‘sufficiently smooth’, we generally mean 
that the function should be differentiable a sufficient number of times at the 
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relevant points of its domain. Sometimes this may imply that the function 
should be differentiable ‘an infinite number of times’. (For example, for the 
sine function f(x) = sinx, we have f’(x) = cosz and f"(x) = —sinz, thus 
the function is differentiable as many times as we please.) 


A polynomial is a function of the form 
p(x) =ap + a,x + ion eee” 


and the polynomial is of degree n if a, # 0. Clearly a polynomial is ‘sufh- 
ciently smooth’ at any point on the real line because it can be differentiated 
as many times as we please. 


zs 

A rational function is a function of the form ae . where p and gq are 
q(x 

polynomials, and such a function is ‘sufficiently smooth’ at every point on 


the real line except the roots of the equation q(x) = 0. At those points the 
rational function is discontinuous. 


1.2.2 Series 


Convergence of series 


Example I.| 
You may have already met the geometric series 
Sn(a)=1+a+a?+---+a", 


where n is an integer and a is some specified number. Such a sum can be 
written in an alternative form if we multiply both sides of the above equation 
by a to obtain 


aS(a) =a +07 ->- be", 
then subtract one equation from the other and rearrange: 
1—a” 


aa) = , (Lad) 


l—-a 


For large values of n, and for any choice of the number a in the interval 
Tr 


1 
—1 < a.< 1, the term = is small, therefore S,,(a) ~ ase More pre- 
—a 


cisely, 


lim S,,(a) : 


N00 ae ee a 


provided that -l<a<l. Hf 


(1.2) 


It is this final statement that leads us to the definition of the convergence 
of an infinite series. We begin by considering a finite sum 
Sq Ag+ Ar + Ao + +p Anal, 1 
als 
which may be conveniently written in the form S;, = ‘Sy A;. We then say 
ore) k=() 
that the infinite series vi Ary is convergent to the sum S if lim S, = 8S. 
Ti CO 
r= 
If the limit lim S,, does not exist, then we say that the series is divergent. 
T= OO 
n—l 
A finite sum such as > A, is often known as a partial sum of the infinite 
series. k= 


The symbol ~ is read as 
‘is approximately equal to’. 
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CO 
For the previous example we could then write bs q* = 
ie=() 
and S;,(a) =1+a+a*+---+a”"! would be the nth partial sum of the 
infinite series. 


if—-l<a<l, 
—a 


Sometimes it is convenient to write an infinite series in a more informal 
manner, and for example we may write 


CO 
S| Ag = Ao + Ai + Ag t++++ Ante, 
k=0 
where the final set of dots indicates that the series continues indefinitely. 


Exercise |.1 
OO 


Show that for —-l1<a<1, R,(a)= >. a’ can be written in the form 
q” k=n 


Rate of convergence 


In some circumstances it is not sufficient simply to know that an infinite 
series converges; we may also need to know how rapidly the series con- 


if 


verges. In the previous example we saw that the infinite sum was 


1 
—1<a< 1, which means that S,(a) — ; is small when n is large. It is 
—a 


instructive to rearrange equation (1.2) into the form 


1 
oo = Paes + ot ale), 


where R,(a) = : 


the remainder term tends to zero determines the rate of convergence of 
the series, and, loosely speaking, in this case we can see that R,,(a) behaves 
like a”. However, intuitive statements like this can be made more precise 
by introducing some additional notation. 


a” 
is known as the remainder term. The rate at which 
—a@a 


The O notation 


The O notation is a convenient device for specifying an order of magnitude. 
For example, sin x is of the same order of magnitude as x when z is small. To 
make such a statement more precise, we define f(x) = O(g()) if there exists 
a constant M such that |f(2)/g(x)| < M for all sufficiently small values of z. 
Thus, for example, it is possible to show that sinz = O(x) and 22 + 327 = 
O(a) if x is small. In the context of our previous example, we could write 
R,(a) = O(a”) for a sufficiently small, to indicate that |R,(a)/a”| < M for 
some constant M > 1. So if we are concerned only with orders of magnitude, 
we could write 


1 
—— =l+a+a°+---+a"'4+O(a"). 


i 
More generally, for a function f(x) = g(x) + anz" + Qn4i 2"! +--+ with 
constant coefficients an, @n41,..-, we can write f(x) = g(x) + O(a”) for suf- 


ficiently small values of x. 
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Example 1|.2 


Given f(x) = O(x?) and g(x) = O(a?) when z is small, establish the orders 
of magnitude of f(x)g(x) and f(x) + g(2). 


Solution 


We know that |f(x)| < k,\z|? and |g(x)| < kelz|?, for some constants ky; 
and ky. Therefore |f(x)g(x)| < kike|x|°, so f(x)g(x) = O(a”). 


Also, we have 
f(x) + 9(x)| < |F(z)| + lg(@)| < Fala? + halal? < 2k la)? 
if x is sufficiently small. Hence f(x) + 9(r) = O(z?). 


Taylor series 


It can be shown that a function f that is sufficiently smooth near a point 
x =a can be represented as a series 


= (k) (a 
2) = S- ag(a —a)* where a; = J a ) (4:3) 
kaa 


Such a series is known as the Taylor series (or Taylor expansion) for f 
near the point x = a, and it is generally convergent for x in an interval 
a—p<x2<a+p (which can be written as |x — a| < p) for some positive 
constant p. The largest possible choice for p is known as the radius of 
convergence of the series. 


A sequence of numbers such as aj, d2,a3,... is often abbreviated to {ax}, 
and in this case these numbers are known as the coefficients of the series. 


Example I.3 


Write down the Taylor series for the function f(z) = fc oS 0. 
Solution 
For this function, f*) (x) = Bees Sef“) (0) = (—1)*R}, and the-co- 
(1+a)etl i 
efficients of the Taylor series are ax, = (—1)*. Thus 
1 ee) 
= 5 (-1)*c*. & 
Lax —. 


Exercise 1.2 


Find the largest interval |x| < p on which the series 
CO 
l—-g?tet—7o +... = lye 
k=0 
is convergent. 


Here we have used the 

triangle inequality 
la + b] < |a| + [6] 

(see Exercise 1.14). 


Note that 
= 7 
and 


f(a) = f(a). 
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1.2.3 Some elementary functions and their 
derivatives 
You will probably have previously encountered most of the following func- 


tions and their Taylor series near x = 0. (We include the geometric series 
again for completeness.) We also give the interval of convergence for each 


series. 
1 CO 
: aaa a Soe ad, 
k=0 
2 n ern 
nL. Liat Soh = an ey 
expr =e oi a ser es = a for x ER, 
k=O) 
3 i O° k p2k+1 
ee (—1)"z (=i y 
eee re te ee a eee erent serene ie R 
me = 3) ee 44° (2k +1)! Oe ee 
2 Wnan ie OS 
¥ (—1)"x (—1)"a 
mesa ae for x € R, 


ze” - —e" 

In(l-—z)=- (2+ See Z es) nin) for —l<@ < 1. 
k=21 

For this course, the natural logarithm of x will be denoted by Inz, although Inz is the logarithm to 

other authors write log, x or just log z. For x > 0 the function In z is defined _ base e. 

to be the inverse of the exponential function, so exp(Inz) = x for x > 0. The 

radius of convergence of the Taylor series for In(1 — x) near x = 0 is 1, but 

the series diverges at x = 1 and converges at x = —1. 


Many other series can be obtained by simple manipulations of the above 
series, as the following example shows. 


Example |.4 


Write down the Taylor series for the following functions near t = 0. 


Cee at _ 


Solution 


(a) Replacing x by —2t in the series for In(1 — x), we obtain 


in{i +2 = = = 
k=] 


(b) It can be shown that we are justified in differentiating the geometric 
series term by term to obtain 


1 d 1 
( ) = 1420 +808 + pratt 


(1-2)? ~ de \1—a - 
=e 
| 


Then replacing x by 3¢ gives 
1 


see 27H? +--+ +n(3t)” 1 +---=S °k(3t)*". OE 


feed 


Exercise 1.3 


Show that sin x In(1 + 2) = x? + O(2°). 
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The remaining trigonometric functions can be defined in terms of sine and 


cosine: 
sin x Coa x 1 
tanx = ——, cotr= — = , secxr = ——, cosecr = — 
COS & sinzx tanzx COS sin © 
The hyperbolic functions 
The hyperbolic sine and cosine functions, sinh and cosh, are defined by Sinh and cosh are 
er J a pant! eer vines and 
Set ee eee ee ee ee re ee Cosn , respectively. 
a it. @ (2n +1)! 
OC gale’ 
= ———, convergent for x € R 
2k + 1)’ : 
k=0 
ef + et op en 
Cosh a se eee ee ee she 
2 2! (2n)! 
ad 2k 
2. , convergent for x € R. 
= (Ze 


The remaining hyperbolic functions are defined by 


sinh x cosh x 1 
tanh xz = ee coke = ee a Tanh, coth, sech and 
aa sinha tanh x cosech are pronounced 
‘tansh’, ‘coth’, ‘sesh’ and 
sech 2 == , cosechz = — * ee 
cosh x sinh x co-sesh’, respectively. 


The graphs of sinh, cosh and tanh are shown in Figure 1.1. 


Figure |.!_ ‘The graphs of sinh, cosh and tanh 


Exercise 1.4 


From the definitions of sinh and cosh, show that cosh? x — sinh? x = 1. 


Example 1.5 
1 
Show that = (cosh! x) ae 
dx x2 — 1 


Solution 
If w = cosh~! 2, then coshw = x. Differentiating with respect to z, using 
dw 
the chain rule, we have sinh w > pa 1. Thus 
2 


dw 1 1 


1 


12 


The Gaussian function 


The Gaussian function 


y = exp[—(x — 1)” /20°], 


where p and o are constants, arises in statistics. 
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Its graph is shown in 


Figure 1.2 for various values of jz and o. 


= 


—2 -l1 0 


Figure 1.2 Graphs of the Gaussian function 


identities and derivatives 


Various identities for the trigonometric and hyperbolic functions arise quite 
often, and you may find the following lists useful. 


y 2 


cos’ x+sin“ x = 1 

sin(x + y) =sinxcosy + coszsiny 

cos(x + y) = cosxcosy — sinxzsiny 

tan.z + tan 

tan(z + y) = a 
1 — tans tan y 


Function f(z) | Derivative f’(x) 


sin x 
COS x 
tan x 
COt Tr 


COS X 
—sinx 
sec? x 

— cosec? x 
sec x tan x 


— cosec £ cot x 


SCC L 
COSECC ZL 


1.2.4 Integration 


A primitive of a given function f is a function F' such that F” 


the notation 


cosh? x — sinh? x = 1 

sinh(« + y) = sinh xcoshy + coshz sinh y 

cosh(x + y) = coshzcoshy + sinh z sinh y 
tanh x + tanh y 

1+ tanh z tanh y 


Function f(x) | Derivative f’(x) 


tanh(z + y) = 


Ina L/x 

sinh x cosh x 

cosh x sinh x 

sain & sech? 2 

coth x — cosech? x 
sech x —sechztanhz 


— cosech x coth x 
1 


Vaz +1 


cosech x 


sinh7! x 


cosh7! x 


tanh! x 


= 7 vRhe 
term ‘indefinite integral’ is often used as an alternative to ‘primitive’, and 


f(x) dx is used to indicate an arbitrary indefinite integral 


of f. Thus indefinite integration arises as the ‘inverse of differentiation’. 
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|.2 Some initial definitions and results 


On the other hand, definite integration originates from the need to eval- 
uate a limit of a sum, for example arising from the area under a curve. 


These two notions of integration are linked by the fundamental theorem of 
calculus (see below), but it is important to appreciate that they are quite 
distinct concepts. 


In an integral of the form / f(x) dz, the function f(x) being integrated is 


called the integrand, and the variable x is called the integration variable. 
Another common way of expressing such integrals is to write the integrand 


laa, 1. / dx f(x). Both notations are used in this course. 


Definite integration is based on the concept of a limit of a sum. Suppose 
that we are given a function f(a) defined on the interval a < x < b. We first 
divide the interval into N subintervals, with end points 


= te < 2 < te < +s Cay SS 
N 
Then we form the sum S = > f (tr)(@p — Lp_1). Often it is convenient to 
a | 
choose subintervals of equal length, which we might denote as dx, and then 
N 


this sum becomes S = =. f (a,) 6x, which is often abbreviated to )> f(a) da. 
r=k 

Finally, we take the limit as the number of intervals increases, and the 

maximum length of the subintervals shrinks to zero, giving 


b N 
/ f(x) dx = im f (&r)(@r — Lr—1). 
- ea | 


It is important to appreciate that this is how the integral is defined, but to 
evaluate an integral we generally try to employ the fundamental theorem 
of calculus. This tells us that if we are able to find a primitive F(x) such 
that F’(x) = f(x), then 


b 
/ f(x) dx = F(b) — F(a). 


Example 1.6 


Find the area bounded by the curve y = Ina, the x-axis, and the lines x = 2 
and x = 3 (see Figure 1.3). 


Solution 


Figure 1.3 


If we approximate the required area A by rectangles, as shown in the figure, 
then A ~ 5° In(x) da. Now we take the limit as the number of intervals 
increases, and the length of the subintervals shrinks to zero, to obtain 


r= Fi In(a) de. 


In an elementary 
introduction to calculus, 
such a sum S is usually 
identified with an estimate 
for the area under the 
graph of f(x) from x =a 
2 = &: 
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To evaluate this integral we may use integration by parts, which gives 


3 3 
A= ii In(xyde = | ln ]; — / dx Try differentiating x ln x if 
2 2 you cannot see why this 


- E Ine — r]; works. 


= 3In3—2In2—1=In(27/4)-1. @ 


Example |.7 


Suppose that we are given a wire that does not vary in any way along its 
length: i.e. its cross-section stays the same, and the material from which it 
is made remains fixed. We say that such a wire is uniform, and we expect 
the mass of any segment of such a wire to be proportional to the length of 
the segment: a piece twice as long would weigh twice as much. 


The ratio of the total mass of the wire to its total length is known as its 
linear density, and if we denote this quantity by p, then the mass M of 
a length L of wire is pL. If M is measured in kilograms (kg) and L in 
metres (m), then the units of p are kilograms per metre (kg m7‘). 


Now suppose that the wire is non-uniform: how should we extend our def- 
inition of linear density? Imagine that the wire is placed along the z-axis, 
and concentrate attention on a segment of wire of length dx containing a 
particular point x. If the mass of this segment is dm, then we define the 


om 
F d it i? t b li = * 
inear density p(x) at x to be Pear Ox 


In order to be more specific, let us consider a wire which is 1 metre long, for 
which p(x) = 2x +1, so that the density increases steadily from the value 
1kgm7! at the origin to 3kgm7! at the other end of the wire. The total 
mass of the wire can then be expressed as a definite integral. 


First we divide the wire into N subintervals, with end points 
= fp Sc << tom Se aris SI ae. 


If the intervals are small, then p(z,) is a good approximation to the linear 
density at all points of the interval 7,_1 <x < x,, so the mass of this short 
segment of wire is 6m, & p(x,)(@, — 1). Adding all of these small masses, 


we obtain an approximation for the total mass of the wire: 
If we choose equal 


n nr 
int ls, th 
as = S¢ dm, fs ¥ pla) (@r = subintervals, then this 
fil pal 


approximation can be 
written as M ~ S> p(x) ox. 


Tr 
= a MF, te, «. WHEE HL, = By — Vy—1. 
r=] 
This approximation becomes more accurate as we increase the number of 
intervals and the length of each subinterval shrinks to zero. Thus in the 
limit we have 


m= fi o(a)ar. 


To evaluate this integral for the particular function under consideration, we 
employ the fundamental theorem of calculus: 


= = xL = 17 me 
M= fl pe)ae= f (20+ 1)ae =| eee 2 


i.e. the wire has a mass of 2kg. IH 


Example 1.8 


Calculate the volume of a sphere of radius r. 


1.2 Some initial definitions and results 


Solution 


First divide the sphere into a number of vertical slices of thickness 6x” (rather 
like slicing an onion), then approximate a slice at position x by a cylinder 
of depth dz and radius y = Vr? — x? (see the shaded region in Figure 1.4). 


Figure 1.4 A sphere divided into elemental ‘cylindrical slices’ 


The volume of such a cylinder is ty?6a = m(r? — x?)dx, so the sum 
S = ¥o ay*é2 =a (re? = 2 ee 


provides an approximation to the required volume. In the limit as the num- 


ber of subdivisions increases and 6x tends to zero, the sum becomes a definite 
ie 


integral and we obtain the required volume, V = 7 / (r? — 2”) dx. Now we 
ae 
may use the fundamental theorem of calculus to obtain 


‘ 
= «| (r? — 27) da =a kaze: -- 0°] = aar®. ® 


Kr 


Integrals over an infinite range 


The integrals discussed so far have been defined over finite intervals, but we 
shall also need integrals defined over infinite ranges. 


Example | .9 


Suppose that we are told that an object is at the origin at time t = 0, and 
that it moves along the z-axis in such a way that its velocity at time t is 
given by v = cost. What is the maximum displacement of the object from 
the origin? 


Solution 
dx 
We know that 3 = t, so x(t) is a primitive (or indefinite integral) 


of v. Thus x(t) = | costdt = sint+c, where c is a constant. We choose c so 


that x(t) satisfies the initial condition «(0) = 0, thus c= 0. So z(t) = sint, 
and it is clear that the maximum value of x(t) is 1. 


Note that this result can also be found using definite integration, by setting 
the lower limit equal to the initial condition x(0) = 0. At any given time T, 


2 
we have «(T) = / costdt=sinT. 
0 
In the previous example the object moved backwards and forwards about 


the origin (consider the graph of x(t) = sint), but what happens if the object 
moves always away from the origin but at a slower and slower pace? 
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Example 1.10 


Suppose that we are told that an object is at the origin at time ¢ = 0, and 
that it moves along the x-axis in such a way that its velocity at time ¢ is 
given by v = 1/(1+1)*. What is the maximum displacement of the object 
from the origin? 


Solution 


As in the previous example, we can use either definite or indefinite inte- 
gration. We choose definite integration here, but you might like to repeat 
the calculation using indefinite integration, using the initial condition to 
calculate the value of the arbitrary constant. 


At any given time 7’, we have 


= z 
1 1 1 
ae) | (+t) | eri Pee 
As time T increases, the displacement x(T7) increases from zero, but we can 


see from Figure 1.5 that it never actually achieves a maximum value. Note 
also that the slope of the graph, v(T’), decreases as J’ increases. 


x(T) =1-1/1+T) 


Figure 1.5 The graph of 2(T) =1-—1/(1+T) 


1 
ee . approaches the value 1 from 
below (but more and more slowly as the speed of the object decreases). So 


sees | Re ive 
lim / —___— qt ) = 1, and we normally write / ——dt=1. EH 
T—0oo 0 (1 + 3 0 (1 -_ 2 


Thus we are led to our general definition of an integral over an infinite range. 


As T increases, the value of x(T) = 1 — 


x 
If the limit jim ( / vit) it) exists and equals a finite value L, then we 
= a 


CO 


CO 
say that the infinite integral is convergent and write / jy a= £. 
a 


If the limit does not exist, then we say that the integral is divergent. 


CO 
In the context of the previous examples, we would not expect / u(t) dt to 


a 
exist if the object fails to slow down, but you might reasonably expect that 
the condition Jim v(t) = 0 would be sufficient to ensure the convergence of 
7 OO) 


CoO 
/ v(t) dt. However, this is not the case, as the following example shows. 
a 


Example I.I1 


Suppose that we are told that an object is at the origin at time ¢t = 0, and 
that it moves along the z-axis in such a way that its velocity at time f is 
given by v=1/(1+¢t). What is the maximum displacement of the object 
from the origin? 


|.2 Some initial definitions and results 


Solution 


As in the previous example, the object moves more and more slowly as it 
moves away from the origin. However, we now have 


fa 
1 
ae a —  dt=In(1+T), 
r= f gat= i +7) 
and «(7’) increases without bound as T’ increases. So here we have a cir- 
cumstance where the object is moving ever more slowly, but eventually it 
will move beyond any fixed point that we care to mention. In this case 


CO 

i 

/ co dt does not exist, and the integral is divergent. 
0 


Sometimes it is possible to establish that an integral is convergent without 
actually calculating the value of the limit. Suppose that u(t) > 0 and that 


o(r)= f o)ar<K 


for some fixed constant K and all JT’ > a. In this case, differentiating with 
respect to T' gives x/(T) = v(T) > 0, therefore x(T) cannot decrease with 
time JT’, but we also know that it cannot exceed the value K. So it must 
approach some fixed value (which is less than or equal to K). 


Example 1.12 

o.) et 
Show that the integral / = dt is convergent. 
1 


Solution 


We use the notation introduced in the previous examples. In this case, 
= 
e 

ai) = —- > Ot 2 1, 80 


ony= [ Sa 


t 


-t 
e 
increases with T. But = ee tt; so 


fy et 7 
ce —dt< a [—e"| i oe 
ee 3 1 ; 


So «(T) is non-decreasing and bounded above, so we can conclude that 
ss sa 
jim x(T) exists and therefore the integral / ‘i= dt is convergent. I 
= O©) i 


The Gaussian integral 


The following integral identity is very useful in many contexts where Gaus- 
sian functions appear: 


/ exp(—ax’) dx = /z where a > 0. (1.4) 


This important relation is obtained by a clever trick, where the square of 
the integral is evaluated using polar coordinates. It will be derived in the 
next chapter in the section on multiple integrals. 
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Exercise 1.5 


(a) The following result is often useful when estimating the size of functions that 
can be expressed as integrals. 


If f(x) is a function such that |f(x)| <M fora<a <b, 


/ " fla) da 


Use this result to show that Inz < x for x > 1. 


then < M(b—-a). 


(b) Using the method of integration by parts, or otherwise, evaluate the following. 


va T 
l l 
(i) / dx (ii) / — da 
1 L i 


(c) Show the following. 


= 
(i) / — dx is divergent. 
1 


. a ee 
(ii) —,- dz is convergent. 
1 4 


Exercise |.6 
The following result is often useful when comparing the size of integrals. 


If f(x) and g(x) are two functions such that f(x) < g(a) fora< 2 <b, 


b b 
then / f(cjde= / g(a) de. 
Use this result to show that / exp(—«*) dx is convergent. 
0 


Hint First split the integral as 


ie exp(—«”) dx = [ exp(—2?) dx + F: exp(—2?) ax. 


1.2.5 Partial fractions 


The following process of adding fractions is probably familiar to you: 
: 3. Be!) alee 5(x — 1) 


e-$ #43 Get — o29eso 
The reverse process is known as expressing a rational function in terms of 
partial fractions, and the procedure is useful in differentiation and integra- 
tion, and when expanding a rational function as a power series. 


Suppose that we are told that the identity 
Sig ~ 1). «-" A B 


illienennicpetiiiicactacuianes 4s (1.5) The symbol = is often used 
(2 = See) ea Se in such a case to emphasise 
the fact that the equality 
applies for all values of x, 
There are various methods that can be used to find the values of A and B; but we shall normally 


perhaps the simplest is to substitute various well-chosen values of 2. use = unless we wish to 
clarify a point. 


holds for some constants A and B, and we want to determine these constants. 


|.2 Some initial definitions and results 


For example, putting x = 0 in identity (1.5) gives the equation 
5) A B 


a ae, 
while putting x = 1 gives 

A 

O=-—+-—. 

- = 


These two equations for A and B are sufficient for us to establish that A = 2 
amd = 3. 


Alternatively, we could multiply both sides of identity (1.5) by (a — 3)(a + 2) 
to obtain 

o(¢ =1), = Ate +2) + B(x — 3). 
Then equating the constant terms and the coefficients of x on both sides of 
this expression gives 

—§5=2A-3B and 5=A-+B. 


Again we have sufficient information to establish that A = 2 and B = 3. 
The latter technique is often referred to as equating coefficients of powers 


a 
| 2 oe ie) 
This method works for a rational function —— where the polynomial P(x) 


Q(x) 


is of lower degree than the polynomial Q(z). 


Example |.13 


| +64 — x 


See I OE Bains. 
Qr—D(xt+2)(¢+1h * 


Express 


Solution 


Since the numerator is of lower degree than the denominator, we attempt 
to express the fraction in the form 


i-je=2 > oe + B C 
(Qn—1)(2@+2)\(a@4+1) 22-1 2x+2 x+1 
for some coefficients A, B and C’. Then we have 


1462 — 2? = A(x +2)(2 +1) 4+ BQr=—1)(@+ 1) +CQe—1)(2 + 2). 


Now we either set x = 5, x = —2, x = —1 in turn, or equate coefficients of 
powers of x, to obtain A= 1, B = —3, C = 2. Thus 
1+62 — 2x? 1 3 2 


(2a —1)(a + 2)(a + 1) — haere eee. 


Exercise |.7 


Express the following in partial fractions. 


iee+3 
) Beer) 
2(2 + 1)(2x + 7) 
(x — 1)(a@ + 2)(a + 3) 
x? + 2x x + 2x Pix} 


Hint: Try writing ea = x + Otay: 


(b) 


(c) 


72 — 
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Repeated and quadratic factors 


A minor adaptation of the method is required when a factor in the denomi- 


2x . 
nator is repeated. For example, to express (og in partial fractions, we 
—2 


write 
i i | . B 
(l—z)* 1-2 (1-2)? 


Multiplying by (1 — x)*, we obtain 2x — 1 = A(1 — x) + B, and then equating 
coefficients of powers of x gives A= —2 and B = 1. 


A similar approach will deal with rational functions in which the denom- 

inator has a quadratic factor that will not factorise into real factors. For All quadratics will factorise 
getetl | if we allow complex 

. ie in partial fractions, we write timeless in the factors. In 

(x* + 1)(x — 2) fact, the technique 

get+etil A Bre+C discussed above can readily 

be extended to this case. 


example, to express 


(at hie = eo ee +1 


The value of A can be determined by multiplying both sides by x — 2 and 
then putting x = 2, giving A = f Then we multiply by (x? + 1)(a — 2) and 
equate coefficients, giving B = = amd C= . Thus 


getaeti a 7 =) 


(a2 +1)(x—2) 5 7-2 ee 


Exercise 1.8 


Express the following in partial fractions. 


(a) 


x 
vg? —5a2+6 


P 
Vink: Seiten teen diet 
4 i es 


(c — 1)? 
eee 


£2 
©) @yza0 +e? 


Example 1.14 illustrates some of the applications of partial fractions, but 
first we remind you of two useful facts. 


+a 
1 , 

® (i eee for ie) = %, 
oat’ 


d 
* / : = In(e +a) +¢- dor e+e >, 


Example |.14 


Use the expression obtained in Example 1.13 to do the following. 


1+ 6x — 2 
(a) Evaluate the integral ces Peres dx for x = ‘. 
1 + Ge =<" 


(b) Express f(z) = (Qn — Ile +2)(e+1) as a Taylor series near x = 0. 
d? 
(c) Calculate a3" 


|.2 Some initial definitions and results 


Solution 


1 
(a) Using the fact that / és dx = |nzx for x > 0, we see that 
/ 1462 — x? = 
(2x — Ve +2)(z +1) 


a Se Qe —1° oye +/oie 


= §In(2x — 1) —-3ln(z + 2) + 2In(z# + 1) +, 


since © > 4. 
(b) First we have 
1+ 62 — 2° eet 2 
(Qe—-1)(v+2)(@+1)° 1-22 1+¢2 liz 


then we recall that the geometric series can be written as 


1 
pe ee for ja) < ¥. 
Replacing x in this series first by — 5a, then by 2x, and finally by —z, 
we obtain 


be 3 : ‘i. ree. 
= —3(1-—ga+4+ (52) — (5 --) 
(Qr —1)(a + 2)(e+1) 2( gu + (52) — (gx) + 
4+%1—a2+27-—2°+---) 


—(1+ 2x + (Qr)* + (2x)? +--+) 
for |x| < 5 


d? 
(c) First we notice that if g(x) = (1+ ax)~', then 2 = —3!a3(1+axr)~4 


So, choosing appropriate values for a, we obtain 


af _ 3! Sin te. 
dx? — \ (1—2a)4 e (irda) (1+ 2)4 


48 ie 18 12 
(1—2e)4  (2+2)4 (1+2)* 


1.2.6 Complex numbers 


A brief history 


Historically, complex numbers arose from attempts to splve equations, such 
as x* = —1, which have no real solutions. Early mathematicians would 
have said that such equations have no solutions, but in the mid-sixteenth 
century, Cardano and his contemporaries were performing calculations in- 
volving \/—1. They discovered that assuming the existence of this so-called 
‘imaginary number’ enabled them to simplify the algebra involved in estab- 
lishing results involving ‘normal’ numbers. Euler introduced the notation 
i = /—1 in 1777, and established the relation e’” = —1 (which involves four 
of the most interesting numbers in mathematics). In 1799 Gauss proved his 
fundamental theorem of algebra, which states that every polynomial equa- 
tion has real or complex roots, and early in the nineteenth century Cauchy 
and others extended the notions of calculus to complex functions. However, 
these are generally matters that are better discussed in a course on pure 
mathematics. Here we are mainly concerned with the applications of com- 
plex numbers; we shall meet them again when we discuss Fourier series and 
transforms. 
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Complex algebra 


For our purposes it is adequate to introduce the imaginary number 7 with Ina more formal 
the property that i? = —1, and to regard the set of complex numbers C as_ definition, of the kind that 


the set of all ‘numbers’ x + iy, where x and y are real, with the normal rules 0" ™ay find in a course in 
pure mathematics, the set 


of algebra. Thus, for example, fC teeeliedednodo be the 
(ts 2i)° = (1+ 24)(1 + 24)(1 + 2%) set of ordered pairs of real 
as i numbers (x,y), with 
= 1+ 3(27) + 3(27)° + (2%) certain specified rules of 
—-14+6i-—12-—8 addition and 
multiplication. The 
= =r complex number 7 is then 
identified with the ordered 


Real and imaginary parts a ee 
It is common practice to write z = x+7y and then refer to the complex 

number z, with x and y as its real and imaginary parts, respectively. We 

use the form z = «+ 72y so frequently for an arbitrary complex number that 

it is not usually necessary to add that we are assuming that x and y are 

real. 


An expression such as x +7y is often known as the Cartesian form of 
a complex number z (by analogy with the Cartesian coordinates of the 
plane). It is also useful to introduce the notation Re(z) = x and Im(z) = y 
for the real and imaginary parts of the complex number z. You should note 
particularly that the imaginary part of z= x+y is y, not zy. You should 
also be aware that two complex numbers can be equal only if both their 
real and imaginary parts are equal. Thus, for example, if « + iy = 2+ iv3, 
then we know that « = 2 and y = /3. This process is sometimes known as 
equating real and imaginary parts. 


We also identify the set of real numbers with the set of complex numbers 
having zero imaginary parts. In other words, we regard the real number 3 
as being identical to the complex number 3+ 07. So a real number may be 
regarded as a special case of a complex number. 


The Argand diagram 


Real numbers are often represented as points on a line, whereas complex 
numbers, having two independent components, require two dimensions. It is 
for this reason that they are often represented as points in a two-dimensional 
Cartesian system, and this representation of the complex plane is usually 
known as the Argand diagram (see Figure 1.6). In this representation the 
x-axis is often known as the real line and represents the set of real numbers. 
We may also refer to the xz- and y-axes as the real and imaginary axes, 
respectively. 


If z=x-+1y, then the conjugate of z is defined to be z* = x — iy (see Some authors use Z rather 


Figure 1.6). than z*. 


The modulus of z = x + iy is defined to be |z| = \/x? + y?, and the identity 
|z|? = z 2* is often useful (as you will see in the exercises). 


It is easy to establish the results 
izi=} ead ee 
for any complex number z. Also, 
2122] = |z1| |z2|,  |21/z2| = |zi|/lza|_ and (z122)" = 2729 


for arbitrary complex numbers z; and 2. 


|.2 Some initial definitions and results 


real line 


Figure 1.6 ‘The Argand diagram 


Example 1.15 
2+ 32 
Given z = a find z* in Cartesian form. 
—1 


Solution 
We can convert z into Cartesian form by writing 
ee oe ea ae 
= ae. 5 
1—2 ia. 
then we can see that 


z* = —$(14+ 5i). 


Alternatively, we can replace every 7 by —i in the expression for z, to obtain 


,.. 2-3 
2 = . 
1+2 
and then 
2—31 1-12 1 
-_ = —5(1+ 52). 
PT. Te 
It would be equally acceptable to write 
(oe 
Seal 
— (314+ 5) Or “Sa 


Exercise 1.9 


a+ bi 
Express 27 = . 


c+di 


in Cartesian form (where a, b, c and dare real, with c? + d? # 0). 


Exercise 1.10 


Rewrite the complex number z = (1 + 37)? — (3 — 7)? in Cartesian form a + ib, where 
a and 0 are real. 


Find the modulus and the imaginary part of z. Also find z*, and calculate z 2*. 
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Polar form and de Moivre’s theorem 
Using polar coordinates, z = x + zy becomes 


z=r(cos@+ isin 6), Some authors abbreviate 


cos @+ isin 6@ to cis 6. 
known as the polar form of z, where 


= 
r=vV/a7+y?%=|z2|>0, cosé = ————, sin#= a 
/ 2 ie y* / 2 wa y? 
On an Argand diagram, the modulus r is the distance of the complex 
number from the origin, and the angle @ is known as its argument (see 
Figure 1.7). 


It can be shown that 
Zz” = r"(cosné + isin né) 


for any positive or negative integer n, a result which is often known as 


de Moivre’s theorem. Abraham de Moivre (or 


a , . Demoivre) was an 
Note that it is possible to define powers of complex numbers z" for any real anglicised Frenchman and 


number n. In fact, de Moivre’s theorem also holds in this case, as will be a contemporary of Sir 


demonstrated in Exercise 1.13. Isaac Newton. He adopted 
various spellings of his 
name while he lived in 


Exercise |. | England. 
Use de Moivre’s theorem and the fact that cos? 6 + sin? @ = 1 to show that 
cos 59 = cos 6(16 cos* 6 — 20 cos? 6 + 5). 


The argument of a complex number 


The polar form of a complex number is not unique, because its argument 
is not unique. For example, the complex number 1 +7 can be written in 
polar form as V2 (cos 7 +7sin F), but it would be equally valid to write 
V2 (cos 27 + isin 27), or indeed V2 (cos(4 + 2km) + isin(F + 2k7)) for any 
integer k. 


Figure 1.7 ‘Two possible values of the argument of a complex number 


An argument of an arbitrary complex number z = x + 2y is defined to be an 
angle @ for which sin@ = y/\/x? + y? and cos@ = x/\/x? + y?. 


If z=r(cos#+isin@) and —7 < 6 <7, then @ is known as the principal Some authors define the 
value of the argument of z. An argument of z is written as arg(z), while principal value to lie in the 
the principal value of the argument is written as Arg(z). interval 0 < @ < 2r. 


Example |.16 


Express the complex number z = V3 +i in polar form. 


|.2 Some initial definitions and results 


Solution 


If # represents an argument of z, then cos @ = 3/2 and sind = 1/2, and we 
see that 7 is a suitable value for 6. Also, |z| = 2, so 
“fe = 2(cos% +isin#). 


Any value of the argument of the form ¢ + 2k7, where k is an integer, is 
equally acceptable, but @ is the principal value of the argument. (When 
determining the argument of a complex number, it is often helpful to plot 
the point on an Argand diagram.) Mf 


We can use de Moivre’s theorem to find roots of complex numbers. For 
example, to find the fourth roots of z = V3 +i, first express it in polar 
form, 


z = 2(cos(Z + 2k) + isin(% + 2mk)) , 


using the general value of the argument, i.e. with k as any integer. Then 
use de Moivre’s theorem to evaluate the fourth roots: 


2/4 — 91/4 (cos ((2 + 2k)) + isin (4(% + 2zk))). 


For & = 0,1, 2,3, this gives four distinct complex numbers: 


_ 91/4 cm saat. ie 1/4 pee 137 
p— 2’ (cos 4 + sin 5), = 2 ee eae), 
1/4 257 257 1/4 i 377 
= 2'/ (cos = + isin “se = 21/ (fos 2, + in 


which are the four fourth roots of z. Other values of k simply duplicate these 
roots — which can, of course, be expressed in Cartesian form if required. 


Exercise |.12 


a) Express the complex number z = 2 — 22 in polar form. 


b) Find the cube roots of z = 2 — 2% in Cartesian form. 


The complex exponential function 


It is possible to extend the notion of Taylor series to complex numbers, and 
to define 


exp(i0) = e”” shee + 


Ea e 6 ie. Ox, ee 


= cos? + isin @. 


This result is known as Euler’s formula. 


Any complex number in the form z = r(cos@+isin@) can then be written 
in the equivalent form z = re’’. 


We can also use Euler’s formula to show that 
eG rs e8 ef ae e 8 
cos 8 = —————-_ and_s sin = ——— 
2 at 


(which reinforces the close relationship between the trigonometric and hy- 
perbolic functions). 


Exercise 1.13 


Use Euler’s formula to prove de Moivre’s theorem. 


25 


You should know that a 
complex number has at 
most n distinct nth roots. 
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Example I|.17 


Express the following numbers in polar form. 


ck wae 


1+2 
Solution 
(fs) Lc eo (cos 4+ 7sin t) In each case a multiple of 
. ae 27 could be added to the 
(b) 1+2/3=2 (cos 3 +tsin 5) argument. 
(c) Since (using the results of parts (a) and (b)) 1+i/3 = 2e’*/ and 
1+i= V2e'"/4) it follows that 
1 +iv3 le te/i2 ™ ne 
7 ee = V2(cos4 +isin4). @ 
Complex numbers are often useful in calculations that would otherwise be 
rather tedious. Suppose, for example, that we wish to evaluate 
/ e™ sin(bt) dt, 
where a and bare real. This integral can certainly be found using the method 
of integration by parts (twice), but the following approach is probably easier. 
Consider the integral f. Ls Re 
a, es re dt wi ~— 
\ et = LSO+ eswn& 
= few dt , okt Olt 2K = endbd + Usrrlot) 
comer. SMe 
I si a. 
= _ + jee ; ak GERACE 
re +C a7. (ox ( S ) 
¢ gat Cos(bt) + isin(bt) tag 
rom a+ 2b 
ae cos(bt) + = a a 2 LC 
a+b a — 1b 
ert 
=a ((acos(bt) + bsin(bt)) + i (asin(bt) — bcos(bt))) +C. 


We can also write 
: L . ak. - fpr 
j= [ertew at = fe (cos(bt) + isin(bt)) dt, = {< we (be)* i | oer Cot) 


so, equating real and imaginary parts of these alternative representations 
of J, we obtain / 


t ; 

t = = | 
Je cos(bt) dt = aR (a cos(bt) + bsin(bt)) + C; | Vint WA thu 
t - | 4 [« 
Je sin(bt) dt = ee) (a sin(bt) ae b cos(bt) ) + Cp \ ak ¢ CG UNS 


(where C = C1 + 1C4). 


So we have found not only the required integral, but also / e™ cos(bt) dt 
with no extra effort. 


In a later chapter we shall meet complex numbers again in the context of 
Fourier series, and the following observations will prove useful. 


|.2 Some initial definitions and results 


For any non-zero integer k, 
r otk? dé = 3 [et**] = + es x) =) ae 2 sin(k7) > 6 (1 6) 
ee tk —~ k; 
eg ae 


while for k = 0 we have 


/ Mao = | dO = 2n. (47) 


To give you some impression of how these results will be used later, consider 
the function 


f (0) co Aye™ + Ase he Aner, 


With the aid of the stated results, it is very easy to express the coefficients 
in this expression in terms of the function f. For example, to determine As, 
multiply by e~°”" and integrate from —7 to 7, to obtain 


f(0)e—>” do = Ay / ed) + A / e"” do 


== Al =i 


+ As / d0+ A, / ef dé, 


= h =i 
ee ee 
which gives A3 = = f(@)e >” dO. In general, 
as a 


: f(d)e—™ do 


a 
AS: Joop 


ta = 1, 2,5,4. 


A 


End-of-section Exercises 


Exercise 1.14 

Prove the triangle inequality, that |a + b| < |a| + |b| for any two real numbers a 
and 0b. 

Exercise 1.15 


Show that 


d d 
7, (sinh o)=coshs and Fa (cosh o = siek x. 


Hence show that @ = Asinhwx + Bcoshwz, where A and B are arbitrary constants, 


is a solution of the differential equation = wd. 
L 


Exercise 1.16 


Express 


1 
n(n + 1) (n+ 1) 


CO 
in partial fractions, and hence show that a ——__— = ]. 
n 
n=t 


Exercise |.17 
Find the Taylor series for sinh(2”) near x = 0. 


Exercise 1.18 


Show that = : 


tas + 22° + O(2°). 


Exercise 1.19 


I 


is odd. Hence evaluate / paw 523 dx. 
24142’ 


Show that f(x) = aa : 
x 
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Exercise 1.20 


d 
Show that = (sinh™ 2) = and sketch the graph of y = sinh™' «. 
x 


1+ 2? 


Show that r+ Vl+2?>0ifzeER. 
d 
Calculate . (In(a + V1+ 2?)), and hence show that sinh” * ¢ = In(x + V14 22). 


Exercise |.2] 


Ag? — 9a +4 
x(x — 1)(a — 2) 


4 2 
Ag* —9x+4 
in partial fractions, and hence evaluate / ee da. 


E 
xpress gl ee 


Exercise |.22 


CO 
For what values of the constant a is the integral / es dt convergent? 
1 


Exercise 1.23 This exercise is optional. 
CO CO 

Show that if / | f(t)| dt is convergent, then / f(t) dt is also convergent. Applied mathematics is 
a a dependent on many results 

from pure mathematics. 
re | This exercise is intended to 
~~ nt 8, give you some indication of 

Show that = dt is convergent. the steps that would be 

a 


required to justify some of 


Exercise 1.24 our assumed results. 


Express the following in Cartesian form. 
1+ 32 2 — 32 1 1 

b 
2 a |) ! ape 


Exercise 1.25 


Express the following in polar form. 
1 


(a) <b) (1+ v8)" (@) (1 +4v8) 


Exercise 1.26 
Show that cosh(i0) = cos@ and sinh(i@) = isin @. 


Exercise |.27 
Show that for any two complex numbers z and w, 


lz +.w|? + |z — w|? = 2lz|? + 2lew/*. 


1.3 Differential equations 


A brief history 


The study of differential equations has a venerable history, and the subject is 
still central to modern mathematics because many physical systems can be 
described by differential equations. As early as 1676, Newton solved a differ- 
ential equation by means of an infinite series, although his results were not 
published until 1693, the year in which a differential equation appeared for 
the first time in the work of Leibniz. Shortly afterwards, Johann Bernoulli 
developed the method of ‘separating the variables’, and he also showed how 


|.3 Differential equations 


a differential equation of the form dy/dx = f(y/x) could be reduced to one 
in which the variables were separable. He and his brother Jakob succeeded 
in reducing a large number of differential equations to soluble forms. 


Developments in the eighteenth and nineteenth centuries were rapid and 
involved some of the greatest names in mathematics: Euler, d’Alembert, 
Lagrange, Laplace, Cauchy, Jacobi, Ampere, Frobenius and Lie, to mention 
but a few. 


From these passing references to its history, it should be evident that the 
study of differential equations is a vast subject, and here we can do little 
more than scratch the surface. 


Introduction 


While it is important to be able to solve various differential equations, it is of 
equal importance to be able to recognise that a given equation is of a certain 
type, for then you can refer to the method of solution in one of the many 
standard texts. You may have seen some of these differential equations 
before in other courses, in which case you can regard the corresponding 
exercises as revision, but they are included here for completeness. 


While the variables, constants and parameters in this section are generally 
real, much of this material applies also when they are complex. However, 
certain complications arise in integration of complex functions, so for the 
moment you may assume that everything is real. 


1.3.1 Initial classification of differential equations 


Essentially, the study of differential equations is a matter of assigning them 
to various categories and discovering as much as possible about each such 
category. Here we shall make a start by dividing them into two major groups: 
ordinary differential equations and partial differential equations. 


Ordinary differential equations arise when we are considering the variation 
of a physical quantity that depends on just one independent variable. If we 
denote the independent and dependent variables by x and y, respectively, 
then an ordinary differential equation may involve functions of x and y, and 
the derivatives of y with respect to xz. Thus 


a d d 
—: + oe +y=0, = + cy = 0 (where c is a constant) 
d*y ae 
and , Pe —k* siny (where k is a constant) 


are all examples of ordinary differential equations. (Quantities such as c 
and k in these illustrations are sometimes known as parameters; while the 
solution may depend on their values, they are generally assumed to be fixed.) 
Sometimes we omit the word ‘ordinary’ when it is clear that we are referring 
to an ordinary differential equation. 


On the other hand, partial differential equations arise when the dependent 
variable is a function of two or more independent variables. If we denote 
the independent variables by x and t, and the dependent variable by y, 
then a partial differential equation may involve partial derivatives of y with 
0? ome 
respect to x and t, as well as functions of x and t; for example, ae = oo 
= 
(where c is aconstant). We shall be concerned with such equations and their 
solutions later in the course. 
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Here we are primarily concerned with the process of finding solutions, not 
with the difficulties that may arise in the definitions of the functions and 
their domains. We assume throughout that unless otherwise stated, the 
functions involved are well behaved: for example, we assume that they are 
differentiable as many times as we please on an interval of the real line. (You 
may take this to mean that the functions are ‘sufficiently smooth’.) 


Solution of a differential equation 


A function that satisfies a given differential equation is known as a solution 
of the differential equation, and, while it may not be obvious how one might 
find such a function, it is easy to verify that a given suggestion is correct. 
For example, we can easily verify that the function y = 3x7 + Ae” satisfies 


d 
the differential equation -~ — y = 6x — 32? for any value of the constant A, 
by writing is 


— — y= (6x + Ae”) — (32? + Ae*) = 6x — 32”. 


A solution may be required to satisfy certain additional constraints, for 
example that (0) = 0, and such conditions will further restrict the choice 
of a suitable function. 


It is not always possible to find the required function explicitly; we may have 
to be content with an implicit definition. For example, it is easy to verify 
that if y = y(x) is a function for which x? + ry + y°? = c, then differentiating 
with respect to x gives 


d d 
oe ty te + hg ew. 
dx dx 
dy y + 22 = ; 
It follows that — = ————,, so y = y(2) is a solution of this differential 
dx x + 5y* 


equation; but we cannot find an explicit formula for y in terms of 7. We may 
also refer to an implicit definition, such as x? + ry + y° = c, as a solution of 
a differential equation. 


Example 1.18 


Show that the function y(x) = 3sin(Az) is a solution of the differential equa- 
d? 

tion — — —)*y, where ) is a constant. 
dx 


Solution 


In order to establish that this is the case, we need only differentiate y twice, 


d 
to obtain first = = 3A cos(Ax) and then 
2 


d*y ' 
- —3)" sin(Ax) = =—A*y; 


as required. Hl 


In the previous example the solution was given as a function, but more often 
such solutions arise implicitly, as in the following case. 


Example |.19 


x 


y—e 
eu wig 


d 
Show that ry = e* + e¥ is asolution of the differential equation = = 
ss 
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Solution 


Differentiating the equation ry = e” + e” with respect to x, we have y+ xy’ 
=e”+e%y’. Rearranging this equation gives y/(eY — x) = y — e”, which 
establishes the required result. 


Notice that in this case we are not able to rearrange the original equation 
into a simple algebraic expression for y in terms of z. Hf 


Order of a differential equation 


The order of a differential equation is the order of the highest derivative that 


d? dy \° 
appears in it. For example, the differential equation — + (H) +y=0 
< x 


is of order two. 


Exercise 1.28 


(a) In each case, form a differential equation of first order by eliminating the 


constant a. 
i ge" (ii) y = sin(a@ +a) 
(b) Form a differential equation of second order by eliminating the constants a and 
b from y = — a 
ae es 
d?y 
(c) Solve the differential equation i i‘ 
- 
(d) Show that a function y = y(x) that satisfies the equation ry + cosy = 3 is a 


d 
solution of the differential equation a — 
dy sny-— x 


Linear differential equations 


A linear ordinary differential equation of order n is an equation that 
can be written in the form 


Wy 5 ciple atch 110), 


dx 
with a,(x) ~ 0, for some function f(x). 


A solution of a linear differential equation of order n containing n indepen- 
dent constants is known as the general solution of the equation. For linear 
equations it can be shown that every solution can be obtained by assigning 
values to the independent constants in the general solution. This is not the 
case for non-linear differential equations. 


The differential equation (1.8) is said to be homogeneous if f(x) = 0, 
and such an equation has n independent solutions y;(x), y2(x),...,Yn(2). 
Moreover, every solution y(a) of the differential equation may be written in 
the form 


ule) = Agyalz) 4+ Ages -+ - ++ Ange (2), 
where Aj, Ag,.. 


(1.9) 
., An are constants. 


In order to find the general solution of equation (1.8) when f(x) #0, we 
need only find a single solution without arbitrary constants (known as a 
particular integral), then add a function of the form (1.9). 


3| 


From the chain rule, the 


derivative of e¥) with 
respect to x is ye. 


With the aid of a computer 
program, it is possible to 
convince oneself that the 
equation zy = e* + e4 
defines a function y = y(x) 
with the negative real axis 
as its domain. 


Part (c) illustrates the fact 
that we might expect the 
solution of a second-order 
equation to contain two 
arbitrary constants. 
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In order to solve a practical problem, it is not usually sufficient to determine 
the general solution of a differential equation; often we need to assign specific 
values to the constants that arise in the general solution. Such values are 
generally determined by certain additional restrictions; loosely speaking, 
one would expect that n pieces of information would be needed to uniquely 
define the solution of a differential equation of order n. 


As a simple illustration, y = Ae’ + Be” is the general solution of the 
second-order differential equation y” = y. The conditions y(0) = 1 and 
y’(0) = 0 are sufficient to determine the two constants A and B (in fact 
Pea 5). This is an instance of initial conditions, which we shall discuss 
shortly. 


1.3.2 Elementary solutions of ordinary differential 
equations 


Separable equations 
A first-order differential equation is said to be separable if it can be written 


d 
in the form — = g(x)h(y), and the general solution of the equation may be 


obtained by evaluating the integrals in the equation 


bree = | g(x) ae. 


Example 1.20 
To obtain the solution of the equation y’ = xy, we rewrite it in the form 
y'/y = x, then integrate to obtain [ow 2 [uly = [oar Evalu- 


ating the integrals gives Iny = 5x" +c, where c is an arbitrary constant. In 
fact, the solution is made simpler by choosing the arbitrary constant to be 
c= 1nA, where A > 0 is also arbitrary, for then Iny = 5x" +n A, and the 
solution is y = Aexp (527). (This is a solution of the original differential 
equation for any choice of the constant A, i.e. it is the general solution.) 


Example 1.21 


x 


d 
The general solution of the differential equation eed ee may be found 
dx 1+y? 


by evaluating / (1 + y”) dy = Je dx to give y+ Sy? =e*+c. This is an 
implicit equation fory. Hi ° . 


The solution of a differential equation that satisfies certain additional con- 
straints is known as a particular solution. 


Example 1.22 


d 
To find a particular solution of the differential equation = + x= 0 for 
2 


which ¥(0) = 3, we first recognise that the equation is separable, then find 
the general solution from the equation | ydy = — | xdzx. Thus the gen- 


eral solution is x7 + y* =, and it only remains to find the value of the 
constant c. Putting y= 3 when « = 0, we have c= 0? + F429. so the 
particular solution is x? + y* = 9. 
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Alternatively, we may notice that the solution may be written in terms of 
definite integrals so that the initial conditions are automatically satisfied. 


Yy - 
We have / yay = — / xdx, and evaluating the integrals gives the re- 
3 0 


quired result directly. This is a useful trick that is worth remembering. & 


Example |.23 


A boating pond is 10m wide and 20m long. A boy stands in one corner; his 
toy boat, attached to a string 10m long, is at an adjacent corner, as shown 
in Figure 1.8. He walks along the long side of the pond and tows the boat, 
keeping the string taut. Find the equation of the path followed by the boat. 


Figure 1.8 ‘The path of a toy boat 


Solution 


We use axes as shown in the figure. When the boat is at the point (x,y), the 
boy has moved a distance x + \/100 — y?, and the slope of the hypotenuse 


of the shaded triangle in Figure 1.8 is —y/,/10? — y?. Thus the gradient of 
the path followed by the boat is 


dy _ y 


de ~~ fio y 


SO 


1 — 
-[ae- | a 


The integral on the right-hand side can be evaluated using the substitution 
u* = 100 — y? for u > 0, which gives 


/ Se i (-ge4) = es te 


Now notice that y = 10, and therefore u = 0, when x = 0. So, using partial 
fractions and the trick mentioned in the previous example, we have 


ac U wu U Fie 1 
_}) ef — a = d 
/ F / ae tee | ( (a+ ex)) “ 


Thus 
—¢ = fut+5In(10-u) —5In(10+ u)], =ut5in — 
" 10+u/’ 


which gives 


1 10 + ./100 — y2 
r=5in( se) -u=5in te ee ag a 
wu 10 — ,/100 — y? 


In this example we cannot rearrange the equation to give y = y(x), but we 


could use a computer program to plot the graph of x = x(y) as shown in 
Figure 1.9. 
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Note that in an integral 
y 


such as f(y) dy, we use 


the same rs for the 
integrand and the upper 
limit. This expression 

should be interpreted as 


(freee)|. 


The ‘shorthand’ is 
frequently useful because it 
saves on substitution of 
variables. 


The integration after the 
substitution uses the 
method of partial fractions 
as discussed in 

Subsection 1.2.5. 


| 
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0 5 10 15 20 ‘3 


Figure 1.9 The graph of « = x(y) (with the independent variable y on the vertical 
axis) 


Homogeneous equations of first order 


The word ‘homogeneous’ is commonly used in mathematics in several dif- 
ferent senses. Here we define a first-order differential equation to be homo- 


d 
geneous if it is an equation of the form a = (2). 
x x 


In order to solve such an equation, we write y = xv and endeavour to find a 
differential equation in v and x that is separable. Differentiating with respect 


d d 
to x, we obtain 5 = = +v, and eliminating y from the original equation 
a ze — 
d Vv vy He 
gives Ze +v = f(v), which can be rearranged into the form = eee 
2 x e 


which is separable. 


Example 1.24 
In order to solve the equation 


dy _ ai+y? 


dx a 


we first notice that it is a homogeneous first-order differential equation, since 
it can be written as 


d ¥ 
ac 8. x 
So we make the substitution y = xv, giving y’ = xv’ + v as above. Elimi- 


nating y gives 
dv i a 
Se = ee 
dx Ss wv | 


and the differential equation becomes 


dv 1 / [2 
x—=-, SO vdv= | —. 
a: »% oF 


Hence Su" =Inx+A for z > 0, and finally replacing v by y/x gives 


y\2 
(=) =2Inxr+C, 
> 


where C' = 2A is an arbitrary constant. 


Exercise 1.29 


In each case, show that the equation is homogeneous, and solve it. 
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2 dn 
Exercise 1.30 


Show that the differential equation 
me eters 


dx y 


is homogeneous, and use the substitution y = xv to solve it. [Hint: Try the further 
substitution u = V1+4 v?.] 


Linear equations of first order 
The general linear equation of first order takes the form 


ay(2) + ao(2)y = F(a). 


For example, the equation 


d 
3-2 + 3x7y = 2 
dx 
is of this form, and in this case it is easy to find the solution if we notice 


that 


d d 
gine 4. 3a7y = —(a°y). 


dx dx 
Then the differential equation becomes 
d 3 5 
, y) —s 4 


and we have only to integrate this expression to obtain the solution 

vy = ax° a 
However, it is not immediately clear how we would proceed if the same equa- 
tion were to arise in a slightly different form, for example, = + 2ry = 62. 


Fortunately, there is a simple method that will transform such an equation 
into a form that can be integrated. 


The integrating factor method 

Any linear differential equation of first order of the form 

d 

— + g(x)y = h(a) (1.10) 


may be solved by multiplying both sides of the equation by the integrating 
factor 


P(x) = exp ( / g(a) ir) | 


The differential equation then becomes 


* (yP(a)) = h(2)P(o), 


which may be integrated to give the general solution 
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Example 1.25 


d 
Find the general solution of the differential equation — + 2xry = 62. 
z 


Solution 


This equation is of the same form as equation (1.10) (i.e. first-order linear), 
with g(x) = 2x and h(x) = 6x. The integrating factor is 


P(x) = exp ( / 2x ir) = exp(z’). 


So the general solution is 


2 = 
[o exp(x ) dx = exp(22) 
= 3+43Aexp(—<’). 


= ape (exp(x*) + A) 


Example |.26 


Find the particular solution of the differential equation 
dq q Differential equations of 
= + a E(t), this kind arise in the 
t theory of electrical circuits. 
where R = 10, C = 10-9 and E(t) = 20sin(100zt), that satisfies q(0) = 0. 
Solution 


The differential equation becomes 


d 
10 + 102g = 20sin(100zt), 


which simplifies to 


d 

= + 100g = 2sin(1007t). 

This equation is first-order linear (see equation (1.10)), with g = 100 and 
h = 2sin(100zt). So the integrating factor is P(t) = exp (/ 100 it) ag, 


and multiplying the differential equation by the integrating factor gives 


d 
= (e9) = 2e1 sin(100zt). 
Integrating, we obtain We discussed integrals of 
— the form / e™ sin(bt) dt on 
ROO eee (sin(100zt) — mcos(100zt)) + C-. oe 
50(1 + 1?) page 26. 


Using the condition gq = 0 when t = 0 gives C = 7/50(1+4+ 77), and finally 


= Sy (sin(100mt) — weos(100mt) + me"). 
Tv 


Exercise 1.31 

(a) Find the particular solution of the differential equation a + y coba-= Be" * 
that satisfies y = —1 when x = 5. 
(b) Find the general solution of the differential equation (x — 1)(a — 2) et +y = 2". 
(c) Show that the transformation v = y'~", where n # 1 is a constant, reduces 


d 
the differential equation _ + yF (x) = y”G(z) to linear form. 
. 
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Linear equations with constant coefficients 


Linear differential equations with constant coefficients are particularly easy 
to solve. We shall concentrate on equations of order two, but higher-order 
equations can be solved in a similar fashion. (For order greater than four, 
the auxiliary equation cannot generally be solved algebraically. ) 


We consider differential equations of the form 


d7y dy 

a5 +b + cy = fle x), 
and we assume that you have seen such equations in a previous course; these 
brief notes are intended primarily as revision. 


where a, 6 and c are constants, 


Procedure 


The above differential equation is said to be homogeneous when f(a) = 0. 
In order to solve the homogeneous equation 


2 

at + pl +cy =0 

where a, b and care real constants and a # 0, we use the following procedure. 

e Write down the auxiliary equation am* + bm +c = 0, and find its roots 
my, and mo. 

e Ifthe auxiliary equation has two distinct roots, then the general solution 
ng= Ac + he", 

e If the auxiliary equation has two equal roots, i.e. m1 = m2, then the 
general solution is y= (A+ Ba)e™™”. 

e If the auxiliary equation has a pair of complex conjugate roots, given 
by my =a+t+i@ and mz =a-—if, then the general solution is 
y = e**(Acos Gz + Bsin Ba). 


(In each case, A and B are arbitrary constants.) 


In order to find the general solution of the second-order inhomogeneous 
linear equation with constant coefficients 


ey 


dy 
afd + 5H 4 oy = Fl), (1.11) 
we first use the procedure above to find the general solution y.(x) of the 
d? 
associated homogeneous equation 5 + = +cy =0. Then we find a 


particular solution yp(x) of equation (1.11). Here y-.(x) is known as the 
complementary function, and yp(x) as a particular integral. 


The general solution of equation (1.11) is then 


y(x) = ye(x) + yp(2). 
For an nth-order differential equation, y will be a solution containing n 
arbitrary constants, so it will indeed be the general solution. 


There are various techniques for finding particular integrals yp, and we refer 
you to any of the standard texts on differential equations for the details. 
Here we illustrate one of these techniques, known as the method of un- 
determined coefficients. In order to apply this method, we need to be 
able to guess the form of a solution and then determine the constants in our 
proposed form by equating coefficients. The following example illustrates 
the method. 
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Note that the term 
‘homogeneous’ as used here 
differs from its application 
to first-order equations on 
page 34. 


If the coefficients a, b and c 
are real, and the auxiliary 
equation has complex 
roots, then these roots will 
always occur in complex 
conjugate pairs. 
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Example |.27 
To find the general solution of the differential equation 


d*y dy 

7 — + 6y = 10, 

we consider first the auxiliary equation m? — 5m +6 = 0, which has roots 
m, = 2 and m2 = 3. Hence the general solution of the homogeneous equa- 
tion (i.e. the complementary function) is y, = Ae?” + Be®”. In order to find 
a particular integral, we guess that a constant function yp(z) = k may be 
suitable, and substituting this into the original differential equation gives 
6k = 10, sok = 2. It follows that the general solution of the differential 


3 
equation is 
yy + yg Ae’ +e = 2. 


This is the general solution of the second-order differential equation because 
it satisfies the equation (this can be checked by substitution), and it contains 
two arbitrary constants. 


If the number 10 in the previous example were to be replaced by 10 + 3x + 2”, 
then we might guess that a+ bx + cx” is a suitable form for a particular 
integral. 


The method may also be extended to differential equations of higher order. 


Example 1.28 
4 


Find the general solution of the equation =, —y= 0. 
x 
Solution 


This equation is homogeneous, so it is not necessary to find a particular 
integral. We consider the auxiliary equation m* — 1 = 0: factorising the 
left-hand side gives 


m* —1 = (m? — 1)(m? +1) = (m—1)(m +1) (mn — 1) (m +: 4), 
so the general solution of the differential equation is 
y = Ae* + Be“ +Ccosx+ Dsinz. 


This is the general solution of the fourth-order differential equation, because 
it satisfies the equation and contains four arbitrary constants. Il 


Exercise |.32 


(a) Find the general solution of the differential equation 


d°y | dy 
—= + — —6y = 0. 
a” 
(b) Given that y; = —xe” is a solution of the differential equation 
d*y dy 
_) 29a 
dx? dx = -s 


find the general solution. 


In each of Exercises 1.33-1.35, a particular integral may be found by using 
the method of undetermined coefficients. Try to guess the form of a solution, 
then equate coefficients to find a particular integral. 
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Exercise |.33 


Find the general solution of the differential equation 


d’y  _dy 
—— 5 Saga ae 
Ee ee _ - - 


Exercise |.34 


Find the general solution of the differential equation 
d*y 


Exercise 1.35 


Find the general solution of the differential equation 


d*y dy 
ee, ee ee oe on 
zo a: 1 i oe 3 


Exercise 1.36 


Find the general solution of the differential equation 


d*y 
—— == {). 
d+ Ty 


Initial conditions and boundary conditions 


We have previously discussed cases where a particular solution of a differen- 
tial equation is defined by the value of the function at a specific point, but 
sometimes we require more than one additional piece of information. As an 
illustration, consider the equation defining simple harmonic motion 

d7y 

dt? 
where w is a specified real positive constant, y denotes the displacement of 
an object, and the independent variable t represents time. The general so- 
lution of this equation is y = Acoswt + Bsinwt, and we require two further 
conditions if we are to determine the values of the constants A and B. Of 
ten such conditions arise in one of two forms: initial conditions, in which 
the value of the function and its derivatives are specified at some chosen 
value of t, or boundary conditions, in which the function is defined at 
two different values of t. Boundary conditions will be dealt with later in the 
course. 


= wy, (1.12) 


Suppose that an object performs simple harmonic motion and satisfies equa- 
tion (1.12), and that we are given the position yo and velocity vg of the object 
at time t = 0; in other words, we are given initial conditions for the motion 
of the object. We can easily determine that A = yo and B = vo/w, and this 
provides us with the particular solution 


y(t) = yo coswt + (ug/w) sin wt. 


Here initial conditions are sufficient to determine a unique solution for any 
choice of yo and vo. This is generally the case; provided that we are given 
an appropriate number of initial conditions, the solution is always uniquely 
determined. An nth-order equation requires n initial conditions. 
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End-of-section Exercises 


Exercise 1.37 


d 
Solve the equation = e**tY, given that y(0) = 0. 
2 


Exercise 1.38 


d 
Solve the equation = = —2zxtany, given that y(0) = 
x 


Nis 


Exercise 1.39 


20y 
3x2 — Qy?’ 
equation is homogeneous, and use the method of partial fractions.] 


d 
Solve the equation = = given that y(1) = 1. [Hint: Show that the 
2 


Exercise 1.40 


d 
Solve the equation — = y(1+Iny—Inz). [Hint: Show that the equation is 
a 


homogeneous, and solve the resulting integral by a suitable change of variable.) 
Exercise |.41 


iy , a 
Solve the equation x (3) — dy —2 —~¢ =U. Hint First solve a quadratic equa- 
s 


dx d 
tion, then remember that s (sinh™* *) : | 
10N, r . es 
dx JV1+<2? 


Exercise |.42 


d? d 
Find the general solution of the equation a oe 4:19y = 1207 + 10¢ +11. 
2 x 


Exercise 1.43 Harder exercise 


(a) Find the family of curves for which the length of the part of the tangent between 
the point of contact (x,y), in the first quadrant, and the y-axis is equal to the 
y-coordinate of the intercept of this tangent with the y-axis. (In other words, 
find curves where the lengths OA and AP shown in Figure 1.10 are equal.) 


Figure 1.10 


(b) Find the equation of the orthogonal family of curves (i.e. the curves such that 
the product of the gradients at the point of intersection with each member of 
the previous family is equal to —1). 


Exercise 1.44 Harder exercise 


A chemical compound C is being formed by the reaction of two substances, A 
and B. During the reaction, substance A combines with substance B in the ratio 
2:1 by mass to form substance C. The rate of formation of C is proportional to 
the product of the masses of the quantities of A and B that remain uncombined. 
Initially, there are 200 grams of A, 100 grams of B, and 0 grams of C in the mixture. 
Given that 100 grams of substance C are formed after 30 seconds, find the amount 
z(t) of substance C formed after time t. 


|.4 Outcomes 


1.4 Outcomes 


After studying this chapter you should be able to: 


understand the terms continuous, differentiable and order when applied 
to functions, and apply the O notation; 

understand what is meant by the convergence of an infinite series; 
differentiate expressions involving elementary functions; 

evaluate the Taylor series for a function, and understand what is meant 
by its radius of convergence; 

integrate elementary functions; 

understand what is meant by the convergence of an integral over an 
infinite range; 

express rational functions as partial fractions, and apply this technique 
to the integration of rational functions and the expansion of rational 
functions in series; 

manipulate expressions involving complex numbers; 

transform between the Cartesian, polar and exponential forms of a com- 
plex number; 

calculate powers and roots of complex numbers using de Moivre’s theo- 
rem; 

understand the terms ordinary, partial, linear, order, general solution 
and particular solution when applied to differential equations; 

solve first-order ordinary differential equations which are separable, ho- 
mogeneous or linear; 

solve higher-order ordinary differential equations which are linear with 
constant coefficients; 

apply initial conditions to general solutions, to fix the arbitrary con- 
stants. 
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Solutions to Exercises in Chapter | 


Solution I.l 


We have 
CO 
H,{a) = > a =g” pa 4 g's Sai ee es. 
b=n 
For —1 <a < 1, using equation (1.2) gives 
a” 
R,(a)= a" tes os ; 
ot 
Solution 1.2 
From Equation (1.1) we have 
; n+l 
1-—=z 
Replacing « by —2x? gives 
- ia ites 
2 ae 
k=0 


1 
which converges to the limit ieee: N — oo if |x| < 1, but diverges if |x| > 1. 
iy 


CO 

Thus the infinite series S| (-1)*a7* converges if |x| < 1, but diverges if |x| > 1, so 
k=0 

the largest interval |x| < p corresponds to choosing p = I. 


This exercise is one of the very rare occasions when we are actually able to determine 
the function f in a simple form from the coefficients {a,}, for in this case 


? + (to dai 


= |i a 
f(x) umn 14+ 22 1+ 22 


N—-co 


Solution 1.3 


From the Taylor series we have sina = x + O(x*) and In(1+ 2) = 2 +O(2?). It 
follows that 


sin zIn(1 +) = [x + O(x?)|[x + O(x?)] = x? + O(c”). 


Solution 1.4 
We have 


a —x\ 2 zr ,-@\ 2 
cosh? x — sinh? x = (<4) — (SS) 


= e2f 1 94 e— 22 e2u —2+4+¢6-22 x 
‘a 4 4 ee 


ae 1 
(a) be = / : dt and ’ <1 for t >1. From the stated result we therefore have 
1 


Solution 1.5 


ee 
Ine = [ —dt<a—-l1l<z. 
ae’ 


(b) The formula for integration by parts may be written in the form 
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(i) Choosing u(z) = Inz and vu(a#) = Ina, we have 


l l 
[a= (ina)? — f 2 ar. 
x 


which can be rearranged to give 


© 
So 
$ 
[ ~Sac= 47) 
i : 
(ii) Choosing u(x) = Inz and v(x) = —x~1, we have 
[a =-="4 [5a- es 
r r 
So 


T 
In x in 7’ 1 1 In T 
ee ee Se fe f -] pe] eee, 
ch ye ( T 1) = : 


which increases without bound as J’ increases, so the integral is divergent. 


(ii) From part (b)(ii), 


In x 
so the integral is bounded. Since also —~- = 0, the integral is convergent. 
. 


Solution 1.6 


For 1<2a<T’, we have 0 < e-@ < re-® so using the stated result gives 
g 
rT os T 1 1 
/ e-™ dr << / ze dx = [—3e-*"| ee ee ae a a, 
1 1 


CO 
Since this integral is bounded, and e-™ > 0, we see that the integral / e~® dx 
1 


CO 
is convergent. Using the hint, the convergence of / e~* dx follows immediately. 
0 


Solution 1.7 


(a) We write 
llz+5 — - B 
(Qc —1)(83a+2) 27-1 32442 
Then putting c = 0 gives —5/2 = —A+ B/2, and putting + = 1 gives 16/5 = 
A+ B/5. Solving this pair of equations for A and B, we obtain A = 3 and 
Ose. 1, oe 
tif + 5 ee, = 1 
(Qn—1)(3a+2) Qr-1 32427 


Chapter | 


(b) We write 


B 
L+2 


# 
eo 


2(a + 1)(22 + 7) 
(x —1)(a@ + 2)(x + 3) 
By substituting suitable values for x, or multiplying both sides by the de- 


nominator (x — 1)(a + 2)(a +3) and equating coefficients of powers of x, we 
find 


A 


a (1.13) 


A=3, B=2, C==1L 


There is a third method which is often useful. Suppose that we multiply both 
sides of expression (1.13) by x — 1, to obtain 


2(2 + Ieee Bix—1)  C(xz-1) 
(x + 2)(a@ + 3) E+2 r+3 / 
2x2x® 


and then put x = 1. We obtain A = = 3 immediately. Similarly, 


3x4 
multiplying (1.13) by z + 2 and putting « = —2, and multiplying by x + 3 and 
putting x = —3, gives ‘B and C at once. 


(c) You may have noticed that in the given expression, the numerator is not of 


lower degree than the denominator, so our first step must be to overcome this 
difficulty. This can be achieved by writing 
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e+e (aa —i)+a) +2 
ge xg? —1 
ee = Ti ae 
= a 
* (2 —1)(x +1)’ 
Proceeding as before, we find 
x? + Qa 3 i 3 
g—1 2a—1) 2a2+1) 
Solution 1.8 
@) = = -= 
g2—5¢@+6 ¢2£-3 «2-2 
(b) If 
2 
x a 
es 1) Q(x) 
then 
P@) .- £ Se De 2 1 
Q(z) (x-1) ~ (x —1)? (o—1)? - s-) (= 2?’ 
sO 
- a 2 A 1 
(x — 1)? e-if e ae 
(c) (x — 1)? 4 1 
(x —3)2(2-2) (a-—3)2 a-2 
(a) x —2 ee: 1 = 1 
(Q+x2\(1+2)2 3(2+22) 3(1+2) (142)? 
Solution 1.9 
at+bi at+tbhi c-—d 
= = — bd bc — ad)1 
c+tdi ctdi c-—di 2 ale ee) 
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Solution 1.10 

z = —16+12i (and notice that it would be equally acceptable to write z = 127 — 16). 
lz] = ./(—16)2 + 12? = 4,/42 + 32 = 20 and Im(z) = 12 (not 121). 

z* = —16 — 12% and zz* = |z|* = 400. 


Solution I.11 


From de Moivre’s theorem we have (cos @ + isin 0)? = cos5@+ isin 50. So, from the 
binomial theorem and using the fact that cos? 9+ sin* 6 = 1, we have 


cos 50 = Re((cos 6 + isin @)”) 
= cos’ 6 — 10cos® @sin? 6 + 5cos @sin* 6 
= cos 6 [cos* @ — 10 cos” 0(1 — cos* 6) + 5(1 — cos* 6)*| 
= cos @ (16 cos* 6 — 20 cos* 6 +5). 


Solution 1.12 
(a) First we calculate |z| = 2? + 22 = 2vV/2, then if we plot the point z = 2 — 2% 


on an Argand diagram, we can easily see that Arg(z) = —{. Thus, in polar 


form, we have 
z = 2V2 [cos(—4) + isin(—4))] . 


Any value of the argument of the form —4 + 2k7, where k is an integer, is 
equally acceptable. 


(b) The general polar form is z = 2/2 (cos(—4 + 2km) + isin(—% + 2km)), for k 
an integer. Hence, using de Moivre’s theorem, we obtain 
zi/3 — gi/6 (cos ($(—4% + 2km)) + isin (¢(-—% + 2km))). 


For k = 0,1, 2, this gives three distinct cube roots: 
zo = v2 (cos(—4) + isin(—Z)), 
{= V2 (cos & + isin 4) 
to = V2 (cos + isin 27) 


The remaining values of k simply repeat these roots. 


Solution 1.13 
For any real number n, 
nm _ (re??)” — 7? C7.) 


z ee’ =r"(cosné +isinné). 


Solution 1.14 
We have 
la + bl? = (a+b)? = a? +B + 2ab < al? + |b]? + 2lal|b| = (lal + [d)), 


and the required result follows. 


Solution 1.15 


= 


d d az —2x ae 
a, (cosh 2) = — (<=) ae ——— = sinha. 
dp 
If @ = Asinhwx + Bcoshwz, then — Aw coshwa + Bw sinhwz, so 
£ 
d*o 


—_ — Ay sinhwx + Bu* coshwa = w*¢. 
dx? 


d d oe Ey —— 
4 (sinhe) = * (“as ca 
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Solution 1.16 


1 i 1 
We can write ————~ = — — , therefore 
nin+1) n n+l 
ae n n+l 


This follows because each of the negative terms, except the last, cancels with the 
first term in the next bracket. Thus 


00 1 N : i 
es Pay tocat i) O-555) oo 
Solution 1.17 
We know that 
nies. Se. 4 ee eee 


SO 


g)2n+1 
sinh(2z7) = ae ad 


Solution 1.18 


From the Taylor series we have 


: 1 
sing = x — a + O(x®) and —— tie =O 
Hence 
sin x 


= [a — 22° + O(x?)| [1+ x + O(2*)| 
=x—t4a°+2°4O0(2”) 

=x“+ 2° + O(z°). 

Solution 1.19 

We have 


f(—2) 


oe 3 = 
1+(-2)? = 14+a42 — 


so the function is odd. Thus 


/ f(x) dx = / nee / f(x) de 
= [sear — f H(-2) ae 


1 1 

= / f(x) dx — / f(x) dx Here we have substituted x 
0 0 

= () 


for —x in the second 
integral. 


For any odd function f(x), and any real number a, it is obvious from the graph 
a 


of f that / flerde =. 


a 


Solution 1.20 


iw= sinh’ x, then « = sinhw. Differentiating with respect to x, we obtain 


dw 
t= hw—. 
Cos we 
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It follows that 
dw if 1 1 


1 
dx  coshw — \/1 + sinh? w — Jl 422 


The graph of y = sinha can be obtained easily by reflecting the graph of 
y = sinh in the line y = z. 


Figure 1.11 The graph of sinh7! x 


We have V14+ 22 > Vx? > —a@ for alla EC R,sor+V14+22>0. 


Now 
= (In (x + 1 +22) = ert 


So we see that = Ee - (In (a+ V1+ a), therefore w and In (a +Jf/1+ o*) must 


differ by a constant, giving 
sinh ‘a2 =In (« +1 +2?) +. 


However, sinh~* x = In (x +vV¥1+ o*) =0 when «= 0, so ¢ = 0, and. this estab- 
lishes the required result. 


Solution |.21 


We can write 
4r* —9r +4 2 1 1 
a 
aix—l1\(a-2) a2 «x-1 2-2 


so for x > 2 we have 


ee a ee 
[wep *-/ G+ t+)e 
= 2Inxz+In(# — 1) + In(a — 2) +e 
= In (a*(x — 1)(x — 2)) +. 


It follows that 


[ Ar? —9r +4 
3 


x(a — 1)(a — 2) dx = [In (x e- le- 2))\3 = In(16/3). 


Solution 1.22 


For a £ 1 we have 


T 
i 1 T i 
— di = ee eS i i ae | 
/ t@ 4. F = ) 


i & 


so lim —dt= 
Too | {* C= 


j if a > 1, while for a < 1 the limit does not exist. 
We also showed in a previous example that / 3 dt is divergent. 
1 


sage 
It follows that / a dt is convergent for a > 1 and divergent for a < 1. 
1 
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Solution 1.23 
First we notice that 0 < f(t) + |f(t)| < 2|f(£)|, so 


fi T aes 
/ (FO) +s) at <2 | slat <2 f | f (t)| dt. 
r 
Thus the integral / ( f(t) +\f (t)|) dt is non-decreasing as JT increases, and it is 


a 
bounded above, so the limit 


. 
tim [ (#0 +f) at 
E ‘i 
exists. But we are given that jim | f(t)| dt exists, so 


a 


Jim FE f(t)dt = tim | (F(0) +s) ae f sei 


CO 
exists. In other words, / f(t) dt is convergent. 


int sem | sint 
Since |---| < < =S = Cont # (0) and / ro) dt is convergent, it follows that / a dt 
1 

eat 
is convergent, therefore / — dt is also convergent. 

1 
Solution 1.24 
(a) i-1 (b) 3+% (1 
Solution 1.25 

1 — 

(a). (Se a The polar forms of these 

a is i complex numbers can be 


obtained by using Euler’s 


(b) (1+ i473)? = (2ein/3)” — Aerin/3 


(c) (1 4 iv/3)'® ‘iad (2eir/3)"° ae 915 -15in/3 = HIS {sit — 915 pin 


formula e?? = cos 6 + isin 6. 


Solution 1.26 
e+e” (cos@+isin§) + (cosé — isin 6) 


h ) ES UM a 
cosh(70) 5 5 cos 0, 


ahah = e = _ (cos 6 + isin 6) . (cos 8 — isin 6) whole 


Solution 1.27 

jz + w|? +|z- ol = (z+w)(z+w)*+ (z-—w)(z-w)* 
= (2+ w)(2" + w*) + (2 w)(2* - wv") 
= 22z* + 2ww” 
= 2|z|? + 2|w|? 


Solution |.28 


(a) (i) If y =e, then Iny = az, and differentiating with respect to x gives 
a=y'/y. Substituting y’/y for a, and using the fact that y #0, we 
obtain the first-order equation 

ylny = zy’. 
(ii) Differentiating gives y’ = cos(x + a); using the identity cos? 9+sin? 6 = 1, 


we have 


(y')* +y? =1. 
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(b) 


Since (x + b)y = x +a, differentiating with respect to x gives (a + b)y’ +y = 1. 
Differentiating again, we have (a + b)y’’ + 2y’ = 0. Eliminating z + b from these 
two equations gives 


2(y')* = (y—1)y". 


d 
Integrating the equation a = 1 twice, we obtain first — — 2+ A, then 
x - 


po Su" + Ax + B, where A and B are arbitrary constants. 
| a . , dy 2 ay . 
Differentiating the equation xy + cos y = 3 gives tes + y—sin ue = 0, which 
x x 
dy y 


can be rearranged to: give —— = ————.. 
dx siny-—2Z 


Solution 1.29 


(a) 


The equation can be written as 
ay i+ . 
dx 1—#’ 


so it is homogeneous. 


Putting y = xv as before, and eliminating y, we obtain the equation 


which can be rearranged to give 


/ : J dy = = 
l+q 149 we x 


Integrating gives 


tan-'v — $In(1 +v*)=Ing+A 


for z > 0, and replacing v by y/zx gives 


. 2 
2tan* (2) —Iln € a ) = Jing 4-6, 
x bs 


that is, 
2tan—'(y/x) — In(x* + y*) = B. 


The equation can be written as 


dy 1fz : f l/y 
dx 3\y ee 
so it is homogeneous. 


Putting y = xv, the equation becomes 


cere 
dx 3 \ vy? 
which can be rearranged to give 
ay t= Qu? 


ee By? 


Rearranging again and integrating, we have 


/ 3u2 cass dx 
_.-": Saree e 


hence $ In(1 — 2v) = —Ina+ 5 In A ( for x > 0, 2u° < 1), which can be rewrit- 
ten more simply as 1 — 2vu? = A/z?, so the required solution is 


g? — Qy® = Az. 
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The choice of the arbitrary 
constant in the form 5 In A 
simplifies the final answer. 


50 


Solution 1.30 


The equation can be written as 


dy. =biyvilagsy 
dg (y/a) | 


so it is homogeneous. Putting y = rv, we have 


I +-7= 


dv _ ~haete a 
dx Vv 


which can be rearranged to give 


dv —v*-1l+tvVJ14+v? 
a 


dx 1) 


and therefore 


v dv 


(1 + v?) V1 te 


[o-] eae 
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The integral on the right can be evaluated by making the substitution u = ¥1+v?, 


and we then have 


es 


—wt1 


Integrating, we obtain 


] 


SO 


Ing = —In(—-w+1)+MZC, 
on. ees '# a Ce 
tl-uo t1-V1i+u2 tae— fe? $y? 


Therefore tx — ,/x? + y? = C, which can be rearranged to give 


yg = CTC = 27). 


Solution 1.31 


(a) The equation is first-order linear, and the integrating factor is 


P(e) = exp ( f cot rd) = exp (In(sina)) = sine 


Multiplying the differential equation by the integrating factor gives 


sin x dy + 4 
— cos 2 = — 
de dx 
and integrating gives 


aint = jae sin x dx = —3e°°"" + A. 


Since y = —1 when z = 5, we have (—1)sin 
A = 2, so the required solution is 
ysinxz = —3e°"* + 2. 


In order to solve the differential equation 
dy 9 
—1l)\(a#—2)— +y= 
(e—I)(@—-2)5 + y=2", 


we first rearrange it into the form 


a ee 
de (a—1)(x—2)” (x—1)(a—2) 


wea) =e" * ein», 


= —3e°2 + A, which gives 


Notice that we must ensure 
that the coefficient of 
dy/dzx is 1 before using the 
formula for the integrating 
factor. It is a common 
error to forget to divide 
the equation 


This is of the same form as equation (1.10) on page 35 (i.e. first-order linear), 


with 
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so the integrating factor is 


P(x) = exp ( / g(x) ir) 
-en [(—5-s4)@ 


= exp (In(x# — 2) —In(a —1)) = — 


Now we obtain a further rearrangement of the differential equation by multi- 
plying both sides by P(x) to give 


x—2\ dy 1 om 

ae | pe ee ey 

g=1j de -|#—-1L- (x —1)? 
which can be rewritten in the form 

d {(x-2 sl 2 ‘ | 

dx \z—1°) z-1 (#-1)?’ 


where we have expanded the right-hand side in terms of partial fractions. 
Integrating gives the general solution 


(pe 


y 1 
ey et Aln(z —1)- 7 +6 


which can be rearranged to give 


1 1 
ai =) 
5 («+ n(x — 1) +e) 


XL 
yf = 
tC 


(c) Differentiating v = y'—" with respect to x gives 


dy dv 
1-n)— =y"— 
and the differential equation becomes 
ys dv 
— Fiei=e ; 
2 © +yF (2) =y"G@) 
This can be rearranged into the form 
d 
= + (L=n)y!F(e) = (1=n)G(2), 
which then becomes 
d 
—_ + (1—n)vF(2) = (1—n)G(2), 


which is of linear form (and can be solved by the previous method). 


Solution |.32 


(a) The auxiliary equation m? +m — 6 = 0 has roots m; = 2 and mz = —3, so the 
general solution is y = Ae?” + Be~°*. 


(b) The auxiliary equation m? — 3m + 2 = 0 has roots m; = 1 and m2 = 2, so the 


general solution is y = Ae* + Be?” — ze”. 


Solution 1.33 


The auxiliary equation m? + 5m + 4=0 has roots m; = —1 and mz = —4, so the 
complementary function is y, = Ae~* + Be~*”. 


The form of the right-hand side suggests that we should try a particular integral of 
the form yp = a+ bx. This gives 


tp +5 oe + 4y, = 5b+ 4(a+ bx) = (4a + 5b) + 462 = 3 — 22 
dx? dx : 
Equating coefficients of « and constant terms gives a = a and 6 = —s. Thus the 


general solution is 


Y =Yct Yp = Ae * + Be** + 4 — 5a. 
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Solution |.34 


The auxiliary equation m? + 4 = 0 has roots m, = 27 and mz = —2i, so the general 
solution of the associated homogeneous equation is y.(#) = Acos(2z) + Bsin(2z). 


We now try a particular integral of the form yp = ax + b, which gives 4(az + b) = 2a, 
from which we see that a = 5 and b = 0. Thus the general solution is 


y(xz) = Acos(2x) + Bsin(2x) + $a. 


Solution 1.35 


The auxiliary equation 2m? + 2m +3 = 0 has complex roots m; = —4 + 1,/5i and 
y Z¢ 2 
i= —+ _ EV Di, so the general solution of the associated homogeneous equation 


is 
y_(x) =e %/? (4 cos(4.V5z) + Bsin(3V5e) ; 

We now try a particular integral of the form yp = ax? + br +c, which gives 
4a + 2(2ax + b) + 3(ax? + ba +c) = 27 4+ 22-1. 


; a + pee — - 
Equating coefficients of powers of x gives a = 3, b= 5 and c= —5:, so the general 
solution is 


y(xz) =e */? (4 cos(4.V5z) + Bsin(3V5e) + 5x + 2r— 2. 


Solution 1.36 


We consider the auxiliary equation m* + 1 = 0, for which the solutions are the four 
possible values of (—1)!/4. Writing —1 in polar form, 


—1=cos(2k + 1)m+isin(2k + 1)z, 


for any integer k, and using de Moivre’s theorem, we have 


(—1)1/4 = (cos(2k + 1)m + isin(2k + 1)n)'/" 


= cos ($(2k + 1)m) + isin ($(2k + 1)z). 


Choosing k = 0,1,2,3, we obtain the four values 


a 1+% a i= 7% 
m, = cos} +isin = = * mg = cos 7 + isin 3 = ~; 

i +3 = 1—12 
m3 = cos 22 + isin 2% = ——— m4 = cos 4 + isin 4 = 


V2 - 2 
(It is easy to verify that each of these is a solution of the equation m* + 1 = 0.) 


The general solution of the differential equation is therefore 
y= (P cos(a/V/2) + Qsin(x/V2)) er/Vv2 
+ (R cos(a/V/2) + S$ sin(x/V2)) er. 


Solution |.37 


The equation is separable, and we have 


‘ks dy = [@ dx. 


Using the initial condition y(0) = 0, we have 


Yy E 
/ er dg= / e*” dy, 
0 0 


SO [-e-¥]¢ =5 ie be which gives 1 — e~¥ = 3(e?* — 1), thus 


zZ 
y= in : 
3 —e22 
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Solution 1.38 


The equation is separable, and we have 


[cot way = -2 f ede. 


Using the initial condition y(0) = 5, we have 


Yy x 
/ cot ydy = -2 f eue, 
m/2 0 


so |In(sin wl 2 =— aL which gives In(sin y) = —x”, thus 
y = sin‘ (exp(—2’)). 


Solution 1.39 
Dividing the top and bottom of the right-hand side of the equation by ry, we obtain 
a 
dx 3a2/y—2y/x’ 
so the equation is Sia Putting y = rv gives y’ = rv’ + v, and the equa- 
Vv 


du 
tion becomes — + v = ———-~. Rearranging this, we obtain 
dx 3 — 2v? 


oC - , 


which is separable. Therefore, using the initial condition v(1) = 1, we have 


[< dx -/ ee 
dv. 
v(2v2 — 1) 


The integral on the right can be evaluated by writing the integrand in partial 
fractions as 


3 — 2? +f /2 3 


= 5 E," 
v(2v2 — 1) J/2v —1 J2u+1 10) 


thus 
Ine = f = ae v2 —- dv 
1 fan — 4 f2v +1 U 


= [In(v2v —1) + In(v/2041) =3in a 


- (oF) 


So the general solution is y? = 2y* — 2x”. 


Solution 1.40 


The equation can be rewritten as 


d 
# _¥(14m(¥)), 
i. -z £ 
where the right-hand side is clearly a function of y/z. 


dv 
Now put y = zv, so y’ = xv’ + v and the equation becomes z— + v = v(1+Inv), 


dx 


dv 
which reduces to ae =vlnv. This equation is separable and gives 


iene 
vlnuv 


The integral on the right can be evaluated using the eas apail u = Inv, and we 
obtain Inz = In(Inv) + nC. Thus x = Clnv or v = e®* (changing the constilik). 
Finally, replacing v by y/x, we have the solution y = eRe , 
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Solution |.41 


First we solve the quadratic equation 
dy : dy 
— ) —2y— —-xr=0 1.14 
(#4) - ayo (1.14) 
to obtain the two possible values 

dy give +e 

dx £ 


This can be rewritten as 


ae et Vit (%) 


so it is homogeneous. Putting y = xv gives 


d 
¢— + ¢=~ hee 


dx 


SO 


a aT dv 
a VI + v2 
Thus Ing + Inc = +sinh™' v or v = sinh (In(exr)). Hence 
y eln(cax) i e— In(cx) 


Fe a nie, 
x 2 


which gives 
209 == (rue — 1) 
It is interesting to note that equation (1.14) also has the solution y = +22, which 


cannot be obtained from the above solution. This is because it is a non-linear 
differential equation. 


Solution |.42 


The complementary function (i.e. the solution of the homogeneous equation) is 
found by solving the quadratic equation 


m? —7m+12=0. 


The roots of this equation are m, = 3 and m2 = 4, so the complementary function 
is ye = Aer + Be™. 


In order to find a particular integral, we try a solution of the form y = ax? + bx + ¢, 
and when we substitute this into the differential equation and equate coefficients, 


we find that a = 1, b= 2 end ¢= =: 
Thus the required general solution is 


= on Ax 2 7 
y= Ae = Be” 2 ae <5- 


Solution 1.43 


d 
(a) Let be the gradient at a fixed point (x,y) on the curve, with x > 0 and 
x 
y > 0. Then the equation of the tangent line at this point is 


dy 
Y-y=—(X - 
dx ( ), 
where X and Y vary along the line. When X = 0, we have 
dy 
Y=y-2— 
Y a 


and the given condition then becomes 


aie dy 
z (#) 7. =e ee 


Solutions to Exercises in Chapter | 55 


This equation can be rearranged to give 


dy 
x? — y =. ed 


and then 
dy y*—x27 
dx Qry 
This equation is homogeneous, so we substitute y = xv to obtain 
— +VU= oe = 
dx 2u 
Then 
dv 1l+uv 
dx 20 ” 


eft is, 
ite? 


Integrating, we obtain — Inx + In(2C) = In(1 + v7), which gives x? + y? = 2Cx 
or 


OE tyes get os 


(a semi-circle of arbitrary radius C with centre at (C,0)). 


In order to find the orthogonal family of curves, we must solve the equation Two lines with gradients T 
and —1/T will always be 
dy 2xry 
<< — orthogonal. 
dx y? — x? 


Substituting y = xv gives 


Vv 2v 
er 7? 
so 
dv — v(1+v%) 
“de 1—v - 


Integrating gives 


[¢- [= 1 —v’ w= | ee es 
(1 +") a a 
so Ina — In(2A) = Inv —In(1+ v7), which simplifies to = = eee and there- 
fore 
en AP = * 
(a semi-circle of radius A with centre at (0, A)). 


The two families of curves are shown in Figure 1.12. 


Figure 1.12 Each curve of the second family is orthogonal to each member of 
the first family 
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Solution |.44 


Suppose that at time ¢t there are z(t) grams of substance C, which must have been 
formed from 2z(t)/3 grams of A and z(t)/3 grams of B. Then 


dz 22 Z 
a a OO 7) (00 a 


where k is a constant which dictates the reaction rate. This can be rearranged to 
give 
2k 2kt / dz 1 


9 9 c— ies 


Since z(0) = 0, we see that c = —1/300, thus 
2kt 1 1 
“9:  300—< 300° 

We can deduce the reaction rate k from the fact that z(30) = 100, which gives 
60k 1 1 1 


“9 ~~ 200 =300 ~—-<600’ 
so k = 1/4000. It follows that 


3008 


CHAPTER 2 


Multivariable and vector 
calculus 


2.1 Introduction 


As in the previous chapter, our intention here is to provide some of the basic 
tools of applied mathematics, and the emphasis is on methods rather than 
justification and proof. 


Some of the material may be familiar, so you can treat it as revision, but 
there is much that you may find new. The reasons for introducing some of 
the techniques may not be entirely apparent at this stage, but their relevance 
will become evident as the course develops. For the moment you should 
concentrate on acquiring the skills that will be needed later, and attempt as 
many of the exercises as possible. 


Some ideas, such as continuity, are included because they are an essential 
component of every mathematician’s background, but we shall not ask de- 
tailed questions on their use. Other notions, such as partial differentiation, 
will be used throughout the course and you will need to gain some expertise 
before you move on to later chapters. Generally, you will be able to judge 
the importance of the topics by the number of exercises we ask you to at- 
tempt. However, the Gauss divergence theorem is an exception and will be 
of crucial importance later. 


Subsections 2.2.1 and 2.2.2 are mainly background reading. Subsections 2.2.3 
and 2.2.4 discuss various techniques that will be essential later on. 


Subsections 2.2.5 to 2.2.8 discuss material which may be familiar, but they 
lead on to ideas that we shall pursue in Block III. 


Subsection 2.3.1 covers material that should be familiar, and, while the 
remaining subsections discuss techniques that are of universal application, 
they lead on to the Gauss divergence theorem, which is the most important 
result of this chapter. 


Section 2.4 introduces techniques that will be used throughout the course. 
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2.2 Multivariable calculus 


Functions of two variables 


Most of the techniques that we discuss in this section can be extended to 
functions of three or more variables, but we concentrate on functions of 
two variables for two reasons. First, while the complexity of the results 
generally increases with the number of independent variables, functions of 
two variables often provide a useful insight into the more general cases. 
Secondly, functions of two variables are relatively easy to visualise because 
they can often be represented graphically in a three-dimensional coordinate 
system. 


A real function of two real variables z = F(x, y) is a mapping from a set of 
points D in the Cartesian plane to the real numbers. The set D is known 
as the domain of the function, and each point (7, y) of D maps to a unique 
value z. The variables x and y are known as the independent variables, and 
z is the dependent variable. 


Consider the real function F defined by F(x,y) = ,/1—(a?+y?). For Remember that the square 
F (x,y) to be real, we require x? + y? <1, so we take the disc defined by T00t symbol indicates the 


this inequality to be the domain of F’. ee 

Here, and generally, we 
assume that the domain of 
a function is the largest set 
on which the definition is 
sensible. 


Figure 2.1. The hemisphere z = \/1 — (2? + y?) 


If we put z = ,/1 — (a7 + y?) and let (2,y,z) correspond to a point in a We use either (x, y) or 
three-dimensional Cartesian diagram, then the function F can be repre- (Z,y,0) when referring to 


sented as a hemisphere, as shown in Figure 2.1. points that lie in the 
(x, y)-plane. 
Generally, a function w = F(a,y) can be represented as a surface, and 


(x,y, F(x, y)) represents an arbitrary point on the surface. 


2.2.1 Level curves 


It is not always easy to visualise a surface in three dimensions, and there 
is an alternative pictorial representation of a function of two variables that 
is sometimes helpful. On an ordnance survey map you may have seen lines 
connecting points of equal height above sea level. Such lines are known as 
contours, and they give some indication of the shape of the surface. ‘The 
same idea can help us to understand the behaviour of functions. 


If we are given a function z = F(x,y), then we may be able to plot the 
contours or level curves corresponding to the equations F'(x,y) = A for 
various values of the constant A. 


2.2 Multivariable calculus 


For the function F(x, y) = ,/1 — (a? + y*) we consider the curves 


V1—(a2?+y?) =A, where0O<A<1. 
In this case it is convenient to rearrange the equation into the form 
geo te y2 wohge AP, 


and we can then see that the level curves are circles of radius \/1— 7, as 
shown in Figure 2.2. Such a diagram can often give some insight into the 
behaviour of an unfamiliar function. 


Bs 


Figure 2.2 The level curves 2? + y? = 1— 


Exercise 2.1 


Sketch the level curves of the functions F(x, y) = 3 — 24% + y? and G(a,y) = y/2?. 


2.2.2 Continuity 


Perhaps the best way to appreciate the meaning of continuity is to examine a 
function that fails to be continuous; we begin with a function of one variable. 


Example 2.1 
Sketch the graph of the function f(x) = 


jar] 


Solution 


Figure 2.3 f(z) =1ie >Oend fey] -—l ise < 0 


The function is defined for all real values of x except x = 0, with 


l- tc >. 
fe)=4 4 = 2 =< Q. 


The graph is shown in Figure 2.3 with a break at the origin. The essential 
point here is that whatever value we assign to f(0), the resulting function 
will be discontinuous at the origin. Hf 
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You saw in Chapter 1 that we define a function f(a) of one variable to 
be continuous at « =a if lim f(x) = f(a). Clearly, in the above example 
ta 


lim f(x) does not exist because f(x) approaches —1 when x approaches 0 
2 


from below, and f(x) approaches 1 when x approaches 0 from above. 


Now we move on to discuss functions of two variables. 


Example 2.2 
Sketch the surface that represents the function F(x, y) = eal 
LY 


Solution 


Figure 2.4 The function F(2,y) = —> 


We see that F(x, y) = 1 when x and y have the same sign, and F(z, y) = —1 
when they do not. The corresponding surface is shown in Figure 2.4, and 
it is clear that the function behaves badly on the axes. First, the function 
is not defined on the xz- and y-axes, but, even if it were to be defined there, 
the sudden changes in value from —1 to 1 would ensure that the function is 
discontinuous. 


A geometric description of continuity is a useful aid to understanding, but 
it cannot be extended easily to functions of more than two variables. A 
slightly different approach, based on comparing the values of the function 
at neighbouring points in the domain, is more productive. 


For a function of one variable, we define f to be continuous at 7 =a if 
lim f(x) = f(a). In other words, we can make f(x) as close to f(a) as we 
4 bee a 0 | 


please by choosing x sufficiently close to a. It is this notion that we wish to 
extend to functions of two variables. 


We say that a function F is continuous at a point (a,b) if F(x, y) converges 
to F(a, b) as (x,y) approaches (a,b) from any direction. (We shall make this 
statement a little more precise when we discuss functions of many variables 
later.) 


In this course the difficulties that arise from discontinuous functions are not 
our main concern, so you may assume that any given function is continuous 
unless we say otherwise. 


Many of the results that come later in the course are crucially dependent on 
the fact that the relevant functions are continuous, but we shall be avoiding 
these complications. 
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2.2.3 Partial derivatives 


In addition to continuity, we often require that a function should be ‘smooth’, 
which in intuitive terms means that there are no sudden dents or sharp 
corners on the surface that represents the function. ‘This is equivalent to 
the requirement that it is possible to construct a flat surface, known as the 
tangent plane, that touches the original surface at the point in question. 
The existence of such a tangent plane is intimately connected to the existence 
of the partial derivatives of the function, so first we introduce the notion 
of partial derivatives and then we use them to define the tangent plane. 


The partial derivative of a function F'(x,y) with respect to x is obtained 
by differentiating the function with respect to x while keeping the variable 
y fixed, and similarly for the partial derivative with respect to y. These 
partial derivatives are often denoted by 


although alternative notations can be useful, e.g. F,(x,y) and F(x, y), or 
just fF, and Fy. 


OF F 
The partial derivatives — and —, when evaluated at x =a and y = 0, 


Ox Oy 
represent the slope of the surface z = F(x,y) at the point (a,b, F(a,b)) in 


the x and y directions, respectively. 


Example 2.3 
OF OF 
Given F (x,y) = x7 + y°, calculate >. and — and hence deduce the slope 
L Y 
in the z and y directions at the point (3, 2). 
Solution 
OF , : . 
To calculate ae we differentiate with respect to x keeping y constant: 
“ 

OF 

— = 22. 

Ox ss 

OF ; 
To calculate Ay we differentiate with respect to y keeping x constant: 
Y 

OF 

= By 

Oy 

O 
Alternatively, we might have written z = x? + y®, in which case — = 2s 
3 

Oz 
and = oe. 

Oy . 


Hence the slopes in the x and y directions at the point (3,2) are F,,(3, 2) = 6 
and F,,(3,2) = 12, respectively. 


Exercise 2.2 


FP FP 
Given F(z, y) = x7y® + sin(re”), calculate nil and ad 


Ox Oy 
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Example 2.4 


Given an equation such as x? + y? — z* = 1, we need to clarify which are 
the dependent and independent variables. 


(a) Supposing that z is the dependent variable, while x and y are the inde- 
z Oz 
— and — 


Ox Oy 
(b) Alternatively, supposing that y is the dependent variable, while x and 
Oy Oy 


z are independent, evaluate — and —. 
Ox az 


pendent variables, evaluate 


Solution 


(a) Rearranging the equation, we obtain z = +,/x? + y? — 1, which is real 
if x7 + y? > 1. Choosing the positive square root, and differentiating 
first with respect to x and then with respect to y, we obtain 


Oz x Oz y 

a and SS 

de JP ty)—1 “Wy G+ P)=1 
(with similar results if we choose the negative root). 


Alternatively, we may differentiate the equation x7 + y? — 27 =1 di- 
rectly. First we differentiate partially with respect to x, and notice that 
the partial derivatives of x? and y? with respect to x are just 2x and 0, 
respectively, while the partial derivative of z? (a function of x and y) is 


O 
222° (keeping y fixed and using the chain rule). Thus we obtain 
. 


O 
Oz 
Be a = E 
: Ox 
Similarly, differentiating partially with respect to y, we obtain 
Oz 
2y — 2z— = 0. 
Y Ay 
Rearranging these equations, we have 
O O 
ae = where z # 0. 
on 2 De. z 


b) Differentiating the equation x? + y? — z* = 1 partially with respect to 
- 


x and z now gives 
O O 
2% + ayo =() snd y= = 22 =i), 
Ox Oz 
and rearranging gives 


wherey #0. 


Example 2.5 


d 1 dz 1 9 1 
w that — ee ee in Seen ee 
se dt (| cP) | A tess ke 


Solution 
We have 


£({ te 347 ee 
dt\Jo (x+t)?/ dt| x+t}], dt\t 14+t/ (+t)? # 


1 


[a (ere)@-L aren ler], # 


Note that we use partial 
derivatives (0) to 
differentiate the integrand, 
and total derivatives (d) to 
differentiate the integral. 
This is because the 
integrand is a function of x 
and t, while the integral is 
just a function of t. 
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Generally, we shall assume that we are justified in ‘differentiating through 
the integral sign’, but if necessary we may appeal to the following result. 


Leibniz’s rule 


Let F'(x,t) be a continuous function with a continuous partial derivative 
OF /ot on the rectangular region a < x < b, t1 <t < tg of the (2, t)- 


plane. Then, for t; < t < tg, 


This result will be used in Block III. 


Exercise 2.3 


1 1 
Evaluate / sin(ax) dx, then use Leibniz’s rule to evaluate / x cos(ax) dx. 
0 0 


(It is not necessary to verify the conditions for Leibniz’s rule.) 


2.2.4 Linear approximation 


It is often useful to obtain an approximation to a given function near a 
given point. For a function of one variable, y = f(a) Say, we may use a 
linear approximation f(x) ~ f(a) + f’(a)(a — a), which is equivalent to ap- 
proximating the function by the tangent line at x = a. 


For a function of two variables we can obtain a similar linear @timate in 
terms of the partial derivatives, and in geometric terms this estimatacan be 
interpreted as approximation by the tangent plane. 


The essential idea is that we find a plane that touches the correspondin® 
surface at a given point, and then we approximate values on the surface by 
values on this nearby plane. ‘To make the ideas more precise, we need the 
equation of a plane. 


(X.Y F ey) Pi 


et 2 
ull 


Figure 2.5 ‘The surface representing F’ near the point P 


In Figure 2.5 the surface corresponding to the function z = F(x,y) passes 
through the point P with coordinates (a,b, F(a,b)), and our aim is to esti- 
mate the value of F(x, y) when (x,y) is close to (a,b). 
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The equation of a plane 


The equation z = A+ B(x —a)+C(y—b) represents a plane passing 
through the point (a,b, A). Differentiating partially with respect to x and 
then with respect to y, we obtain 


Oz Oz 
aoe and re 


In other words, the slope of the plane in the x direction is B, while the slope 
in the y direction is C. 


The tangent plane to the surface z = F (x,y) at the point (a,b, F(a,b)) is 
defined as having the same slope as the surface in the x and y directions, so 
its equation is 
OF OF 
z= F(a,b)+ (4# —a)— + (y—- 0), 
(a,b) + (2 — a) 5 + (yb) 

where the partial derivatives are evaluated at (a,b). Thus our linear approx- 
imation for the given function is 


F(z,y) ~ F(a,b) + (@ — a) + (9S i) 


This result is often written in terms of the delta notation for small incre- 
ments, in which case it becomes 


6 x a + Ge (2.1) 
Ox Oy 


where 6F = F(z,y) — F(a,b), 6c =x—a and dy=y-—b. Equation (2.1) 
tells us approximately how much the function F’ changes when we move 
from a point (a,b) to a neighbouring point (a+ dz,b+ dy). This approximate 
change in the function is depicted in the inset of Figure 2.5. 


2.2.5 Local extrema 


You will probably recall that the equation (x — a)? + (y — b)? = & defines a 
circle of radius 6 with its centre at (a,b). The interior of this circle is defined 
by the inequality (x — a)? + (y — b)? < &. 

A point (a,b) is said to be a local maximum of a function F(x, y) if there 
is a region (x — a)? + (y — b)? < 6? such that F(z, y) < F(a, b) at all points 
(x,y) in the region. 

A point (a,b) is said to be a local minimum of a function F(z, y) if there 
is a region (x — a)? + (y — b)? < 6? such that F(z, y) > F(a, b) at all points 
(x,y) in the region. 


Local maxima and minima are known collectively as local extrema. 


Stationary points 


A point (a,b) is a stationary point of the function F (x,y) if the partial 
= FP 

derivatives a and ai are both zero when evaluated at (a,b); in other 
x y 

words, F;,(a,b) = F,(a,b) = 0. 


Stationary points of functions of two variables play much the same role as 
such points for functions of one variable. This stems from the fact that it is 
possible to show that (for a well-behaved function) every extremum must be 
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a stationary point. Thus, if we wish to find a local maximum or minimum 
of a given function, we may restrict our search to its stationary points. On 
the other hand, not all stationary points are local maxima or minima. A 
stationary point that is not an extremum is called a saddle point. 


Figure 2.6 shows examples of typical stationary points. 


local maximum local minimum saddle points 


Figure 2.6 ‘Typical behaviour of a function near a stationary point 


Exercise 2.4 


OF | 
(a) Find the partial derivatives —— and —— of the function 


Ox Oy 
F(a, y) = e"¥ + x” +siny. 


(b) Given /a+./y+ \/z = 1, and regarding z as the dependent variable, calculate 


(d) (i) Show that the surface z = F(a, y) = 1+ (x — 2y+ 1)(y — 2x + 1) passes 
through the point (1,2,—1), and find the tangent plane at that point. 


(ii) Show that F(z, y) has a stationary point at (1,1). Show that the tangent 
plane at the point (1, 1,1) is parallel to the (x, y)-plane. Putting s = x — 1 
and t = y — 1, or otherwise, find an expression for F(x,y) — F(1,1) in 
terms of powers of x — 1 and y — 1. 


(iii) The lines y = 2— 2 and y = z intersect at the point (1,1). What can you 
say about the sign of F(z,y) — F(1,1) on the line y = 2 — a and on the 
line y = x? What do you conclude about the stationary point at (1,1)? 


(e) Show that the function F(x, y) = (x? + y’)e” has two stationary points, and 
that the stationary point at the origin is an extremum. 


2.2.6 Higher-order partial derivatives 


From Exercise 2.4, you may have observed that the partial derivatives of a 
function F'(x,y) are themselves functions of two variables. This means that 
they too can be differentiated partially with respect to x and y. 
For example, if F(x,y) = x* sin y, then 

OF OF 9 


— =2rsiny and ——=2* cosy. 


Ox Oy 
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Differentiating again, 


= (ee) = Z8ny 5 (Fe) = 2x cosy 

Ox \ Ox ’ Oy \ Ox 

“ (=) = 26 eosy and = (=) = —z’ sin y. 
Ox \ Oy Oy \ Oy 


These are known as the second-order partial derivatives or second partial 
derivatives of the function F(2,y) and, respectively, they are commonly 
abbreviated to 


OF - OF -- OF 0? F 
Ox?’ = =OyOx’ Ozxdy Oy?” 
or alternatively Fox, Fy2, Fay and Fyy. 
OF 
It is no accident that in this example we have a = ——; this is generally 
OyOxr OxOy 


true provided that the function F(x, y) is sufficiently smooth for the second- 
order partial derivatives to exist and to be continuous. (This result is often 
known as the mixed derivative theorem; we shall not prove it here.) 


It is, of course, possible to take still higher derivatives, e.g. Prax, Fyaz- 
The mixed derivative theorem then states that the order in which these 
derivatives occur does not matter; so, for example, Fycen = Eager = Pooys = 


Ff tio: 


Taylor approximations 


It is often the case that the linear approximation to a function is not suf- 
ficiently accurate to provide the information that we require; we may need 
a higher-order approximation. Suppose that near (0,0) it is possible to 
approximate the function F(x, y) by a quadratic expression so that 


F(z,y) ~ A+ (Br+Cy) + (Px? + 2Qzy + Ry’), 
where A, B, C, P, Q and R are constants. 
This is known as a quadratic approximation for F(x, y) near (0,0). 
More precisely, we suppose that 
F(x,y) = A+ (Ba+ Cy) + (Px? + 2Qzy + Ry’) 
+ (higher powers of x and y). 
Putting x = y = 0 on both sides of this equation gives A = F'(0,0). 
Differentiating both sides of the equation partially with respect to x gives 
ad = B+ (2Px + 2Qy) + (higher-order terms); 
then putting 7 = y = 0 gives B= F,,(0,0). Similarly, differentiating par- 
tially with respect to y and putting z = y = 0 gives C = F, (0,0). 


Differentiating both sides of the approximation twice with respect to x gives 
O° F 


Ox? 
then putting 7 = y = 0 gives P = + Fre (0, 0). Similarly, Q = 5 Fry (0, 0) and 
R = 4Fy,(0,0). 

Thus, for a function F (x,y) that is sufficiently smooth near (0,0), the 
quadratic approximation to F'(x,y) near (0,0) becomes 


F(x,y) = F(0,0) + Fy(0,0)a + F,(0, 0)y 
+ 5 (Fra (0,0)2* + 2Fyy(0,0)xy + Fyy(0,0)y7) . 


2P + (higher-order terms); 


(2.2) 
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This is known as the Taylor approximation of order two for F(x, y) 
near (0,0). 
(Notice that the linear terms in this expression are identical to the linear 


approximation that we discussed earlier. ) 


It is possible to extend this process and obtain an infinite series for F(x,y) In this course we shall not 
in ascending powers of x and y. Such a series is known as the Taylor series generally require 


for a function of two variables. approximations of order 
higher than two. 


Example 2.6 


Find the Taylor approximation of order two for F(a,y) = (a+ 2y + lje”4 
near (0,0). 


Solution 
We have 
Fe =e +(2@+2y+1)ye% and F, = 2e°% + (x + 2y+ I)ze™. 
Then 
Fire = (2+ y + zy + 2y")ye™, 
Foy= (1 + 2e + 49d a?y + cy + 2ry*)e™, 
Fy = (4+ a4 2xy + 27)ze™. 
Thus F(0,0) = 1, F,(0,0) = 1, F,(0,0) = 2, Fre(0,0) = 0, Fyy(0,0) = 1 and 
F,,(0,0) = 0. Substituting these values into the general expression (2.2), the 
Taylor approximation is 
(x + Qy + 1)e™ ~ 1+ (x + 2y) + $(2zy) 
=1l+2+2y+ zy. 
Alternatively, we could obtain the same result more easily by noticing that 
eV =1+ Hay) + A(cy)?+---, 90 
(xa + 2y+ 1)e™ ~ (x@+ 2y +1) (l+ay+ (ry)? +--:) 
—l+2+2y+2y+ (higher powers of x andy). 


A similar argument applies to approximations near an arbitrary point (a, b). 
For a function F(x,y) that is sufficiently smooth near (a,b), the Taylor 
approximation of order two for F'(x,y) near (a,b) is 
F(x,y) ~ F(a,b) + F(a, 6)(@ — a) + Fy(a, b)(y — ) 
+ 5 (Fea (a, b)(x — a)” + 2Fery(a,b)(x — a)(y — 6) 
+ Fyy(a, b)(y — b)*) 

(Again, the linear terms in this expression are identical to the linear approx- 
imation to F(x, y) near (a, b).) 


The corresponding infinite series in powers of x — a and y — b is known as 
the Taylor series of F' at (a,b). 


We shall not require terms of the Taylor series beyond those corresponding 
to the Taylor approximation of order two. However, if we let dv = x — a, 


dy =y—), po = (a,b), dp = (6a, dy), and define an operator U = baa _ we 
x 
then we can write the Taylor series at (a,b) in the neat form . 


os | 
F (po + 6p) = F(po) + | GU" 
k=1 


(where all the partial derivatives are evaluated at po). The first three terms of 
this series give the Taylor approximation of order two for F near (a,b). 
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The classification of stationary points 


Suppose that F'(a2,y) has a stationary point at the origin, so that the first- 
order partial derivatives are zero at (0,0). Then the Taylor approximation 
of order two for F(x, y) near (0,0) becomes 

F(a,y) ~ F(0,0) + 5 (Fee (0,0)x* + 2Fyy(0,0)ry + Fyy(0,0)y") . 


Near an extremum at (0,0), it is the sign of F(x, y) — F(0,0) that determines 
if the point is a local maximum or a local minimum, so we are led to consider 
the sign of expressions of the form 

Pa? + 2Qazy + Ry’ 
(where the constants P, Q and R are not all zero). It turns out that the 
sign of this expression is dependent on the sign of PR — Q?. 
Suppose that PR — Q? > 0. Then P cannot be zero, so we may write 
(Px + Qy)’ + (PR- Q*)y? 

P y) 


where the numerator is positive (except at the origin). 


Px? + 2Qzy + Ry’ = 


Thus, at points other than the origin, Px? + 2Qxry + Ry’ takes the same 
sign as P. 


Now putting P = Fy,(0,0), Q = Fzy(0,0) and R = Fy,(0,0), we see that 
the function F(a, y) has an extremum at the origin if 


F,(0,0) = F,(0,0) =0 and PR-—Q?*>0, 


and this extremum is a local maximum or minimum according as P (or R) 
is negative or positive, respectively. 


The above result can be easily extended to a function F(x,y) that has a 
stationary point at a point (a,b) other than the origin. In this case we need 
only observe that the function F(z +a,y+ 6) has a stationary point at the 
origin, and, applying the previous result, we obtain the following. 


A test for extrema 

The function F(z,y) has an extremum at the point (a, b) if 
F,(a,b) = F,(a,b)=0 and PR-Q*>0, 

where P = F,,(a,b); Q = Fata oma = 7.) 


This extremum is a local maximum or minimum according as P (or R) 
is negative or positive, respectively. 


If PR — Q? <0, it can be shown that F(z,y) does not have an ex- 
tremum at the point (a, }). 


A stationary point that is not an extremum is known as a saddle point. 


Near an extremum at (a,b), F(x, y) — F(a,b) does not change sign, but 
near a saddle point there are always points for which F(z, y) — F(a, b) 
is positive and points for which it is negative. 


If PR — Q? = 0, then the above test for an extremum fails. The point in 
question may or may not be an extremum, and this may be determined by 
examining a higher-order approximation, although we shall not discuss the 
general process here. 
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Example 2.7 
Find and classify the stationary points of the function w = x + y? + 32y. 
Solution 
We have 

O O 

Se & Ge? 4 By and mee ee 

Ox Oy 
so at a stationary point, x? ++y = 0 and y* +2 = 0. Solving this pair of equa- 
tions, we obtain z = 0 and y = 0, or x = —1 and y = —1, giving stationary 
points at (0,0) and (—1,—1). We also have 

Oew Oew Oew 

—~ = 62, =g d —> =6y, 

a . Oxdy sy Oy? d 
SO , 

O?w 0*w O7w 
PR-Q? = —= | = 36ry — 9. 
. Ox? Oy? (=) sis 


At (0,0) we have PR — Q? = —9 < 0, so the stationary point is a saddle 

point.-At (—1, —1) we have PR — Q? = 27 > 0, so the stationary point is an 
2 

ali = 6x evaluated at (—1,-1), 

| i 

giving P = —6 < 0, so the extremum is a local maximum. & 


extremum. Also at this point we have P = 


In the functions in the following example, PR — Q? = 0, so the test fails. 
However, in these cases, we can use elementary analysis to determine the 
nature of the stationary point. 


Example 2.8 


(a) Given F (x,y) = x+y’, show that PR — Q* = 0 and that the function 
has a local minimum at (0,0). 


(b) Given F(a, y) = x°y?, show that PR — Q? = 0 and that the function 
has a saddle point at (0,0). 


Solution 
(a) Since F,(0,0) = F,(0,0) = 0, (0,0) is a stationary point. We have 
ot Pw _ ( dw 
Ox? Oy? OxOy 
so at (0,0) we have PR—Q* =0. However, F(z,y) = oe >) = 
F(0,0), so there is a local minimum at (0,0). 


Since F,,(0,0) = F,(0,0) = 0, (0,0) is a stationary point. We have 


0? w 02w ow \* 
at yt ~ (any) ~ 6a ME=y) — Ba?) 


so at (0,0) we have PR — Q? = 0. However, F(a, y) — F(0,0) = ay?, 
which is positive when x and y have the same sign, and negative when 


they have opposite signs. Thus there cannot be a local extremum at the 
origin, and the point is therefore a saddle point. 


2 
) = (2s?y!)(azety?) — (a6%y")? 


Sg 


Exercise 2.5 


a) Find the Taylor approximation of order two to the function w = e*? near (1, 2). 
Find the Tayl imati f order two to the functi _ 12 


(b) Find and classify the stationary points of the function w = sin x cosh y. 
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Later we shall discuss partial differential equations, and this is an opportune Recall that a partial 


moment to introduce a solution of such an equation. differential equation is any 
equation involving partial 


derivatives. 


Exercise 2.6 


Show that any function of the form F (x,y) = f(x) + (y), where f and g are 
0°F 
0. 


differentiable, satisfies the equation 


OxOy 


2.2.7 Functions of many variables 


Suppose that we wish to investigate the behaviour of a function F of the n 


variables £1, 2%2,...,%n-. It is often convenient to write x = (21, 22,...,2n) 

and to regard F’ as a function of the vector a. We shall discuss vectors in 
the second part of this 

Let x = (%1,22,...,2%n) and a = (aj, 42,...,@,). Then we define a func- chapter. 

tion F(a) to be continuous at a if lim F(x) = F(a), and note that x may 


approach a from any direction. 


The linear approximation is easily extended to functions of many variables. 


If a = (a1, a2,...,An), we write 6F = F(a) — F(a) and 6x; = 2; — a; for 
1=1,2,...,n. We then have 
OF OF OF 
OF ~ —6d — 6 +--+ ——dEn, 
Ox1 ih Ox2 a OLn s 


where all the partial derivatives are evaluated at a. 


It is also possible to extend ‘Taylor approximations of higher order to func- 
tions of many variables, although, as one might expect, the complexity of 
the expressions increases with the number of variables and with the order 
of the approximation. 


Exercise 2.7 
Given x = (21, 22,23) and 
F(a) = A+ (Bix, + Boxe + B3x3) 
+ (Cy 24 + Coon + C3322 + Cyox122 + Co3r0x%3 + C312%321), 


calculate the coefficients in this expression in terms of partial derivatives of F’. 


Taylor approximation of order two 


For an arbitrary function F(a) of three variables, the Taylor approximation Exercise 2.7 illustrates 


of order two near x = a is = = es girionis en 
fia} gh +e sere + (%2— a2) a + (x3 — 23) a 
- ; (x: ~ an) Sa (xq — on) + (x3 — oa) 
+ 2(x21 — a1) (x2 — an) gor 
+ 2(x2 — a2)(x3 — oa) 


O° F 
i 2(x3 ae a3)(x4 anit a) | , 
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where all the partial derivatives are evaluated at x =a. It is easy to see 
how this result can be extended to more than three variables. 


Stationary points 


The notions of stationary points and extrema can be easily extended to 
functions of more than two variables. 


Near the point a = (aj,@2,...,@n), it is the nature of the difference 
6F = F(x) — F(a) that characterises the behaviour of the function F’. It 
this difference is positive (or zero) for all x close to a, then the function will 
have a local minimum at a; alternatively, if the difference is negative (or 
zero), then the function will have a local maximum at a. However, if in every 
neighbourhood of a there are points at which dF’ = F(a) — F(a) is positive 
and points at which it is negative, then there cannot be an extremum at a. 


There is an alternative approach to stationary points which will be relevant 
in later chapters. The idea is to examine the behaviour of a function as we 
approach a stationary point from a nearby point. In Block III we shall see 
that this alternative approach can be extended to a situation in which there 
is no direct analogue of the notion of a partial derivative. 


Since we are concerned only with points close to a, it is convenient to write 
5 a+ €&, where E= ee with G+ Gt +E = - and € 
is small. (Essentially, this means that we approach the point a along the 
straight line joining a to a+ €.) We expect dF to be small when ¢ is small. 
We might expect that 6F = F(a+e€) — F(a) is approximately proportional 
to €, or, more specifically, that 


F(a + €&) — F(a) 


is bounded (i.e. not infinite) when ¢ is small. For well-behaved functions, 
and for an arbitrary point a, this is generally the case, and we have 


OF ae: 


OF 
F(a+<c&)—F(a)~e (sé + Dae °? 5 elie o-én) = OTe}. 


Example 2.9 
Given F(x) = 2129 + 273 + £32) and a = (1, —1, 2), find F(a+ e€) — F(a) 
and show that F(a + e€) — F(a) = O(e) for all €. 
Solution 
We have 
F(a + e€) = (1+ e€,)(—1 + e€2) + (—1 + e€5)(2 + €€5) 
+ (2 + €€3)(1 + €&1) 
ss sh el E par fo rieahi te tes 4 ab) 
+ €*(E& + bo€s + £361) 
= —1+e(& + 3€,) + €7(E1&> + £03 + £361) 
and F(a) = —1, so 
F(a + e€) — F(a) = €(&; + 3&2) + €° (6189 + £263 + £361) = Ole). 
The previous example illustrates what we expect to happen for an arbitrarily 


chosen point a. If ¢ is small, then €? is smaller, so the difference F'(a + e€) — 
F(a) is approximately proportional to e. 
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However, for a given function F’, there may be exceptional points a for which 
F(a+ e€) — F(a) = O(e7) for all € 

(in other words, dF is smaller than we might expect). These are, of course, 

the points for which 
Co a —, oF 
Or, Oro ait OLn 


so that the linear approximation for 6F' is zero. These exceptional points 
are in fact the stationary points of F’. 


0, 


Example 2.10 


Given F(a) = 71%29 + x9%3 + £321 and a = (0,0,0), find F(a + c€) — F(a) 
and show that F(a + e€) — F(a) = O(e*), so F has a stationary point at 
the origin. 


Solution 
In this case we have 
F(a + e€) — F(a) = €7(€,€) + €9€3 + &3€) = O(e7), 
and we conclude that the point a = (0,0,0) is a stationary point. 


We could, of course, have reached the same conclusion by calculating the 
partial derivatives 


= 7 and = ee es 
— =%2+ 2 = a oe 8 ee : 
OF e oF 
We Can S66 that. = = = = — =0 at (0,0,0), which confirms that it 


Ox4 Ox2 Ox3 


is a stationary point. 


We could, therefore, have defined the stationary points of a function F' to 
be those points a for which F(a + ce€) — F(a) = O(e7), and indeed it is this 
formulation that is most easily extended to the cases we shall consider later 
in the course. 


Generally, the classification of stationary points (into maxima, minima and 
saddle points) involves the determination of the sign of a quadratic expres- 
sion in many variables, but we shall not discuss that process for functions 
of more than two variables. 


2.2.8 The chain rule 


Often we encounter situations where we wish to transfer between coordinate 
systems. For example, we may find it convenient to transfer from Cartesian 
coordinates (x,y) to polar coordinates (r,@), in which case we may wish 
to represent partial derivatives with respect to x and y in terms of partial 
derivatives with respect to r and @ (or vice versa). The device which makes 
this possible is known as the chain rule. 


Suppose that we are given a variable w which is dependent on two variables, 
x and y. Then the linear approximation gives 


Ow Ow 
dw ~ —dxr + —dy, Zon 
An By oY (2.3) 
where the accuracy of the approximation increases as 6x and dy both ap- 
proach zero. But suppose now that the variables x and y are themselves 


dependent on a parameter, t say. For example, we could be given w = e7Y 
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where x = cost and y = sint, and wish to examine the behaviour of w as 
t varies. (In this case the point (x,y) would lie on a circle, so we would 
examine the behaviour of w = e*Y for points on this circle.) 


Suppose that dt is the increment in ¢ that gives rise to the increments 0x 
and dy. Then (dividing both sides of equation (2.3) by dt) 


bw Owdx | Ow dy 

ét ~ <Ox bt | Oy bt’ 
and in the limit as 0t approaches zero, the approximation becomes equality 
and we have 

dw Owdxr Owdy 

dt ~~ Ox dt ° Oy dt’ 


This result is known as the chain rule for functions of two variables. 
(It is usually abbreviated to the chain rule when there is no danger of 
confusion. ) 


Suppose now that x and y are dependent on two variables, wu and v say 
(rather than a single variable t). Then the previous argument applies with 
one minor difference. Keeping v fixed, the derivatives in the chain rule now 
become partial derivatives, and we have 


Ow Owodx Owoy 
Ou Ox Ou" Oy du’ 
Similarly, if we keep wu fixed, then we have 
Ow Owodx Owdody 
dv Ox dv Oy Av’ 


These are the relationships that enable us to change from one coordinate 
system to another. In Section 2.4, we demonstrate how to switch between 
derivatives involving Cartesian and polar coordinates. 


Exercise 2.8 


dw MY 
(a) Given w = e”¥, where x = sint and y = cost, calculate 4 first by substituting 


for x and y, and then by using the chain rule. 


(b) Given a function w = w(x,y), we make the substitutions = au + bu and 
y = cu+ dv (where a, b, c and d are constants) to express w as a function of 
u and v. Express the partial derivatives 


ow ww wy Ow 
Ou’ Ov’ Ou?’ Oudv Ov? 


in terms of partial derivatives with respect to x and y. 


2.2.9 Multiple integrals 


In Chapter 1, we introduced the idea of definite integration in a single vari- 
able as the limit of a sum. In this section we generalise this idea to integra- 
tion over two or more variables. We begin with the case of integration over 
two variables. 


Integrals over plane surfaces 


We consider a function f(r) = f(x,y) defined on a region S of the Carte- 
sian plane. We first divide the region into a large number, N, of small 
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A pair of functions 
(x(t), y(t)) defines a 


parametric curve. 


Note that we use partial 
derivatives (O) when 
differentiating w with 
respect to x and y, but 
total derivatives (d) when 
differentiating x and y 
with respect to t. This is 
because w(x, y) is a 
function of two variables, 
but x(t) and y(t) are 
functions of a single 
variable. 


These small elements are 


not necessarily rectangular, 
although it is convenient to 
define them by a grid as in 
Figure 2.7. 


elements 6Am, as shown in Figure 2.7(a), so that the sum of these elements 
approximates the region S. 
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(b) 
Figure 2.7 Regions divided into rectangular subregions 


For each such region we choose a vector r corresponding to a point lying 
in )A,,. Then we define the integral of f over S to be the limit of the sum 


N 
» fay ones 
m=1 


as N tends to infinity (and the size of the elements shrinks), so that 


N 
[ f(r)dA= yim) 2 f (Tm) 5Am- 


It is not essential to form the small regions using a Cartesian grid as in 
Figure 2.7, although this may enable us to write the integral in a more 
recognisable form. 


The simplest case is when the region S is a rectangle, as shown in Fig- 
ure 2.7(b). In this case we have dAm = 0x; Oyj, and the sum over m can 
be found by first keeping j fixed and summing along a horizontal strip (i.e. 
over i), then summing over j. This process corresponds to first integrating 
f(x,y) with respect to x while keeping y fixed, and then integrating the 
result with respect to y. Thus the integral becomes 


[ima f° ([ He.) ar) dy 


(summing first across an arbitrary horizontal strip, and then vertically). 
Alternatively, we could reverse the order of summation to obtain the same 
value: 


[imaars f (f° Fe.)ay) dar 


(summing first up an arbitrary vertical strip, and then horizontally). These 
forms are known as repeated integrals. 


In our first example we shall consider a particularly simple case, namely 
where f(x,y) =1. 


Example 2.1 | 
1 
Evaluate the repeated integral / ( / 
0 x 


Solution 
We have 
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Changing the order of 
summation can produce a 
different value, but only in 
very exceptional cases, 
which will not be 
encountered in this course. 
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Note that in later chapters we use the notation / dx f(x) rather than 


[fo dx, because writing [ae [a f(x,y) in place of [( [0 y) iy) dx 


avoids the need for brackets. 


Example 2.12 
Evaluate the integral | f(r)dA where f(r) = f(x,y) = x7y and S is the 


rectangle with its sides parallel to the axes for which one vertex is at the 
origin and the diagonally opposite vertex is at the point (1,2), as shown in 
Figure 2.8. 


Solution 
We have 
1 2 1 2 
[amas | (| yy ix = | x (/ yay dx 
S 0 0 0 0 
‘ 2 
=| x [gy"] de 


A slight adaptation of the method allows us to integrate over more compli- 

cated regions, as in the following example. 

Example 2.13 

Evaluate the integral / f(r) dA where f(r) = f(z,y) = x7y and S is the 
S 


region bounded by the curves y = x? and y = 1— 2”, for x > 0, and the 
y-axis (so S is defined by the inequalities y > «7, y < 1— x? and x > 0; see 
Figure 2.9). 


Solution 


Figure 2.9 An elemental ‘rectangular’ area in the region S 


Writing the integral as a repeated integral, where we integrate first along a 
strip from one curve to the other, as illustrated in Figure 2.9, we have 


[ie dA = [ve (vay) ae 


1//2 
/ x*(1—2x*) dx = /2/60. & 
) 


| 
NS) Te 
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Figure 2.8 


It is often useful to draw a 
diagram of the integration 
region, as this helps in the 
construction of the 
integration limits. 
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You may wish to repeat this calculation with the order of integration re- 
versed, i.e. perform the z integral first, and then the y integral. However, be 
warned that it will be necessary to alter the boundary conditions. It may 
be instructive to first examine the solution to Exercise 2.9, which solves a 
similar double integral in two ways. 


Exercise 2.9 


Evaluate the integral / x? dA where S is the region in the first quadrant bounded 


Ss 
by the hyperbola zy = 16 and the lines y = zx, y= Oandz = 8 (soy<27,0< 2 <8, 
y > 0 and y < 16/z). 


Integrals involving polar coordinates 


Recall that a point (x,y) in the plane can also be labelled by polar coordi- See page 24 of Chapter 1. 
nates (r,0), where x = rcos@ and y = rsin@. 


The preceding method can be extended to integrals over areas defined in 
terms of polar coordinates. 


Example 2.14 


Figure 2.10 The cardioid r = 2(1 — cos@) 


The curve shown in Figure 2.10 is known as a cardiozd, and is defined in terms 
of polar coordinates by the equation r = 2(1 — cos@) where —7 <0 <7. 


Find the area of the region bounded by the curve. 


Solution 


Figure 2.11 The region bounded by the cardioid divided into elements that are 
approximately rectangular 


In order to find the area S bounded by the curve, we divide the region into 
small elements as shown in Figure 2.11(a). The element shown is approx- 
imately rectangular, so its area is approximately r 60 dr. Summing over all 
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such elements produces an integral of the form / / r dr dé, and writing this 
as a repeated integral we have + 


T 2(1—cos 6) 
[[rarae = f / rdr | d@ 
S —T 0 


= | 4(1—cos6)*d0=67. & 
0 


The preceding example illustrates a general rule which allows us to transform 
between integrals of functions using Cartesian and polar coordinates. For 
any function f and integration region S, 


[[ feuavay = |] rf(rc0s0,rsin®) a0 ar, 


where the region S and the function f are expressed on the left-hand side 
in Cartesian coordinates, and on the right-hand side in polar coordinates. 
Transforming between coordinates is often essential for performing an inte- 
gral, as the following example illustrates. 


Example 2.15 


Show that for a > 0, the integral of the Gaussian function satisfies 


Hey / i exp(—ax?) dx = = 


Solution 


We first make the change of variable u = ,/az, so 
i 1 

Ifa) = — | exp(—a* de = — = F(T). 

Va J—o0 va 


So it is necessary to evaluate J(1). We do this by considering its square: 


wo? (fier) (fia) 
f(fcrre)n 


We transform this integral to polar coordinates (r,@), remembering that 


r? = 7* + y? and that the region of integration is the whole two-dimensional 


plane 0 <=. f= 22,0 < 7 = ©: 


(I(1))? = [ ( [er ao) dr 


Hence we have shown that J(1) = /z, so I(a) = i as required. 
V a 


| 
bo 
= 
=, 
6 
o 
7 
3 
ol 
3 


|| 

| 

~ 
ie 
5, 
% 

| 

~ 


Exercise 2.10 

Evaluate the integral / r? dr d@, where S is the region defined by 0 < r < 1 and 
Ss 

O< 0 < e732. 
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This result was quoted in 
Chapter 1 as 
equation (1.4). 
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Triple integrals and multiple integrals in general 


The notion of a double integral can be generalised easily to integrals of 

functions of three variables. We consider a function f(r) = f(x,y,z) where We assume that you have 
r = xi+yj + zk is defined on a region R of three-dimensional space. We met vectors before; we 
then divide the region R into a large number, N say, of small regions, each shall refresh your memory 
of volume AV;, and the triple integral is defined as a limit of a sum: i eee 


ITI. Ly eal » Pinte, HOV. 


Such integrals are often known as volume integrals. 


In this course multiple integrals generally arise from theoretical considera- 
tions, and we rarely need to evaluate them except in cases where the sym- 
metry, or other factors, makes the calculation particularly easy. However, 
as an aid to understanding, it may help to see how we might undertake such 
a calculation. For this purpose we shall choose a very simple region. 


Example 2.16 


> )——__} 
. im 


Figure 2.12 The rectangular box R divided into elemental subvolumes 


In this example we suppose that R is the rectangular box defined by the 
inequalities 


(ee er ee oo ee ee ae 
(as shown in Figure 2.12), and that 
f2y2)=e > ee 
We then have 


= [ (fe eaparep Ce) iy) its 


we s(1 2 oe +o"). 


Integrals over surfaces that are not planar 


We define integrals over surfaces that are not planar in a similar fashion, 
first dividing the surface into a large number of small elements, then forming 
a sum, and finally taking the limit. In practice we try to choose a coordinate 
system which allows us to write the integral as a repeated integral, as in the 
following example. 
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Example 2.17 
Given f(r) = 3z — \/x? + y?, evaluate the integral / f(r) dA where S is 
S 


the curved surface of a cylinder of height h/2 and radius R with its axis on 
the z-axis and its base on the plane z = h/2. 


Solution 


Figure 2.13 ‘The curved surface of a cylinder divided into elemental rectangles 


In this case it is helpful to choose cylindrical polar coordinates (r, 6, z) where 
x=rcos@ and y=rsin@. Then f(r) = 3z— /22+ y? = 3z—7r, and an 
area element on the surface is 0A = R06 0z (see Figure 2.13). 


On the curved surface of the cylinder we have r = R, so 


Lie nae 2 ([ _R)R io) se 
a i ([- id) 62 ee 


h 

= an | (32R — R°) dz 
h/2 

= thr(9h-4R)R. OW 


Most of the integrals of this kind that will arise in this course are particularly 
simple. For example, if f(r) = 1 and S is a sphere of radius p with its centre 
at the origin, then 


[im iA= [aa = Arp’. 


Similarly, if f(r) = |r|? + 2, then 


. [ teyaa=(e +2) f dd =anp?(o? +2), 
S S 


End-of-section Exercises 


Exercise 2.1] 


O7w Ow 
Find the partial derivatives —~ and —— if w = 


Ox? Oy? Vx? + y2 
Exercise 2.12 


Verify the mixed derivative theorem (see page 66) for the function 
- 


Se) = 


Recall that a sphere of 
radius p has surface area 
Anp*. 
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Exercise 2.13 
Show that each of the following functions is a solution of 
07 w r < 
Ox? Oy? 


the two-dimensional form of Laplace’s equation in Cartesian coordinates. 


0, 


(a) w=e* cosy (b) w= 2° — 3zy? 

Exercise 2.14 

Find and classify the stationary points of the function 
F (x,y) =1— 22x — 12y + 10x? + 2xy + 5y”. 

Exercise 2.15 


Find the point on the plane 2x — y + 2z = 16 that is nearest the origin. 


Exercise 2.16 
Given w = 1/2? + y? + 22, calculate 
Ow Ow aod OPw 
OxOyOz’ OzOxOY OyOzOx 


Exercise 2.17 
7 cos 0 
Evaluate the repeated integral / ( / psn dp dé. 
0 0 


Exercise 2.18 


The region S in the (x, y)-plane is bounded by the three lines x = 3, y = 2/3 and 
y = 0. Write the integral of the function exp(x?) over the region S as a repeated 
integral in two ways, namely with the x integral performed first, and then with the 


y integral performed first. Use one of these to evaluate / exp(x”) dA. 
S 


2.3 Vector calculus 


In this section we shall be concerned mainly with the behaviour of func- 
tions whose domains are sets of vectors, so we now revise some elementary 
properties of vectors. 


2.3.1 Vector algebra 


We assume that you are familiar with the fundamental notions of vector 
algebra. This is essentially the material that you are likely to have met in 
previous Open University courses, but for students who have gained their 
knowledge of the subject from a different source, we provide a very brief 
summary and explain our notation. 


A brief summary 


Scalars are quantities that are defined by their magnitude, such as mass 
and length, and we use the term to distinguish them from vector quantities, 
which have both magnitude and direction, such as velocity. Vectors are often 
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This equation will be 
discussed briefly in 
Section 2.4. 


This exercise 1s optional. 
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denoted in printed material by bold type, for example v or F’, or by under- 
lining in handwritten material, v or F (although several other conventions 
are in common use). They are often specified by pairs of numbers in two 
dimensions, and triples of numbers in three dimensions: for example, (2, y) 
in two-dimensional Cartesian coordinates, or (x, y, z) in three dimensions. 


The following statements are made in the context of three-dimensional vec- 
tors, but similar results apply in two dimensions. The vector (x,y) in two 
dimensions is often equated with the vector (x,y,0) in three dimensions. 


A vector r = (x,y,z) may be multiplied by a scalar a to give 

ar = (ax, ay, az), 
and two vectors in Cartesian form may be added to give 

ri tre = (£1, 91,21) + (2, yo, 22) = (v1 + Za, yr + Y2, 21 + 22). 
The zero vector O corresponds to (0,0,0) in Cartesian coordinates. 


The magnitude or length of a vector r = (x,y,z) is denoted by |r| and 
defined to be |r| = \/a?2 + y?+ 27. A unit vector is a vector of magni- 
tude 1. The unit vectors in the directions of the Cartesian coordinate axes 
are labelled i, j and k, so a three-dimensional Cartesian vector r = (2, y, 2) 
may also be written in the form r = zi+ yj+ zk. The scalars z, y and z 
are known as the Cartesian components of the vector r. 


We write a ‘hat’ symbol over a vector to denote the fact that it is a unit 
vector. Thus a unit vector in the direction of a given vector r is denoted by 


in| 


7, where r = 


The scalar product 


The scalar (or dot) product of two vectors a and 6 is 
a-b= |a||b| cos 8, (2.4) 


where @ is the angle between the vectors. Using a unit vector b, the scalar 
product a- b represents the length L shown in Figure 2.14; this is known as 
the projection of the vector a in the direction of b. Note that the scalar 
product of two vectors is a scalar. 


Figure 2.14 The projection of a vector a in the direction of b 
In Cartesian components, the scalar product of vectors a = aji+ agj + a3k 
and b = 611+ boj + b3k can be written as 

a-b=aj,b, + agbe + agzbs. CRS 


Note that it is not immediately obvious that equations (2.4) and (2.5) are 
equivalent. However, you should be prepared to use either definition as the 
need arises. 


The following list includes some useful properties of the scalar product. 


8 | 


Strictly speaking, (2, y, z) 
is a representation of the 
vector r rather than the 
vector itself, but we shall 
not be concerned with such 
fine distinctions. 
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Some properties of vectors 


t=) je bbe’, i jie, 
a-b= |a||b|cosé 

= (a1i + aaj + agk) - (bi + b2j + 63k) 

= ab, + agbe + a3b3 

=, 
la| = \/af + a5 + a5 = J/a-a, 
|Aa| = |A||a| for any scalar 4, 
(a+b)-c=a-c+b-e, 
(A + w)a = Aa+ ya for any scalars 4 and w. 


Example 2.18 


For a fixed vector v, and any unit vector n, show that v-n has maximum 
value |v|, when 7 is chosen to be parallel to v. 


Solution 
Using equation (2.4), 
v-n = |v| |n| cos @ = |v| cos 0. 


For any given v, this has a maximum value |v| when cos@ = 1. This corre- 
sponds to 6 = 0, i.e. to choosing 7 to be parallel tov. I 


Exercise 2.19 


(a) Find the cosine of the angle @ between the vectors a = (1,2,—1) and 
b= (2,1, 1). 


(b) Find the magnitude of the vector c = 3i + 2j — 5k. 


2.3.2 Scalar fields 


In applied mathematics, real functions of two and three real variables are 
often known as scalar fields. The term ‘scalar’ is used to distinguish them 
from the ‘vector fields’ that we shall discuss shortly. We begin with a physical 
situation which may help you to understand the general cases that come 
later. 


If a drop of ink falls onto a large sheet of wet blotting paper, then the stain 
will spread. The density F' of the stain (i.e. the mass of ink per unit area) 
at a particular instant of time might be represented quite well by a function 


Pie; y) a Aexp [—a(a* - y)] : 


where A and a are positive constants (and the blotting paper is assumed to 
lie in the (a, y)-plane with the centre of the stain at the origin). In this case 
F (x,y) is an example of a two-dimensional scalar field. 


Suppose that time is kept fixed, and imagine that we follow a curve drawn 
on the paper through the ink stain: we shall notice a change in the ink 
density as we traverse the path. At any particular point on the path, we 
intend to determine the direction in which the ink density increases most 
rapidly. 
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If we represent F(x, y) as a surface, as in Figure 2.15(b), then the height of 
the surface above the point (x,y) represents the density of the stain at that 
point. 


F(x, y) 


(b) . 


(a) 


Figure 2.15 (a) The ink stain in the (x, y)-plane; (b) the surface representing 
F(2,y) = Aexp[—a(2? + y”)| 


The gradient 


At most points on the surface representing F’ we should be able to deter- 
mine both the direction in which the slope is greatest and the magnitude of 
this greatest slope. We can then associate a vector having this magnitude 
and direction with each point (2, y). Since the magnitude and direction of 
greatest slope will vary from point to point, the magnitude and direction 
of the vector representing this will also vary from point to point. In other 
words, this vector will be a function of position (x,y). This vector func- 
tion is known as the gradient of F’, denoted by VF’; we sometimes write Read V as ‘grad’. 


VF (x,y) to emphasise its positional dependence. Some authors write grad F’ 
instead of VF’. 


For our ink stain illustration, WF always points in the direction for which 

the ink density increases most rapidly. In Figure 2.15 we show VF ' at two 

points, (x,y) and (2z’, y’), by the darker arrows. If we were to draw the level The lighter arrows point 
curves for the surface representing F’, then we would obtain a diagram very uP the lines of greatest 
similar to Figure 2.2 (page 59), and we would see that the gradient at each slope, and the darker 


point is perpendicular to the tangent to the level curve that passes through ee aes their Pao 
that point. onto the (x, y)-plane. 


We now formally define the gradient of a function, and state some properties 
associated with it. 


The gradient 


The gradient of any function F(x, y) is defined as 


OF. OF 
fF = —i+ —-j. 
= ies ay? 


At each point (z,y), WF’ is a vector which points in the direction of 
maximum slope (uphill) of F’; its magnitude |W F| is the value of the 
maximum slope at that point. 


Additionally, if nm is any unit vector, then VF'- 7 gives the slope of F’ 
in the direction of n. 


The proof of these properties is presented as an optional example. 
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Example 2.19 This example is optional. 


Prove the aforementioned properties of VF’. 


Solution 


We first need to find the slope of the surface representing F’ in an arbitrary 
direction, say along the direction of the unit vector n = cos#i+ sin 6j. 


We consider the change in the value of F' between two adjacent points on the 
surface, P and Q say, with coordinates (x,y,z) and (x + 6x, y + dy, z + 62), 
respectively, as shown in Figure 2.16. 


Figure 2.16 The points P and Q on the surface representing F' 


We now keep P fixed and consider the slope of the line PQ. If we let 
6s = ,/dx2 + dy?, then we see from Figure 2.16 that the slope of the line 


se We know (from equation (2.1)) that 


joining P and Q@ is given by ‘ 
S 


deca ba yoy, 


where the partial derivatives are evaluated at P. Therefore 


() () ‘) 
slope of PQ = “ ~ Fp _ iy =(f,c0s6 +2410. 


This result holds in the limit as 6s — 0, so 
slope in the direction of n = (F,i+ F,j)-n= VF +n. 


So, keeping P fixed so that F;, and F, are constant, we can calculate the slope Of course, this argument 
in any direction simply by varying the direction of N. From Example 2.18, fails if F2 = Fy = 0, since 
we see that the maximum slope has value |W F'|, and this occurs when VF’ then there is no direction 
om of maximum slope. 

is parallel ton. & 


We mentioned earlier that at every point VF' is perpendicular to the tan- 
gents of the level curves of F’, but we shall not prove this here. 


In the case of the ink stain example, we have F(x, y) = Aexp |—a(a* + y*)| ,; 


SO 

OF 2 2 

— —2a Ax exp|—a(x +y )] 
and 

OF 

Du = —2aAy exp [—a(a* + Tm) y 
therefore 


VF =-—2aA exp|—a(x* ~ y)] (xi + yj). 
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This is a vector in the same direction as a vector from the point (x,y) to 
the origin, which is what we would expect from the symmetry shown in 
Figure 2.15. 


Similar ideas can also be applied to a three-dimensional scalar field F(x, y, z), 
for which we define 


oF. OF, OF 
VF = —i+ —j+—k 
az Oy O 
O O 
It is sometimes useful to regard V = i— + j— +k-— as an operator that 
Ox Oy Oz 


can be applied to functions, in which case it is known as a vector differ- 
ential operator. 


Example 2.20 
Calculate VF for the function F(x, y,z) = x*ysin z. 


Solution 
We have 
OF OF OF 
aoe 2ry sin Z, on =r’sinz and eS ry COS Z, 
* OF, OF, OF 
r= 44 —j+- 6 
Ox = dy? i Oz 


2 


= (2cysin z)i + (x* sin z)j+(2°ycosz)k. I 


Example 2.21 


Calculate the gradient of the function F(a, y) = x? — xy + 3y at the point 
(2, 1). 


Solution 


The point (2, —1) refers to a point in the domain of F (i.e. the point (2, —1, 9) 
on the corresponding surface). F,(z,y) = 2x — y and F,(z,y) = —a + by, 
so F,(2,-1) =5 and F,(2,—1) = —8. It follows that VF(2, —1) = 5i — 8). 


A schematic representation of VF’, together with the level curves of F’, is 
shown in Figure 2.17. The arrow at each point represents the vector VF’ at 
that point. Note that at each point, VF' is perpendicular to the tangent to 
the level curve of F’. 
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Figure 2.17 The level curves and gradient of F(z, y) = ac? —agyt+3y? 


V is known as del or 
nabla in some texts. 
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Exercise 2.20 


Calculate the gradient of the function F(zx,y) = (x? + 2y7)sin[(a + y)z] at the 
point (1,1). 


2.3.3 Vector fields 


A vector-valued function F(x, y) = Fi(a,y)i + Fo(x, y)j in two dimensions, 
or F(2,y,z) = Fy(2,y, z)i+ Fo(2, y, z)j + F3(x, y, z)k in three dimensions, 
is often known as a vector field, and we commonly make abbreviations like 
F = F\i+ Foj when the choice of independent variables is clear. 


We saw in the previous subsection that the gradient of a function is a vec- 
tor field; but not all vector fields arise as gradients. Physical examples of 
vector fields include fluid velocity (where the vector at each point is equal 
in magnitude and direction to the fluid velocity), and the force on a charged 
particle in an electric field (where the vector at each point is equal to the 
force on the particle at that point). 


The spatial variation of a vector field is inevitably more complicated than 
that of a scalar field: for example, in three dimensions, it involves the nine 
partial derivatives 

OF, OF, OF, OF» OF) OF» OF OF OF 

Be: Oy” Of Ce Oy’ Ge ae Be Oz 
These partial derivatives completely describe the spatial variation of F’, 
and certain combinations of them occur quite often. The combination that 


interests us here corresponds to the scalar field known as the divergence of 
F,, defined as the scalar product of the operator V with the vector field F: 


Ga Se : s Some authors use the 
ie iy + ay? = aoe (Pelee Fgj4 Bek) notation div F' in place of 
V -F. 
_ OF, , OF , oF 
Ox Oy ig 


Note that although F is a vector field, V - Fis ascalar field. The importance 
of V - F will become clearer when we discuss the Gauss divergence theorem; 
however, this theorem involves integration over surfaces and volumes, so first 
we need to say a little more about integration. 


Exercise 2.21 

1 . 
(x2 + y2 + 22)1/2 
(b) Given f(x,y, z) = 277 + y? 4+ 27, calculate V+ (Vf). 


(a) Given F(z, y,z) = (xi + yj + zk), calculate V- F. 


2.3.4 Integration 


We now return to the topic of integration. We begin with the problem of 
how to determine the length of a curve, before moving on to discuss the 
integration of vector-valued functions. 
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The length of a curve 


We continue with the idea of integration being the limit of a sum. 


~ QW, 
a ~, 
& x & 
s 


rf 
ie?) 
, os 


Figure 2.18 A curve divided into short subintervals 


Suppose that we wish to determine the length L of the curve y = f(z) 
between x = a and x = b, as shown in Figure 2.18. One approach might be 
to divide the curve into a large number (say N) of small segments that are 
almost linear. The sum of the approximate lengths of these small segments 
would then give us an approximation to the total length of the curve. If 
we increase the number of such subdivisions, we might expect to obtain an 
increasingly accurate approximation to the total length. In other words, the 
length of the curve could be represented as the limit of a sum, i.e. as an 
integral. 


For the length 6s; of the kth segment, we have bs%, ~ bx? “bs dys. So 


N N N 2 
L~ "bsp = > 1/502 + dy2 => ban 1+ (3) 
k=! Ral k=] 


and in the limit as N — c we have 


b 2 
Lz / =. (3) az. 
\ dx 


Example 2.22 


Calculate the length of the graph of the function y = x3/? from the origin 
to the point (1,1). 


Solution 


The required length is 


1 dy 2 1 5 
i= | 1+ (2) dx = | 1+ ($21/2) dx 
0 dx 0 
1 
=3/ V4 + 9x dx 
0 


= (4 + 92)9/?] 
= 7-(13V13—8). I 


The same idea can be extended to a curve that is not the graph of a function 
(i.e. not of the form y = f(x)), for example the curve shown in Figure 2.19. 
Such curves often arise in the parametrised form x = z(t), y = y(t) with 
ty <t < te, and are called parametric curves. A smooth curve is defined to 
be a curve for which the functions x(t) and y(t) have continuous derivatives. 
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° ai ° = ‘ ° 37r 37 
Figure 2.19 The curve x(t) = sin 2t, y(t) =sint with —3' <t< 4 
The curve in Figure 2.19 can be specified in parametrised form by putting 
x(t) = sin 2t and y(t) = sint with —27 <t < %. 


In general, if s denotes the distance along the curve from any fixed point, 
we have 6s? ~ dx? + dy? as before, therefore 


6s\? ba \* oy : 
aes ~ ee -+- as . 
ot ot ot 


and in the limit as 6¢t — 0 we have 
ds\" _ (dx\" | (dy\" 
aij eg dt) ~ 


d dx\?  (dy\? 
Thus = = (F) + (F) , so for the length L of the curve shown in 


Figure 2.19 we have 


37/4 
—37/4 
37/4 
= V 4cos? 2t + cos? t dt ~ 6.917 84. The integral was evaluated 
—3r/4 by numerical methods. 


Such line integrals arise in a number of areas of applied mathematics. For 
example, it is relatively easy to calculate the work done when an object 
is pushed up a straight path against a constant force F’, as shown in Fig- 
ure 2.20. We simply multiply the length of the path d by the component F’ 
of the force along the line of the path. 


Figure 2.20 A constant force acting along a straight path 


Imagine now that we push an object up a winding cliff path, and that we 
wish to find the total work done. In theory we could divide the path into a 
large number of almost linear sections, then calculate the work done when 
the object is pushed along the many short straight paths, and thus obtain 
an approximation to the total work done. We expect the approximation to 
improve as the number of subdivisions increases. 


In the following subsection we shall make this idea a little more precise. 
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Line integrals of vector fields 


Suppose that we are given a parametric curve C' defined by a vector r(t), 
so the points on the path are determined by the values of the parameter tf. 
Now suppose that a particle is propelled along the curve by a force F(r), 
and that we wish to calculate the work done by the force. If a constant force 
F acts on a particle and displaces it through a vector distance ds, then the 
work done by the force is F'- ds. So in order to calculate the total work 
done by the variable force F(r) as it moves along a curve, we approximate 
the path by a large number of straight-line subpaths on which the force and 
the direction are approximately constant (see Figure 2.21). We then have 


N 
total work done ~ 2 F(rpz) + org. 
k=1 


r(t) 


Figure 2.2! A subpath dr, on the curve r(t) 


This approximation will improve as N increases (and the maximum length 
of the subpaths decreases), leading us to define a line integral as 
N 
J =_ jm Pig * es 
N—-oo 


n=! 


which we shall write in a more recognisable form shortly. 


Such integrals can also be used to represent quantities other than the work 
done. 


Suppose now that C is a curve in three dimensions, and that it can be 
parametrised by writing 

r(t) = 2(t)it+ y(t)j + z(¢)k, 
with the beginning and end points of C corresponding to t=a and t= 6, 
respectively. On C we have 6r = 6x1 + dyj + 6zk, so x = ai _ us =the 


where dt is the increment in t that gives rise to the increment dr, therefore 
d d d 
or ~ (Fi ee: a) t = (a'(t)it y/ (Hj + 2’ (Hk) ot. 


If F(r) = Fy(a2,y,z)it+ Fo(a, y, z)j + F(x, y, z)k, so that on C the compo- 
nents of F are functions of t, then 


N 
2 FP, )* OF » 


where t; = a and ty = b, so Not = b—a. 


89 


In two dimensions, 


r(t) = (#(t), y(t)); 


in three dimensions, 


r(t) = (a(t), y(t), z(t). 


The i component of F is 
F, (a(t), y(t), 2(t)) = Fi), 
and similarly for the other 
components. 


90 Chapter 2 Multivariable and vector calculus 


In the limit as N — ov, 


b b 
dx d dz dr 
I= / (AO + x(t) 4 A) | dt = if B. (F) dt 


(which is just a ‘normal’ integral with respect to f). 


The line integral J is usually written as 


[Fo -dr, 


where it is understood that the curve C’ has been parametrised so that 
r = x(t)i+ y(t)j + z(t)k is a function of a parameter t. 


The following example (in two dimensions) may help to clarify the notion 
of a line integral. 


Example 2.23 


Evaluate the line integral / F(r)-dr where F(r) = x7i+ 2’y?j and C is 


C 
the part of a parabola defined by the equation y = 1 — x? from A, the point 
(0,1), to B, the point (1,0), as shown in Figure 2.22. 


Figure 2.22 Part of the parabola y = 1 — x” 


Solution 
First we parametrise the curve by writing 2(t) = t and y(t) = 1—??, with 
0O<t<1. Then ¢t =0 and t = 1 correspond to the points A and B, respec- 


o = —2t. We then have 
dt 


Fir) = eit vy? j — i+ ra = Py 


d 
tively. Also notice that = ='1 atid 


and 
dr wd —— : 
gt Se (a(t)i oe y(t)j) = i— 2tj, 
SO | 
1 
i F(r)-+dr =| (t7i + t?(1 — t7)?j) - (i— 2t§) dt 
C 0 


1 
~ (t?—2¢°(1-¢t*)’)dt=}. I 
) 


The evaluation of such integrals in three dimensions is very similar. 


Exercise 2.22 


(a) Calculate the length of the curve x = t?, y=?t? from t=0 tot =4. 
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(b) At time ¢ seconds, the position of a particle is given by x =t, y= t?//2, 
z = t?/3. How far does the particle travel in the two seconds from t = 0 to 
j= Ze 


(c) Evaluate / F -dr when F(z, y,z) = v*yi+ zj + ry’*k and C is the straight 
‘s 
line joining the origin to the point (1, 2, 3). 


2.3.5 The Gauss divergence theorem 


We can now move on to discuss the Gauss divergence theorem, which ‘in- 
volves an integral over a volume and an integral over the corresponding 
surface. 


The flux of a vector field across a surface 


The concept of the flux of a vector field across a surface arises in many 
areas of applied mathematics, but perhaps the easiest to visualise is the 
flow of a fluid or gas across a surface. Imagine the flow of a fluid in three 
dimensions, where at each point P with position vector r = (x,y, z) the fluid 
has velocity F'(r); i.e. the fluid velocity is described by the vector field F'(r). 
Now imagine that we immerse a wire mesh in the fluid, and that we wish to 
estimate the volume of fluid passing through the mesh in one second. The 
situation is depicted in Figure 2.23(a), where the mesh is shaped rather like 
the visor of a motorcycle helmet. 


Figure 2.23 A surface approximated by a ‘wire mesh’, and the volume of fluid 
passing through an aperture in one second 


First we consider just one small aperture in the mesh, and calculate the 
volume of fluid passing through it. Suppose that the aperture is at position P 
with area AA, and that fluid passing through this aperture has velocity F’. 
Also suppose that the vector 7 is a unit normal to the surface at P. 


A vector perpendicular to the tangent plane at a point P on a surface is known 
as the normal to the surface at P. At every point on a smooth surface 
there are two normals, one ‘into the surface’ and one ‘out of the surface’. If 
the former is denoted by n, then the latter will be —n. For a closed surface 
(i.e. a surface without a boundary) we usually take the outward normal to be 
positive. 


In general, the fluid velocity and the unit normal will not be parallel, and 
in Figure 2.23(b) the shaded volume represents the fluid passing through 
the area AA in one second. This is a parallelepiped with base area AA and 
length F (see Figure 2.24). 


9 | 


This is not a real wire 
mesh, but rather a 
mathematical construct 
with infinitesimally thin 
wires, and ultimately with 
an infinitely fine mesh. 


If the aperture is small 
enough, then the fluid 
velocity will not vary 
across it. 
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Figure 2.24 


From the diagram we see that the volume of this parallelepiped is equal to 
the volume of a box with base area AA and height F'-n, where n is the 
unit normal at P. Thus it is the component of velocity perpendicular to the 
mesh that is relevant, and the volume of fluid passing through the aperture 
in one second is approximately F(r)-n AA. 


If we sum over all apertures, and form the limit of the sum as the size of the 
mesh shrinks to zero and the number of subdivisions increases, we obtain 
an integral 


/ Cnr 

S 

This integral is the flux of the vector field F' across the surface S, which in 

this case represents the volume of liquid crossing the surface in one second. 

It is often convenient to abbreviate 2 dA to dA, in which case the integral 

is written | F(r)-dA. 

S 

Example 2.24 

If the velocity of a fluid is given by the vector field 
F = cay?2*i+ x? y2z7j+ ry’ zk, 


calculate the volume of fluid per second passing through the rectangular 
area S given by r=1,0<y<1,0<2<1. This situation is depicted in 
Figure 2.25. 


Figure 2.25 


Solution 


The unit normal to S is given by n = i, sodA =idA =idydz. Hence the 
volume of fluid passing through S is 


v= FQ.y2)-dA 
S 
1 


1 
=| (| (y227i + y27j + y*zk) - idy dz 
0 0 
1 1 1 
=|} (| yPetdy) ae = | a2 dz = 5. 
0 0 0 


We say that the flux of fluid across the surface S is 5. If the fluid velocity is 
in metres per second, and the coordinates along the z- and y-axes represent 
metres, then S has area 1 m?, and + m° s~! of fluid passes through S. 


The vector dA = ndA has 
magnitude dA and 
direction normal to the 
surface. Since surfaces are 
in general curved, the 
direction of n and dA will 
depend on r. 


Note that we could have 
chosen n = —i, in which 
case we would be 
calculating the volume of 
fluid passing through S in 
the —i direction, which 
would simply give a change 
of sign in the final result. 
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More generally, consider the flux of an arbitrary vector field F’ across the 
surface S of a small rectangular box with sides of length dx, dy and dz, as 
shown in Figure 2.26. 


F(x, y, Z) F(x + x, y, Z) 


Figure 2.26 An elemental rectangular box 


The flux across the right-hand shaded area of this box is approximately 
F(x + 6x,y,z)+idy dz, whereas the flux across the left-hand shaded area is 
approximately —F(z, y,z) -idy dz (because here the outward normal is —i). 


Applying a similar argument to the other two pairs of parallel surfaces, we 
obtain an approximation for the flux across the entire surface S: 


[ Fedax (F(x + 62, y,z) — F(x, y, z)) -idy dz 
S 
a 


(F(x, y + dy, z) — F(z, y,z)) + dx dz 
+ (F(a,y,2+ 62) — F(z, y, 2)) + kdax by 


=7y (Seen — Fi(a, y, 2) 


Ox 
Fo(x,y + dy, z) — Fo(a, y, 2) 
+ sates isi iia ait tae te necrgneS 
oy 
F3(x,y, 2+ 6z) — F3(a, y, =) 
fe eee |, 
Oz 
where 6V = 6x dy 6z. If we now take the limit as the volume shrinks to zero, 
we obtain 
1 OF, OF) OF3 
lim — | F-dA=—+4+—4——=V°-F 2.6 
Jim, sp | a Oy wae ea 


where V is the vector differential operator discussed on page 85. This result 
has an interesting consequence when F' describes the velocity of a fluid. 


The integral / F -dA represents the volume of fluid crossing the surface in 
S 


one second (both into and out of the region). For an incompressible fluid, 
where the amount of fluid entering the region equals the amount exiting, we 
would expect this volume to be zero (unless there is a source of fluid inside 
the region or a sink allowing some fluid to escape). In other words, for F 
to represent the velocity of an incompressible fluid, we require V - F' = 0 
throughout the region. 


The Gauss divergence theorem 


Equation (2.6) is valid for an infinitesimally small volume element dV. The 
Gauss divergence theorem generalises this result to a finite volume R. 
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Gauss divergence theorem 


For any (sufficiently smooth) vector field F', the volume integral of 
V - F over a region R is equal to the (outward) flux of F across the 
surface S enclosing R. This statement can be expressed as 


| vera = | P-aa. 
R S 


Figure 2.27 Schematic representation of the Gauss divergence theorem 


This result is easily established from equation (2.6). First we use a Cartesian 
grid to divide the region R into a large number of box-shaped volumes dV; 
with corresponding surfaces 05;. Then 


/ V-FdV~)S(V-F)iiVi, 
ft i 


where (V- F); means V - F evaluated at some point within the small vol- 
ume 6V;. Using equation (2.6), we obtain 


¢ 3 0S; 


where / F' - dA; represents the flux of F' out of the box-shaped region 6V;j. 
69; 


Figure 2.28 ‘Il'wo adjacent elemental boxes 


Now, adjacent boxes have a rectangular face in common, and the outward 
normals on this common face have opposite signs (as shown in Figure 2.28). 
The flux of F out of this face from one box is equal in magnitude to the flux 
of F into the adjacent box, but these fluxes have opposite signs, so their 
contributions to the sum cancel. This cancellation occurs for all the faces 
except those at the boundary of the region R, so we have 


| v-rav~ | F-aa 
R hed 


As the number of elements in the subdivision increases, and their size de- 
creases, this approximation becomes more accurate, and in the limit we 
obtain the required result. 


The Gauss divergence theorem is a very important result. It relates an 
integral over a region to an integral over the boundary of that region. In 
this sense it is a generalisation of the well-known result in one dimension 
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od 
/ ah dx = f(b) — f(a). In fact, the theorem can be generalised to any Here, the points a and b 

q ax form the ‘boundary’ of the 
number of dimensions. At a practical level, it enables us to calculate the line interval [a,b]. 


flux of a vector field by calculating a volume integral of a scalar field. In 
general, this is much simpler than calculating integrals of vector fields over 
surfaces. 


Example 2.25 


Verify the Gauss divergence theorem for the vector field 


F = cy i+ yzj + 2ek 
and the region # given by the unit cube O< x <1,0<y<1,0<2< 1. 


Solution 


First we integrate V - F over the region R. Since 
O(xy? O O 
Ox Oy Oz 


we have 


[verw=f[ (f ([ (w+? +2)dr) ay) ae 
“(fueroe)s 


1 
=| (2+2)dz= 3. 


Next we integrate the field F’ over the surface of R. Since R is a cube, 
we need to evaluate six integrals, one for each of the faces at x = 0, x = 1, 
Goat, vet. 2=0, 2 = I. 


For the face at x = 0, the unit normal out of the cube is given by n = —i, 
so the surface area element is dA = —idydz. Hence F(0,y,z)-dA = 0, and 
we obtain 
Sy= F(0,y,z)-dA=0. 
=) 


For the face at x = 1, the unit normal out of the cube is given by n = i, so 
the surface area element is dA = idydz. Hence F(1,y,z)+dA = y’ dy dz, 
and we obtain 


1 1 
S,= | F(l.y.2)-dA= | (/ yPdy) de =}. 
ad | 0 0 


For the face at y = 0, we have n = —j, so dA = —jdxdz. Hence we have 
F(x,0,z)-dA =0, and we obtain 


S3 = / F(z,0,z)-dA=0. 
y=0 


For the face at y= 1, we have dA = jdzdz, so F(z,1,z)-dA = zdzrdz, 
and we obtain 


1 1 
s.= | F(z.1,2)-dA= | (/ zd) dz =}. 
qeoed O 0 


For the face at z = 0, we have dA = —kdzdy, so F(x,y,0)-dA =0, and 
we obtain 


$5 = | F(z,y,0)-dA = 0. 
z=() 


96 Chapter 2 Multivariable and vector calculus 


For the face at z = 1, we have dA = kdzdy, so F(z,y,1)-dA =adzdy, 
and we obtain 


1 1 
So= | F(x,y,1)-da= | (| v de dy = 5. 
z=1 ) 0 


Adding these surface integrals gives 
$1 + So + 93+ S4+55+S6 = §, 


the same result as obtained by finding the volume integral. Hf 


Exercise 2.23 


Verify the Gauss divergence theorem for the vector field F = xsin(zy)i+ xzj over 
the region 0 <.¢ < 1,0< 9 = 1,82 = 1. 


2.4 Changing coordinates 


In Subsection 2.2.8, we showed that if the Cartesian coordinates x and y are 
functions of two independent variables u and v, then any function w(z, y) 
can be differentiated with respect to u and v, using the chain rule. First, 
keeping v fixed, we have 


dw _ dwe , dw dy 
Ou Oxdu Oy du 
Now keeping wu fixed, we have 


Ow  Owdx ; Ow Oy 

Ov Ox dv Oy dv 
In this section we use these results to transform between derivatives in polar 
and Cartesian coordinates. 


Polar to Cartesian coordinates 


Recall that in polar coordinates, points (x,y) in the plane are labelled by 
(r,@) where x = rcos@ and y = rsinJ@. 


Now suppose that we have a function w(z,y) and we wish to express its 
partial derivatives with respect to r and @ in terms of its partial derivatives 
with respect to x and y. 


First we find 


= cose, a O a2 Sep aq = 70088, 
r 


O 
Or 
and then, using the chain rule, 


Ow Owdx Owdy Ow 


Br an Oe ee 


and 


Ow _ Ow Ox Ow Oy 


O O 
nce + rcs 6 — = ¢%— — y—. 


B60. On 60 «bee Oar a Oe 


2.4 Changing coordinates 97 


Similarly, we can calculate the second-order partial derivatives with respect The algebra may be 
to r and @. unpleasant, but the process 
is straightforward. 


Example 2.26 


For an arbitrary function w(2,y), calculate the second derivative with re- 
spect to # in terms of the second derivatives with respect to x and y. 


Solution 


We have just shown that 
Ow oe Ow 


00° Oy 7 x’ 
since w is arbitrary, it is convenient to consider the operator 
O O O 
06 By On. 
This allows us to express —s in terms of partial derivatives with respect 


to x and y in the form 
ew O Ow Ow 
ae 0 46 7) "Oy ven) (eo Yas) 
It remains to simplify the expression on the right: 
Oew O O Ow Ow 
wat = (8 ¥5z) (755 -¥ae) 
O Ow O Ow O Ow O Ow 
#3 (8) (8) =r (8) 
sO"w Ow 0? w Ow 0? w > O*w 
se Oy? = & +55 | of (= tos | = Ox? 


20°w Ow ,0°w Ow Ow 


Example 2.27 
2 
Given z = x* + y’, calculate a ee and as in terms of r and @. 
Solution 
We have 


ke Da ee ee 
or @rcr <dytr | 
= 2r(cos? 6 + sin? 6) = 2r, 
G2 O20% .O2ny 
ied abe! 9 
70 = 52 00 | By OO 2x(—rsin 0) + 2y(r cos 6) 
= 2r7(— cos @sin 0 + sin 0 cos @) = 0, 
0*z 
Oroe 


In this example we could have noticed that z = 2? + y? = r?, from which 
the results follow immediately. 1 


= 0. 


Cartesian to polar coordinates 


Again we have the relationships x = rcos@ and y = rsin@, but while pre- 
viously we considered x and y to be functions of r and @, now we suppose 
that r and @ are functions of x and y. 
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2 


Suppose that we wish to express in terms of partial derivatives with 


2 
respect to r and #. One might suspect that it would be useful to write 
r=vV/a?+y2 and @=tan '(y/z) for x > 0, 
and proceed as before. However, using this approach the algebra becomes 
unpleasant, and there is a better way. 


We begin with the equations x = rcos@ and y = rsin#@. First we differenti- 
ate them partially with respect to x, to obtain 


Or aes. Or . le, 
a a and O=— sme + reeset = 


and these equations can be solved to give 


a = cos? ane ve = —_— 


Ox Ox r 
Similarly, differentiating the equations x = rcos@ and y = rsin@ partially 
with respect to y, and solving the resulting equations, we obtain 


Or OO cosé 

oo and a a 
These relationships are exactly what we need in order to find expressions 
for = and = in terms of — and aril From the chain rule, we have 

Ox Oy Or 00 

Ow OwOor Owd0é Ow sinddw 

ac oe 
and 

Ow OwOr Ow0é _ Ow cos#dw 

Oy. aca. Olay Wee fee 


Ow . 
If we now replace w by a in the first of these expressions, we obtain 
x 


00 _ ogg 2 (Ow) _ sind d (dw 
a ais Ox r 00\ 0x)’ 


, 1 Ow 
and using the above expression for Dan , this becomes 
x 


cos? — — 


Or r 00 


nn = 


Ox? — Or Or r 00 oly 
Finally, using the rule for differentiating a product, we have 
2 0?w 2sin@cos@ 0?7w — sin? 6 d?w 

Ox? Or? r Oro00 r2 66? 

sin” 0 Ow . 2 sin 6 cos @ Ow 

rar i 0607 
In the following example and exercise we assume that a function z(r,@) be- 
comes a function z(x,y) after a change from polar to Cartesian coordinates. 


Thus the independent variables change from r and @ to x and y, while the 
dependent variable remains as z. 


O2w a) ( Ow | iz al Ow | 


(2.7) 


Example 2.28 
Oz 


O 
Given z = rsin 26, show that Pct + oa = &. 


Ox Oy 


ela 
Note that “a and oa are 


not simply the inverses of 


Ox Ox 
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Solution 

We have 
02 - @C20r-. 22.00 sin 0 
wae + Aan sin 20 cos 6 + 2r cos 26 (- 


= sin 20 cos 8 — 2cos 26 sin 6 


and 
O O 6 6 
= ee 5 Be + i 5, = sin 20 sin 6 + 2r cos 20 —— 
= sin 20 sin 6 + 2 cos 26 cos 0, 
therefore 
Oz Oz . 
t— +y— =rcos6 (sin 26 cos # — 2 cos 26 sin 6) 
Ox Oy 


+ rsin 0 (sin 26 sin 6 + 2 cos 26 cos 6) 
= rsin20(cos*@+sin*0)=z. I 


Exercise 2.24 


Pw , 
Given a function w(z,y), find an expression for ae terms of r and 6, where 
J Ow Ow 
x =rcos@ and y=rsin@. Hence find an expression for —> + —— in terms of 


2 2 
partial derivatives with respect to r and 0. Ox Oy 


Laplace’s equation in polar coordinates 


The operator 
a ome a 
Ox2 * Oy? ” a 
is known as the Laplacian operator; you will meet it again in Block II. 


The partial differential equation 

O7w = O7w _ O7w > 

Ox? Oy? Oz? 
is the three-dimensional Cartesian form of Laplace’s equation, which 
arises in many areas of applied mathematics. For various reasons, it is 
often useful to transform this equation into spherical or cylindrical polar 


coordinates, and to do that we require methods similar to those developed 
above. Applying the results of Exercise 2.24, the two-dimensional form of 
ee 
aplace’s equation —~ + — > = 0 becomes 
- ‘ Oz? Oy/* 

Ow 107w 10dw = 

Or2 72 962 or Or 
(This form of Laplace’s equation may appear to be more complicated than 
the original. In fact, as we shall see later when we discuss separation of 
variables in Block I, this form can be advantageous.) 


Exercise 2.25 


Transform the partial differential equation 


Ow Ow 
at + oe == 0), 


where a and b are constants, by making the change of coordinates u = bx — ay and 
v = ax + by, so the new independent variables are u and v. 
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Exercise 2.26 


Transform the partial differential equation 

Poy Cw 

Oxr2 Oy? 
by making the change of coordinates u= x+y and v = x — y, so the new indepen- 
dent variables are u and v. 


End-of-section Exercises 


Exercise 2.27 


Find the cosine of the angle 9 between the vectors a = 2i—j+ 3k and b =1+4 2j—k. 


Exercise 2.28 
Given a = i+ 2j+ 3k, b = 21+ 2] —k and c=i—j-+ 2k, calculate 
(a-c)b—(a-Dbjc. 


Exercise 2.29 


(a) Show that the curve y? = 4ax, where a is a constant, can be defined in param- 
etrised form by writing x = at? and y = 2at. Find the length of the curve 
from the origin to the point (a,2a). [Hint: Express the length as an integral 
in terms of t, then use integration by parts. 


(b) Show that the surface 


oe y" 


4 9 


can be specified in terms of two parameters u and v in the form 7 = 2(u+ v), 
y = 3(u—v), z = 4uv, and show that the point P with coordinates x = 2, 
y = 3, z = 0 lies on the surface and corresponds to choosing u = 1 and v = 0. 
Find the equation of the tangent plane to the surface at P, and find a vector 
perpendicular to this plane. 


Exercise 2.30 
(a) Find the equation of the tangent plane to the surface 
z= f(x,y) = (a+ 2y+ 1)(@—y + 2) 


at the point (0,0,2) on the surface. Show that the point corresponding to a 
vector of the form 


r = 2k + A(i+ 3k) + pj + 3k), 
where 4 and p are arbitrary scalars, lies on this tangent plane. 


(b) Given o(2,y,z) = z— f(x,y), evaluate V¢ at the point (0,0, 2). Show that any 
vector r — 2k (which is parallel to the tangent plane at (0,0, 2)) is perpendicular 
to Vé¢ evaluated at the point (0,0, 2). 


Exercise 2.31 


(a) Calculate the gradient of the scalar field F(x, y, z) = x? — 3xy+ y® + 2az + 2 
at the point (1, —1, 2). 


(b) Calculate the divergence of the vector field F' = sinh(xy)i+ cosh(yz)j + e**k. 


2.5 Outcomes 


Exercise 2.32 


Find the work done by the force F = (x + yz)i+ (y+ 2z)j + (z+ cy)k in moving 
a particle from the origin to the point P with coordinates (1,1, 1): 


(a) 
(b) 


along the straight line joining the origin to P; 


along the curve z = t, y= ¢*',2=?°. 


Exercise 2.33 


(a) 


Bi 


The function ¢ is defined by ¢(x, y, z) = x*yz, and the vector field F' is given 
by F = Vo. A point P with coordinates (cost, sint, sin 2t), 0 <t < 27, is an 
arbitrary point on a curve C. Show that the point P, corresponding to t = 0 
coincides with the point P» corresponding to t = 27, and show that 


[Fo ar =o 


d 
Hine Calculate > (cos* t sin t sin 2t). 


Generalise the result of part (a) by considering an arbitrary vector field F’ = 
Vo for some function ¢(x,y,z), where P, with coordinates (x(t), y(t), z(t), 
t; <t< te, is an arbitrary point on a curve C. If the point P, corresponding 
to t; coincides with the point P2 corresponding to tz, show that 


[ Fo) a L¢ 


5 Outcomes 


After studying this chapter you should be able to: 


understand the meaning of continuity for functions of many variables; 
understand the meaning of and be able to calculate partial derivatives; 
determine the tangent plane to a surface at a point; 


calculate first- and second-order Taylor approximations for functions of 
several variables; 


calculate and classify stationary points of functions of two variables; 
apply the chain rule for functions of several variables; 

integrate over two-dimensional planar regions in Cartesian and polar 
coordinates; 

integrate over non-planar regions when given an appropriate coordinate 
system; 

understand the concepts of scalar and vector fields; 

calculate the gradient of a scalar field and the divergence of a vector 


field; 


calculate the length of a curve specified either by a function or in para- 
metric form; 


calculate line integrals of vector fields; 

calculate the integral of a vector field over a surface to determine the 
flux; 

verify the Gauss divergence theorem in elementary cases; 

use the chain rule to find derivatives with respect to new coordinates; 
appreciate the derivation of Laplace’s equation in polar coordinates. 
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You may assume that the 
curve C' is smooth. 
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Solutions to Exercises in Chapter 2 


Solution 2.1 


The equation y = x? is recognisable as the equation of a parabola with the y-axis 


as its axis of symmetry. The equation x = y” also represents a parabola, but with 
the x-axis as its axis of symmetry. The equation x = a+ y? represents a translation 
of the latter parabola along the x-axis. These observations allow us to sketch the 
following curves. 


- 
ze ." A=15 
A=3 A=2.5 
2 A=1 
YA Seal a= 9.5 
Sy 2 = Se ‘ 
a 7 
=? 4 ~ | A 
/ } A= 2s . £2215, 
Figure 2.29 The curves 3-22 +y?2=2 Figure 2.30 The curves y/x? = 2 


The level curves in the first diagram are typical of a ‘well-behaved’ function, while 
the fact that the curves meet at the origin in the second figure means that the 
function is not defined there. 


Solution 2.2 
OF 
Ox 


OF 
= 2ry® + e¥ cos(xe¥) and apa 327y? + xe¥ cos(re”). 
Y 


Solution 2.3 
We have 


‘3 pita die [_see2)) . (1 seen), 


a 


and differentiating partially with respect to a gives 


d . d (1 Pe 
aa (| sin(ax) ir) ee (Za - cosa) se _ — ae 


But, using Leibniz’s rule, 


+ ( / Bex) ir) i / | = (sin(az)) dr = / panes 


therefore 


sin a 1 — cosa 


1 
/ 7 eon( os) da = — 
0 


a a? 


Solution 2.4 
OF 


OF 
(a) 3 ye*Y + 2x and = = ge™¥ + cosy. 
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O 
(b) Differentiating partially with respect to x gives ky l/ a4 Sz t/ a =; 86 
x 


Oz = Jz om Oz = Vz 
a — (x > 0). Similarly, Dy = ae (y > 0). 


(c) Ifcos(xyz) = 1, then xyz = 2km for some integer k, so z = 2k7/(axy). It follows 
that a = —2kn/(x*y) and 5 = —2kr/ (xy). 


(d) (i) Putting « = 1 and y = 2 gives z = —1, thus the point (1,2, —1) does lie 
on the surface. The partial derivatives are 


a 

5 = (y—2e +1) — 2(@ — 2y + 1) = by — 40-1 
and 

Oz 

ee aes 

Oy ei 7 


Evaluating the partial derivatives when x = 1 and y = 2, we obtain 
F,,(1,2) = 5 and F,(1,2) = —4. Thus the equation of the tangent plane is 


z= 5(a@ — 1) —4(y— 2) - 1. 


(ii) Putting « = 1 and y = 1 gives z = 1, so the point (1, 1,1) does lie on the 
surface. At this point, the partial derivatives are F,(1,1) = F,(1,1) =0, 
from which we deduce that this is a stationary point. At (1,1), the 
equation of the tangent plane reduces to z = 1, so it is parallel to the 


(x, y)-plane. 
Putting s = x —1 and t = y — 1, we obtain 
F(x,y) — F(1,1) = (s — 2t)(t — 2s) 
=e 40626 
= —2(4 — 1)? + 5(x — 1)(y— 1) — 2(y — 1)*. 
You may have noticed that this expression contains no linear terms in 


powers of x — 1 and y — 1; this is a consequence of the fact that there is 
a stationary point at (1,1). 


(iii) Putting y = 2— 2 gives F(2,y) — F(1,1) = —9(x — 1)? < 0, and putting 
y = x gives F(x, y) — F(1,1) = (a —-1)? > 0. 


We conclude that every neighbourhood of the point (1, 1) contains points 
at which F(x, y) — F(1,1) is positive and points at which it is negative. 
Thus the stationary point at (1,1) cannot be an extremum. 

(ce) We have 
OF 


OF 
ae = (w+ y )e® + Que” = (2 +1) +y° — le and ae 


so stationary points occur at the intersections of the circle (x4 +1)? +y? =1 


with the line y = 0, namely at (0,0) and (—2,0). 


Since F(x, y) — F(0,0) = (x? + y?)e” > 0, the stationary point at the origin is 
a local minimum. 


Solution 2.5 


es x ies x See oe — x _ 
(a) w, = ye ae ap, = PE , BO a = ye", Wey SE Ee oe oo, = 


x7e*¥. The quadratic approximation is therefore 


w(x, y) ~ w(1, 2) + wz(1, 2)(@ — 1) + w,(1, 2)(y — 2) 
+ 5 (Wee(1, 2)(x — 1)? + 2wey(1, 2)(a — 1)(y — 2) 
+ Wyyl1, 2g — ay") 
=e? (1+2(@—1) + (y—2) +2(@ — 1)? + 3(@ — 1)(y — 2) 
+ 4(y—2)?). 
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Shortly we shall see a 
simple criterion that will 
allow us to reach this 
conclusion more easily. 


You can find these by 

substituting y = 0 into the 
equation (x +1)? +y? =1 
to obtain x = 0 or x = —2. 
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Ow Ow ; 
(b) Since os cos x cosh y and ie sin xsinhy, the stationary points occur at 
Y 


the solutions of the pair of equations cos x cosh y = 0 and sinz sinh y = 0. Since 
cosh y > 1, the first of these equations gives cos xz = 0, so 


ait + & 
r=sZtnr= 5 mq for n an integer. 


For these values of x we know that |sinz| = 1, so the second equation gives 


= 0. 

The second partial derivatives are 
07 w h Ow ah 0? h 
—— =—sinxcoshy, z~ —=cosxsinhy, => =sinzcoshy. 
Ox? 7 Andy 4 Gy? " 


Using the previous notation, at the point (5 + n7,0) we have 
PR — Q? = —(sinxcoshy)” — (coszsinh y)? = —1, 


so the stationary points are not extrema (but saddle points). 


Solution 2.6 
OF ac 9 O° F 
If F(x, y) = f(z) + gly), then iF (y), So anda. 0 


(This simple observation will be of interest later, when we discuss the general 
solution of partial differential equations such as F,, = 0.) 


Solution 2.7 


Putting x = 0 (i.e. 2] = 22 = 23 = 0), we obtain A = F(O). Then differentiating 
partially with respect to each of the variables in turn, and putting x = 0, we have 


OF OF OF 
B=|— By =|— B3 =| — 
: (i) : (oe i ae) 


In order to find the remaining coefficients, we differentiate partially twice and then 
put x = 0, giving 


1 (06°F 1 (06°F 1 (06°F 
see a) ce’ (<3 a (3 


0° F 0° F 0° F 
Cio = | -——_ Cos = a 
- (5a ) any “ (Sa soda Cai (soa ) te 


Solution 2.8 


es, so 


(a) We have w =e 


: (cor 2s Tee. 

Using the form w = e*”, we have 
Ow 2. Ow dx dy 
—-=ye"y", —=zre"¥, —=cost, — =-sint. 
pel = Oy dt dt 


Hence, using the chain rule, 


dw Ow dx Ow dy 


—— = — — + — = ye*4 cost — xe*’ sint 
i fee” eee ce 


age (cos? = sin? saat 
which agrees with our first calculation. 


(b) We have 
Ox Oe — ; Oy Oy 


Ou eS a ay 


Solutions to Exercises in Chapter 2 


Then, using the chain rule, 
Ow OwOdx Ow Oy Ow Ow 
du Ou du Oyu “dx dy 
and 
Ow Ow Ox Pci Ow Oy Ow Ow 
Ov Oxdv. Oydv Oz Oy 


These relationships are valid for any function w, and it is convenient to write 


al =a—-+c— and z : “ 
du O O Ov Ox Oy 
Then 
oe hase 3 ae as se ne 
Ou2 -\ Ox = Oy Ox Oy)  —- Ox? Oxrdy : Oy?” 
O7w O O Ow Ow Ow Ow 
= ( hy 5 ) OSE +45) baa + (bc + ad)a f Oa y2 
0? w 0 6) Ow Ow > O*w 0? w Wem 
Solution 2.9 


(b) 


Figure 2.31 The region S$ divided into (a) vertical strips and (b) horizontal strips 


There are two alternative approaches. We could divide the region S into two regions, 
J and K, as shown in Figure 2.31(a), where: 


J is bounded by y = 0, x = 4 and y=9; 
K is bounded by xz = 4, y = 0, + = 8 and y = 16/z. 


We then have 


[eaaq [raat | a? dA, 
Ss J K 


where the areas J and K written as repeated integrals are 


x 4 
Sus A= ‘ (/ dy) ax = | x dx = [424], = 64 
J 0 0 
16/x 8 : 
ae ar? A= fos / yn éa = / l6rd7 = come = 384. 
K 0 4 


Thus 


a? dA = 64 + 384 = 448. 


WD 
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Alternatively, we could divide the region S' into two regions L and M, as shown in 
Figure 2.31(b), where: 


L is bounded by y = 2, y = 16/z and y = g; 
M is bounded by y = 0, « = 8, y= 2 and pea. 


We then have 


[ear= fears | x? dA, 
Ss rs M 


where the areas L and M written as repeated integrals are 


4 16/y 4 
/ z’* dA = / / x? dx | dy = / [Sex| ids dy 
L 2 . 2 ? 


4 3 
4096 y 
ae — — | dy = 108 
I ee 3 


Thus 
/ z* dA = 108 + 340 = 448, 
Ss 


as above. 


Solution 2.10 
We have 


nm /2 | m/2 
[Parao= f (/ r? dr w= / do = 7/6. 
S 0) O 0 


Solution 2.11 


We have 
es x 3 Beamon. 9? 
a aye ar ee 
Similarly (and taking advantage of the symmetry), 
Ow yy? — 2? 
Oy? e (a? + y2)5/2° 
Solution 2.12 
We have 
oF ae os OF Qry 
Ox (a? + y?)? oy (ie tee 
Then 


OF OF 2Qy(3x7 — y*) 


Solution 2.13 


O O 
(a) If w =e* cosy, then Fs = e* cosy and = —e* sin y. 
x y 
0? o 0? 0? 
So — =e” cosy and ae = —e* cosy, which gives — + a == {). 
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0? 0? 0? 0? 
So as = 62 and =a = —6z, which gives — + a ==). 
Solution 2.14 
We have 
F, = —22+20e¢+ 2y and 2, = -—2> 27 > 109. 


The single stationary point occurs at the intersection of the lines 10” + y = 11 and 
xz + 5y = 6, and thus is at (1,1). Since F,, = 20, Fyy = 10 and Fz, = 2, it follows 
that (using the standard notation) 


PR-— Q* = 200 —4= 196 > 0, 


therefore the stationary point is a local minimum. 


Solution 2.15 


Let P, with coordinates (x,y,z), be an arbitrary point lying in the given plane. 
Then the distance from P to the origin is d = \/x?2 + y2+ 22. The algebra is 
simpler if we minimise w = d? = x? + y? + 2? rather than d, and we note that on 
the plane we have y = 2x + 2z — 16, so w = x7 + (2a + 2z — 16)? + z?. As P ranges 
over the plane, there will be a point which is closest to the origin, and at this point 
the function w(z, z) will have a minimum. At this stationary point, 


O 
a = 27 + 4(2e + 2z — 16) = 102 + 8z — 64 =0 
ie 
and 
Ow 
ae = 4(27 + 2z — 16) + 2z = 82 4+ 10z — 64 = 0, 
z 
therefore 7 = z = 32/9, so y = —16/9. Since there is only one stationary point, we 


know that this must be the required minimum, and it occurs at (32/9, —16/9, 32/9). 


We can confirm this conclusion by performing the evaluation 


0?w 07 w ( O07 w 


2 
———- — | -—-] =10x 10-8* = 36 > 0, 
on" OF" ae) 
2 
Ww 
so the stationary point is an extremum, and since 22 = 10 > 0 this extremum is 
fy 


a minimum. 


Solution 2.16 
We have 


Ow Ow Ow 3ryZz 


OxO0yOz . OzOxOYy = OyOzOxr (x2 + y? + 22)8/2° 


Solution 2.17 
We have 


7 cos 6 ae a 
/ (/ psn w= | [50° sin 6°" " do = | cos? @sin 6 dé. 
O O 0 4 


Now we can either make the change of variable u = cos@, or note that sin@d@ = 
—d(cos 6): 


us cos T 
/ (/ psn de = -3/ cos* 6. d(cos 0) = 4 [— cos” a0 a 
0 0 0 


wl 


107 


Note that here we have 

chosen to eliminate y, but 
we could equally well have 
chosen to eliminate x or z. 
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Solution 2.18 


(b) 


Figure 2.32 


Dividing the region into horizontal strips, then summing over the strips (see Fig- 
ure 2.32(a)), we obtain 


If. exp(x?) dA = [ (/, exp(a?) de dy, 


which cannot be evaluated directly in terms of the elementary functions. 


Alternatively, dividing the region into vertical strips, then summing over the strips 
(see Figure 2.32(b)), we obtain 


// exp(2*) dA = [ fe exp(x”) iy aa = [ ly exp(a?)]°/° dx 


= [ [4x exp(a”)] dx 


é [exp(2”)], 


= sie = 1), 


The point to notice here is that the order in which we perform the integration may 
well affect the difficulty of the calculation. 


Solution 2.19 


(a) We have 
ja} = V1+44+1=V6 and |b| = V4+1+1= V6, 
and 
avb= (1,2,-1)-(2,4,1) =e 
So 
-b 
cea NB 2 


(b) The magnitude of the vector c is V3? + 2? + 5? = V8. 


Solution 2.20 
We have 
F, = 2x sin|(x + y)r] + (a? + 2y?)mcos|(x + y)z] 
and 
F, = 4ysin|[(x + y)m] + (x? + 2y?)mrcos[(x + y)z]. 
It follows that F,(1,1) = F,(1,1) = 3a, therefore VF = 3z(i + j). 
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Solution 2.2 | 


(a) We have 
V-.F= oF of La g? + y? + 27 3 
Ox Oy Oz (a2 + y2 + 27)8/2  (a?%@ + y? + 22)1/2 
= 2 
(a2 4 y2 + 22)1/2° 
(b) We have 
A alee oh geet ges OF k = 227i + 2yj + 2zk, 
Oy Oz 


therefore V - ii = 6, 


Solution 2.22 


ao 
(a) —=2ta = 3t7, so the length of the curve is given by 


el @ le [v@ (2t)2 + (302)? dt 
HS @ 


[ t(4 + 9t?)1/? dt 
0 


# |(4+ of?)9/2] 
& (37/37 — 1). 


(b) We have a’(t) = 1, y'(t) = V2t, 2/(t) =t?, so the distance travelled is 


2 
r= | V1+ 2t? + t* dt 
0 


= [area 


2 
= |é + ae = a 


(c) The straight line joining the origin to the point (1,2,3) can be parametrised 
asx=t, y = 2t, z= 3tfroomt=0tot=1. Thenr = (i+ 2j+ 3k)t and 


1 
/ Fir)\.de= / (2t3i + 3tj + 4t°k) - (i+ 2j + 3k) dt 
e 0 


1 
z= / (6t + 14t°) dt = 65. 


Solution 2.23 
First we integrate V - F over the region R. Since 


O(x sin(7y)) n O(zz) 
Y 


>. 7s 
Ox 


= sin(7y), 


we have 


[verw=[ OU sin(y) )dr) ay) a 
Ene) 
-{? = dz = 


Next we integrate the field F' over the surface of R. Since R is a cube, we need to 
evaluate six integrals, one for each of the faces at x = 0, x = 1, y=0, y=1,z=0, 
pm 4. 
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For the face at x = 0, we have dA = —idy dz, so F(0,y, z)-dA = 0, and we obtain 


= F(0,y,z)-dA =0. 
xt=0 
For the face at x = 1, we have dA = idydz, so F(1,y, z)- dA = sin(my) dy dz, and 
we obtain 


1 1 
a = F(l,y,z)-dA= / (| sin(7y) iy) a = = 
et 0 0 sa 


For the face at y = 0, we have dA = —jdzdz, so F(xz,0,z)-dA = —xzdzrdz, and 
we obtain 


1 1 
S= | F(x,0,2)-dA =~ [ (| v2 dr dz = —4. 


For the face at y = 1, we have dA = jdzdz, so F(z,1,z)-dA = xz dzdz, and we 
obtain 


1 1 
si = | F(x.1,2)-dA= | (| vedr) de =} 
gel @) O 


There is no contribution to the surface integral from the faces z = 0 and z = 1, 
since in each case F' - dA = 0. 


Adding the surface integrals gives 


2 
S1+ 52+ 53+ 54=—, 


the same result as obtained by finding the volume integral. 


Solution 2.24 


Multivariable and vector calculus 


We have 
O7w w 
Dit By ay 
4 05. (sino + 1) tm noe eS) 
r o8 r OO Or r 06 
ce 2sin@cosé 0?w cos’ 9 O"w | cos! Gdw  2sin@cosé Ow 
Or? r Oroe r2 667 OOF r2 aly 


Hence, using equation (2.7), 
Ow , Ow ah ‘a2 1 rw 1 Ow 
Or2 Oy? Or? 2 OG? Or” 


Solution 2.25 


From the chain rule we have 


Ow a Ow Ou Ow Ov Ow ee 


icine * aM ae sti 


oO jee ee Ju — 


and 
dw _ dwdu , dwdv dw , dw 
Oy Oudy dv dy Ou Ov’ 
therefore 
Oe -.O0 — Ow Ow Gy Ow). Jy i. 0@ 
a Gia a (oe +a ae) +o an +on ) =a +06 a. 
The given partial differential equation therefore reduces to 
O 
— =0. 
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Solution 2.26 
We have 
Ou Ou | Oe Ov 
ere “Wie "ae OS 
It follows from the chain rule that 
Ow Owdu  Owdv Ow Ow 
dx «Oudz Ovdx Bu dv 
and 
Ow  Owdu  Owdv_ Ow Ow 
dy Oudy Ovdy du dv 
These relations are valid for any function w, and it is convenient to write 
QO oO O 
dx du’ dv 
and 
- - # O 
dy Ou dv’ 


in which case we have 


oe=(2+2)(F+2)-F Ow 
Ox? Ou Ov Ou Ov Ou? Oudv 
and 
Ow a- 2 Ow Ow O7w 07 w 
oy (= - 5) (se -S) ~ Ou2  ~ dudv 
It follows that 
aw Bw _ Pw 
Og? . Oy Oudv’ 
so the given partial differential equation reduces to 
O7w 
dudv 


Solution 2.27 


We have 
cos § = 2 = _ Qi-j+3k)-G+2j—-k) 
ja] |B| 22 + (—1)?2 + 32,/12 + 22 + (-1)? 
_ 2-3-2 
V14V6 
2 
~ Sea 
Solution 2.28 
We have 


a-c= (i+ 2j+3k)-(i-j+2k) =5 
and 
a-b= (i+ 2j + 3k) - (2i+ 2j —k) =3, 


SO 


07 w 


Ov? 


O07 w 


Ov?” 


(a-c)b— (a- b)c = 5b — 3c = 5(2i + 2j — k) — 3(-—-j + 2k) 


= 7i + 13] — 11k: 
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Solution 2.29 
2 
(a) Since ¢ = | we have « = at? =a (=) = = so y* = 4az. 
a a 


Since x = at” and y = 2at, we have 


Thus the length of the curve is 
1 1 
Ss = / VJ (2at)? + (2a)? dt = 2a f V1+t? dt. 
0 0 
1 
Putte f= / V1+t? dt and integrating by parts, we have 
0 
1 t (7? 
ik i+#| - | sat 
0 0 V1+#? 
1 2 1 
| ee: 1 
=vi- | a+ | at 
0 V1+#?? 0 VI+¢ 
. 1 
1 
=vi-1+ | sat 
o v1+#? 
= V2-—I+ [sinh t]- 


0°? 
SO 


S =2al =a(v2+In(1 + v2), 


(b) Just as a curve can be parametrised by a single parameter as in part (a), so a 
surface can be parametrised by two parameters, and here we have x = 2(u+ v), 


y = 3(u—v) and z = 4uv. Solving the first two of these equations gives 


cy co Y 
2u=-+= d 2v=—-—-— =. 
a ae 
Then substituting into the third equation gives 
oa x2 oy? 
= es 


showing that the parametrisation in terms of u, v and w does specify the given 


surface. 


Putting w=] and v=, weobtain z= 2, y= 3 and — =70,.36: therpomt 


(2,3,0) does lie on the surface. 


gy? 
Differentiating z = 23 ee partially gives 
Oz 2 aa Oz iy 
or----2 oo 9’ 


80. 23(2;3) = | and 230,37) = —¢. The equation of the tangent plane at the 


point (2, 3,0) is therefore 
z= (2-2) — $y -3), 


which can be written as 3x — 2y — 3z = 0. 


There are several ways of finding a vector perpendicular to the tangent plane, 
but the following is probably the simplest. First we notice that the plane 
3x — 2y — 3z = 0 passes through the origin, and if we let (x, y, z) be any vector 


from the origin to a point in this plane, we can see that (3, —2, —3) - (x, y, z) 


3x — 2y — 3z =0. So the vectors (x,y, z) and (3, —2, —3) are orthogonal, thus 


the vector 3i — 2j — 3k is perpendicular to the tangent plane at (2, 3,0). 


It is also worth noting that the equation of the plane may be written as 
f(x,y, z) = 3x — 2y — 3z = 0; and then the vector Vf = 3i — 2j — 3k is per- 
pendicular to this plane. This is because Vf is always perpendicular to the 


(tangent plane to the) level curves of f. 
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Solution 2.30 


(a) We have 


a 
SF (r+ 2y +1) + (o-y+2) = 20+ y 43 


5 et ee 


The equation of the tangent plane at (0,0, 2) is therefore 
z= fx(0,0)a + fy(0,0)y + f(0,0) = 3x + 3y + 2. 


The point corresponding to the vector r = 2k + \(i+ 3k) + w(j + 3k) has co- 
ordinates (A, 4,2 + 3+ 3), and these satisfy the equation z = 3x + 3y + 2, 
so the point does lie in the plane. 
(b) We have 
Og, _ OG, , Ob a os. 


Sime to Sek aa: Se 
Vo et By) Dz . Dy? j+k, 


whose value at (0,0,2) is —3i — 3j+k. 
It follows that 
Vo: (r — 2k) = (—3i — 33 +k) - Ai+ nj +3(A 4+ w)k) = 0, 
so the vectors are perpendicular. 
In order to clarify what is happening here, note that the level curve ¢(z, y, z) = 
z — f(x,y) = 0 corresponds to the surface z = f(x,y). So at a given point on 
the surface, vectors which lie in the tangent plane to z = f(x,y) will always 
be perpendicular to V@. 
Solution 2.31 
(a) We have 
F, =22—3y+2z, F,=—3e+3y*, F,=2r+1, 
SO 
VF = (2x — 3y + 2z)i+ (—3a + 3y?)j + (2x + 1k. 
At the point (1, —1,2), we obtain VF = 91+ 3k. 


(b) We have 
Og Ay nS 
= ycosh(ry) + zsinh(yz) + xe** 
Solution 2.32 


(a) We can parametrise the straight line joining the origin to the point P by writing 
g=t, y=tand z =T Wiest). Then 


F(r) =(a+yz)it (yt xz)jt+(z+cry)k 
=(t+?)i+ (¢+)j+(t+t’)k 
= (t+ ¢)i+j+k) 


dz a 
re a ae +[k=i+j+k, 


so the work done is 


F(r)-dr = f (t+ A) +i +h “it jrei a 


1 
=3/ (t+ 07) dt 318 [be dl = 5 
0 
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(b) The particle moves along the curve x = t, y = t?, z = t? from the origin to the 
point P as t increases from 0 to 1. On this curve, we have 


F(r) = (x@+ yz)it+ (y+ ez)jt+ (z+ay)k = (¢4+ P)i+ (# + t*)j + 2°k 
and 
aw daz ay . 


dz 2 
ae ait a +k = i+ 2t¢j + 3t*k, 


so the work done is 


1 
/ F(r)-dr = / ((t+¢°?)i+ (t* + t*)j + 2t°k) - (i+ 2tj + 3t7k) dt 
C 0 
1 
= / (t +t? + 2t(t? + t*) + 6t°) dt 
0 


1 
=f (t + 2t? + 9t°) dt = 3 
0 
Generally, the value of a line integral depends upon the path of integration; this is 
an example of one that does not. For the general vector field 
Fir) = F,(a2,y, z)i+ Fo(a, y, z)j+ F(a, y, z)k, 
it can be shown that the integral / F(r) -dr is independent of the path C between 


two given end points if there exists a function (x, y, z) such that 


_ 99 _ 9 _ Ae. 
a a oe ee 


in other words, if F(r) = V¢. In this example it is easy to verify that we can use 
b= s(x? +y? + 27) + zyz. 
Solution 2.33 


(a) We have 
(cos 0, sin 0, sin0) = (1,0,0) = (cos 27, sin 27, sin 47), 
so P,; and P are the same point. 
Now 
F=Vo=>= a + oo + ook = = Qeyzi + x*zj +2 yk, 


SO 


F(r)-dr = / (Qeyzi + x?zj + x*yk) - (dri + dyj + dzk) 
“i c 
27 
; ; aie? dz 
=| (Qar(t)y(t)z(t)i + x(t)? z(t)j + x(t)?y(t)k) - | —it+ i j+—k) d 
0 dt dt? dt 
27 
= / (2costsint sin 2ti + cos’ tsin 2tj + cos’ t sint k) : (— sinti+costj + 2 cos 2t k) dt 
0 
27 ' 
= / (—2 cost sin? t sin 2t + cos® t sin 2t + 2 cos” t sin t cos 2t) dt 
0 


2n 
= / = (cos* tsint sin 2t) dt 
9 «Coat 
== [cos* tsint sin 2t] - me 
(b) Using the chain rule, we have 
o, Q, an dx oy, a 
Fi(r)-dr —i+—~—j+ — —ji+ —j k } dt 
[ (r) = | (i+ fi+ Ze) oe a 


rat Oddx Oddy Odddz 
=f (Sete eta) 


-f(Ban(wole=a anatraaty 
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Laplacian operator 99 
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linear approximation 63, 64 
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linear differential equation 31 
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local extrema 64 
local maximum 64 
local minimum 64 
logarithm function 10 
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maximum 64 

method of undetermined coefficients 
minimum 64 

mixed derivative theorem 66, 79 
modulus of a complex number 22, 24 
multiple integral 73, 78 


nabla 8&5 


of 
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natural logarithm 10 roots of a complex number 25 
non-linear differential equation 31, 54 
normal vector 91 saddle point 65 

scalar 80 
O notation 8 scalar field 82, 85 
odd function 5 scalar product 81, 82 
operator 6 second-order partial derivative 66 
order notation 8 separable equation 32 
order of a differential equation 31 series 7 
order of magnitude 8 sine function 10, 25 
ordinary differential equation 29 slope of a surface 83, 84 
orthogonal functions 55 smooth curve 87 


smooth function 6, 61 
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surface area 79 
volume 14 
stationary point 64, 68, 71, 72 
sum of a series 7 
surface area 
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partial derivative 61, 66 

first-order 61 

higher-order 65 

second-order 66 
partial differential equation 29 
partial fractions 18, 44 
partial sum 7 
particular integral 31, 37 
particular solution 32 
polar coordinates 24, 76, 77, 79, 96, 99 
polar form of a complex number 24 
polynomial 7 
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projection 81 


tangent function 11 

tangent plane 61, 63, 64 

Taylor approximation of order two 67, 70 
Taylor series 9, 10, 67 

transform 6 

triangle inequality 27 

trigonometric functions 10-12 


quadratic approximation 66 triple integral 78 


quadratic factor 20 
uniform wire 14 


radius of convergence 9 unit vector 81, 91 

rate of convergence 8 

rational function 7, 18 vector 80 

real function 5 | vector differential operator 85 
real line 5, 22 vector field 86 

real numbers _ 5, 22 volume integral 78 

real part of a complex number 22 volume of asphere 14 


remainder term 8 
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